Next Article in Journal
Splicing Modulators Are Involved in Human Polyglutamine Diversification via Protein Complexes Shuttling between Nucleus and Cytoplasm
Next Article in Special Issue
Liposome-Based Carriers for CRISPR Genome Editing
Previous Article in Journal
Visualization of Cellular Membranes in 2D and 3D Conditions Using a New Fluorescent Dithienothiophene S,S-Dioxide Derivative
Previous Article in Special Issue
Expected and Unexpected Effects of Pharmacological Antioxidants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Screening and Analysis of Possible Drugs Binding to PDGFRα: A Molecular Modeling Study

by
Matteo Mozzicafreddo
1,*,
Devis Benfaremo
1,2,
Chiara Paolini
1,
Silvia Agarbati
1,
Silvia Svegliati Baroni
1,2 and
Gianluca Moroncini
1,2
1
Department of Clinical and Molecular Sciences, Marche Polytechnic University, 60126 Ancona, Italy
2
Clinica Medica, Department of Internal Medicine, Marche University Hospital, 60126 Ancona, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(11), 9623; https://doi.org/10.3390/ijms24119623
Submission received: 28 February 2023 / Revised: 3 April 2023 / Accepted: 12 April 2023 / Published: 1 June 2023
(This article belongs to the Special Issue Recent Advances in Drug Discovery)

Abstract

:
The platelet-derived growth factor receptor (PDGFR) is a membrane tyrosine kinase receptor involved in several metabolic pathways, not only physiological but also pathological, as in tumor progression, immune-mediated diseases, and viral diseases. Considering this macromolecule as a druggable target for modulation/inhibition of these conditions, the aim of this work was to find new ligands or new information to design novel effective drugs. We performed an initial interaction screening with the human intracellular PDGFRα of about 7200 drugs and natural compounds contained in 5 independent databases/libraries implemented in the MTiOpenScreen web server. After the selection of 27 compounds, a structural analysis of the obtained complexes was performed. Three-dimensional quantitative structure–activity relationship (3D-QSAR) and absorption, distribution, metabolism, excretion, and toxicity (ADMET) analyses were also performed to understand the physicochemical properties of identified compounds to increase affinity and selectivity for PDGFRα. Among these 27 compounds, the drugs Bafetinib, Radotinib, Flumatinib, and Imatinib showed higher affinity for this tyrosine kinase receptor, lying in the nanomolar order, while the natural products included in this group, such as curcumin, luteolin, and epigallocatechin gallate (EGCG), showed sub-micromolar affinities. Although experimental studies are mandatory to fully understand the mechanisms behind PDGFRα inhibitors, the structural information obtained through this study could provide useful insight into the future development of more effective and targeted treatments for PDGFRα-related diseases, such as cancer and fibrosis.

1. Introduction

Platelet-derived growth factor receptors alpha and beta (PDGFRα and PDGFRβ) are membrane receptors that play a key role in a variety of diseases, such as cancer [1,2,3,4,5], immune-mediated pathologies, such as systemic sclerosis (SSc) [6,7,8,9,10], and viral infection [11], by exerting function on the transduction of extracellular signals into the cell. PDGFRs are transmembrane glycoprotein dimer molecules consisting of an extracellular ligand-binding region divided into five Ig-like domains bound to an intracellular receptor tyrosine kinase (RTK) domain through a single transmembrane alpha helix. The tyrosine kinase domain is a typical type III RTK formed by a bilobal kinase, a juxtamembrane (JM) domain, and an activation loop (A-loop) [11,12]. The JM domain is inserted into the active site, while the A-loop (that starts with a conserved amino acid sequence DFG) controls access to it and partially blocks the ATP and substrate site in the inactive state of the enzyme. In the active state of the RTK, the DFG motif is located in a hydrophobic pocket close to the ATP binding site. When an inhibitor binds this pocket, the ATP binding site undergoes a conformational change that reversibly or irreversibly prevents its normal function [13].
Considering that the human PDGFRα has become one of the most important therapeutic targets for the above-mentioned diseases [4,5,6,7,8,14,15,16,17], its modulation/inhibition is a critical aspect that determines the investigation of new ligands or new information to design new effective drugs. As previously reported [18,19,20,21,22,23], inhibiting PDGFRα (overexpressed and/or abnormally activated) can reduce many types of cancer, including prostate, ovarian, breast, pancreatic, and liver cancers, by inactivating downstream signaling pathways that regulate cell proliferation, migration, and angiogenesis [24,25,26]. On the other hand, these signal pathways can induce proliferation and chemotaxis of myofibroblasts, collagen production, and adhesion in endothelial cells [27]. In pathological states, which are also correlated to viral infections, they are upregulated in patients with a progressive fibrotic disease such as SSc [28,29,30].
For these reasons, the structure-based drug design (SBDD), which promotes the development of novel drugs with a potential affinity for therapeutic targets, has turned out to be an essential tool for rapid and cost-efficient lead discovery. SBDD is a growing, iterative, and powerful approach that involves the structural evaluation of targets in the drug discovery process. It has the ability to reduce the time and cost of developing new drug lead molecules with potential therapeutic effects [31]. A virtual ligands database screening, followed by molecular docking, molecular dynamics (MD), and three-dimensional quantitative structure–activity relationship (3D-QSAR) analysis, as previously reported [2,3,5,32], provides an excellent workflow (Figure 1) to assess an efficient SBDD.
Here, we analyzed 27 molecules selected from the results of structure-based virtual screening following this workflow and characterizing the most important features of an ideal ligand/inhibitor for the intracellular RTK domain of human PDGFRα.

2. Results

2.1. Virtual Screening of the Small Compounds Database

The structure-based virtual screening of the databases (as descripted in the Methods section) against the human intracellular PDGFRα (hiPDGFRα) provided an early version of the best compounds list (Tables S1–S5). This was evaluated by taking the best ligands of each of the 5 drug-like chemical libraries, resulting in the 27 ligands reported in Table 1. Among these libraries, the Drugs-lib had the most involvement, showing ligands with higher affinities, although the natural products selected showed Kd values in the sub-micromolar range, indicating medium–high affinities.

2.2. Molecular Docking

The molecular docking analysis, which was used to check the virtual screening results, provided affinity values reported in Table 1 and binding geometries of the ligand/hiPDGFRα complexes. These values are in the nanomolar/sub-micromolar range and confirm the ranking obtained by the structure-based virtual screening (Figure S1).
To validate the docking process, we compared the Imatinib/hiPDGFRα predicted complex with the available crystal structure, obtaining a ligand–ligand superimpose root mean square deviation (RMSD) value of 0.198 Å. This value is significantly lower than 1 Å, indicating a high reliability of the method.
All complexes involved the hiPDGFRα active site, although the docking box was set for the entire protein, indicating no other binding site. The four best ligands interacted with hiPDGFRα through the formation of five hydrogen bonds with Glu644, Thr674, Cys677, and Asp836, showing very similar geometries. As reported in Figure 2A, Bafetinib also interacted with the protein, forming six hydrophobic interactions with Val607, Glu644, Ile647, Thr674, Leu825, and Phe837, and one π-cation interaction with Lys627. This π-cation interaction was also formed in the Flumatinib/hiPDGFRα complex, which additionally showed nine hydrophobic interactions (with Leu599, Val607, Lys627, Glu644, Val658, Thr674, Leu825, Asp836, and Phe837) and a halogen bond between Glu644 and the trifluoromethyl group of the ligand (Figure 2B). This π-cation interaction, the halogen bond, and the distances and number of hydrophobic interactions between the ligand and the hiPDGFRα seem to be crucial for binding affinity.
Although the Imatinib/hiPDGFRα crystal structure shows seven hydrogen bonds (with Glu644, Met648, Thr674, Cys677, His815, Val816), a salt bridge (Asp836), two π-cation interactions (with His816 and Tyr676), and seven hydrophobic interactions (with Leu599, Lys627, Glu644, Val658, Ile672, Thr674, and Asp836), the predicted Imatinib/hiPDGFRα complex, after the molecular dynamics analysis, shows only five hydrogen bonds previously described and nine hydrophobic interactions (with Leu599, Val607, Lys627, Ile672, Thr674, Tyr676, Leu825, Asp836, and Phe837). This could be explained by the difference in temperature between the two types of tests (293 K for the crystallographic analysis and 300 K for the molecular dynamics).
Regarding the binding between natural ligands such as curcumin and quercetin and hiPDGFRα, this involves only a minor number of hydrogen bonds and hydrophobic interactions, even if their strength and number establish a medium–high affinity (Kd,pred = 0.2 μM and Kd,pred = 0.7 μM, respectively). This still determines an excellent competition against ATP that has an affinity for the wildtype enzyme in the order of sub-millimolar (Km = 179.6 μM) [12], underscoring the potency of their enzymatic inhibition capabilities.

2.3. Molecular Dynamics

The stability of the best 4 predicted complexes, obtained through molecular docking, was analyzed using molecular dynamics simulation for 50 ns. The first four nanoseconds of these simulations are represented in Figure 3 since no further RMSD variations were established. The calculated RMSD of each ligand and complex trajectory with respect to the backbone indicates that the complexes achieve sufficient stability (RMSD < 0.17 nm, 1.7 Å) within less than 4 ns, except for Flumatinib (Figure 3G), which shows RMSD values of about 0.22 nm.
The RMSD trend of Bafetinib (Figure 3A) and Flumatinib (Figure 3G) with respect to the other ligands seems to require a few more picoseconds and nanoseconds, respectively, to reach final stability. This behavior could be explained by the establishment of the π-cation interaction.
The average short-range Coulombic and Lennard–Jones interaction energies calculated during the molecular dynamics analysis are reported in Table 2. As noted in the molecular docking section, the electrostatic contributions between Bafetinib and Flumatinib are essentially equivalent, while there is a significant difference in terms of intermolecular pair Lennard–Jones potential, indicating a more efficient binding of Bafetinib to the active site.

2.4. Three-Dimensional Quantitative Structure–Activity Relationship (3D-QSAR) and ADMET Analyses

The three-dimensional QSAR analysis produced comparative molecular fields analysis (CoMFA) models considering the CoMFA potential or field, the best of which are summarized in Table 3 based on optimal principal components (OPC). Among these, MW, nHA, nHD, LogP, TPSA, rotatable bonds (RBs), molar refractivity (MR), length, and max length were taken into account and modulated to achieve the best model. Some of them are reported in Table 4 and Table 5.
Figure 4 shows the relationships between experimental and calculated activities by CoMFA models. The fitting statistical results, also reported in Table 3, indicate that the steric model (STE) is the more consistent and robust one, even if the electrostatic one (ELE) shows a relatively good squared correlation coefficient r2 and internal predictive coefficient q2. This would confirm the result of interaction energies obtained by Gromacs, which describes the best complex in terms of affinity as the one with the higher Lennard–Jones energy contribution, the forces involved in the steric effect.
ADMET predictions, in particular, some properties and relative descriptors of the 27 selected molecules, are reported in Table 4 and Table 5. A major part of these compounds satisfies the conditions regarding the acceptance of the first two indexes (Lipinski and Pfizer) related to absorption, permeability, and toxicity. Only Ditercalinium does not satisfy the conditions of these indexes since it presents the heaviest and longest molecule and the highest MR and LogP, which are out of range. Several compounds could cause skin sensitization and respiratory toxicity—although no compound would cause acute toxicity during oral administration.

3. Discussion

The PDGFRα is a critical factor that plays an essential role in regulating cell proliferation, survival, and chemotaxis [24,25,26,27]. Several studies have identified PDGFR as a potential target to treat pathologies, such as cancer and progressive fibrotic diseases [4,5,6,7,8,14,15,16,17], by searching for synthetic or natural compounds acting as ligands/inhibitors. Currently, however, no one has considered involving such large compound libraries for this purpose, such as those implemented in the MTiOpenScreen web server.
The details of binding interactions between human intracellular PDGFRα and the compounds selected after the structure-based virtual screening were revealed by molecular docking and the subsequent MD simulations. Molecular docking analysis, in addition to complexes’ affinities, was able to provide important critical issues regarding the geometries. The predicted affinities for the 27 selected complexes vary from the nanomolar order (1.01 nM for the Bafetinib/hiPDGFRα complex) to the micromolar one, a typical range of values for a good drug inhibition and reliable QSAR analysis. Some of the most stable complexes show a π-cation possible interaction with Lys627, halogen bonds with Glu644, and numerous hydrophobic interactions.
MD analysis reveals that the complexes achieve their stability in a very short time and remain stable for all the MD simulations, demonstrating no substantial unfolding of the structures. The type of interaction energies calculated during this analysis shows that the steric contribution (the Lennard–Jones energy) is predominant in the most stable complexes. This behavior is confirmed by CoMFA models obtained through the three-dimensional QSAR. The steric model field (or at least a combination of the two models) produces a very consistent activity predictor model—suggesting the descriptors included in this field and the structural information obtained by molecular docking analysis—as a basis for designing new effective drugs. These predicted models and structural information confirm the molecule ranking reported in Table 1 and identify Bafetinib as the best candidate for a hiPDGFR ligand/inhibitor.
Although this analysis was carried out with bioinformatics tools and needs experimental evidence and confirmation, we can assert that it has the potential to create a platform for developing new drugs or repositioning already-known drugs. Though the predicted model is an excellent starting point for developing an ideal hiPDGFRα ligand/inhibitor, all 27 selected compounds have the potential to effectively bind and inhibit hiPDGFRα activity based on predicted geometries and affinities. This is partially supported by the fact that compounds such as Imatinib, Nilotinib, Crenolanib, and Nintedanib, which are included in this set, have already been experimentally shown to have activity against hiPDGFR [1,4,8,9,11,15]. Moreover, Bafetinib is already known to inhibit the Bcr/Abl fusion protein tyrosine kinase and the Src-family member Lyn tyrosine kinase for treating Bcr-Abl+ leukemias, including chronic myelogenous leukemia (CML) and Philadelphia+ acute lymphoblastic leukemia [33,34]. Further studies are necessary to determine whether these compounds should proceed to the pre-clinical and clinical phases of evaluation. As reported in a previous study on Nintedanib [35], it will be crucial to carefully adjust the dosage of these new drugs in a real-life setting to ensure their optimal efficacy and safety for patients. This requires close monitoring of potential side effects and close collaboration between healthcare providers and patients to make necessary adjustments based on individual patient factors, such as age, weight, and comorbidities.

4. Materials and Methods

4.1. Virtual Screening of the Small Compounds Database

The structure-based virtual screening of compounds, such as drugs, food constituents, and natural products, against the human intracellular PDGFRα was performed using the MTiOpenScreen web server [36,37] that is based on Autodock Vina-based docking genetic algorithm and analyses 5 compounds libraries: iPPI-lib, diverse-lib, Drugs-lib, Food-lib, and NP-lib. These drug-like chemical libraries contain modulators of protein–protein interactions, diverse chemical compounds, purchasable approved drugs, food constituents, and natural products, respectively, currently for a total of about 7200 molecules. After the preparation of the receptor (pdbID: 6jol) [12] removing Imatinib and water molecules, the PDB file was uploaded on the server setting the grid center coordinates (64.868, 43.016, 0.324) and the size of the search space (30 × 28 × 36 Å), keeping all the other parameters as default.

4.2. Molecular Docking

Twenty-seven molecules were mainly selected on the basis of the best binding energies obtained from the results of the MTiOpenScreen web server. To check these results, we carried out a further docking analysis using the SwissDock web server, as previously reported [38], based on the EADock DSS algorithm [39,40], and setting the docking type as accurate without a definition of the region of interest to include the whole protein in the search space. The final geometry of the complexes obtained was rendered by PyMol software (The PyMOL Molecular Graphics System, Version 2.5.2, Schrödinger, LLC., Cambridge, MA, USA), and the 3D representation of the interaction was obtained by the Protein-Ligand Interaction Profiler (PLIP) web server [41].

4.3. Molecular Dynamics

The molecular dynamics simulation of the best 4 complexes obtained (including the Imatinib/receptor complex as reference) was performed using Gromacs software version 2022.4 [42], the CHARMM36 force field updated July 2022 [43], the SPC216 water model, and the CGenFF server for generating ligand topologies [43]. Each dynamics simulation started with solvation, neutralization adding three chlorine ions, and minimization of the system using the steepest descend algorithm for 50,000 steps. Next, the NVT and NPT equilibration stage was performed in a 100 ps run, stabilizing the temperature at 300 K (using V-rescale Berendsen thermostat) and the pressure at 1 bar (using C-rescale pressure coupling), respectively. Finally, MD simulation was run for 50 ns, and the rms tools of Gromacs were used to determine the root mean square deviation (RMSD) of each ligand trajectory in relation to the protein backbone and, also, to calculate the protein backbone RMSD relative to the energy minimized conformation. The Coulombic and Lennard–Jones interaction energies were calculated using the capability of Gromacs to decompose the short-range non-bonded energies.

4.4. Three-Dimensional Quantitative Structure–Activity Relationship (3D-QSAR) and ADMET Analyses

To simplify the design of selective human intracellular PDGFRα inhibitors, the three-dimensional quantitative structure–activity relationship (3D-QSAR) analysis was carried out using the comparative molecular fields analysis (CoMFA) [44] available on the web portal www.3d-qsar.com [45]. To achieve this analysis, we first performed an alignment of molecules’ conformations and a molecular interaction fields (MIFs) calculation of the 27 selected molecules. We set all parameters as default. The two CoMFA potentials obtained, called steric (STE) and electrostatic (ELE), were calculated by means of the Lennard–Jones and Coulomb law definitions, respectively.
The ADMET (absorption, distribution, metabolism, excretion, and toxicity) analysis was carried out using the screening section of ADMETlab 2.0 web server [46] and uploading the SMILES notation of the 27 selected molecules. In addition to the most common molecular descriptors, the results provide four acceptance indexes called Lipinski, Pfizer, GSK, and Golden Triangle. The Lipinski rule is related to absorption and permeability, while the Pfizer rule is related to toxicity. Compounds satisfying the GSK and Golden Triangle rule may have a more favorable ADMET profile.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24119623/s1.

Author Contributions

Conceptualization, M.M. and G.M.; methodology, software, and validation, M.M.; writing—original draft preparation, M.M. and G.M.; writing—review and editing, M.M., D.B., C.P., S.A., S.S.B. and G.M.; supervision, M.M. and G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received external funding from the Marche Biobank Project POR FESR 2014–2020 to G.M.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Winger, J.A.; Hantschel, O.; Superti-Furga, G.; Kuriyan, J. The structure of the leukemia drug imatinib bound to human quinone reductase 2 (NQO2). BMC Struct. Biol. 2009, 9, 7. [Google Scholar] [CrossRef] [PubMed]
  2. Salem, M.S.H.; Abdel Aziz, Y.M.; Elgawish, M.S.; Said, M.M.; Abouzid, K.A.M. Design, synthesis, biological evaluation and molecular modeling study of new thieno[2,3-d]pyrimidines with anti-proliferative activity on pancreatic cancer cell lines. Bioorg. Chem. 2020, 94, 103472. [Google Scholar] [CrossRef] [PubMed]
  3. Keretsu, S.; Ghosh, S.; Cho, S.J. Molecular Modeling Study of c-KIT/PDGFRalpha Dual Inhibitors for the Treatment of Gastrointestinal Stromal Tumors. Int. J. Mol. Sci. 2020, 21, 8232. [Google Scholar] [CrossRef]
  4. Iqbal, N.; Iqbal, N. Imatinib: A breakthrough of targeted therapy in cancer. Chemother. Res. Pract. 2014, 2014, 357027. [Google Scholar] [CrossRef] [PubMed]
  5. Yu, X.; Zhu, X.; Zhang, L.; Qin, J.J.; Feng, C.; Li, Q. In Silico Screening and Validation of PDGFRA Inhibitors Enhancing Radioiodine Sensitivity in Thyroid Cancer. Front. Pharmacol. 2022, 13, 883581. [Google Scholar] [CrossRef]
  6. Paolini, C.; Agarbati, S.; Benfaremo, D.; Mozzicafreddo, M.; Svegliati, S.; Moroncini, G. PDGF/PDGFR: A Possible Molecular Target in Scleroderma Fibrosis. Int. J. Mol. Sci. 2022, 23, 3904. [Google Scholar] [CrossRef] [PubMed]
  7. Moroncini, G.; Grieco, A.; Nacci, G.; Paolini, C.; Tonnini, C.; Pozniak, K.N.; Cuccioloni, M.; Mozzicafreddo, M.; Svegliati, S.; Angeletti, M.; et al. Epitope Specificity Determines Pathogenicity and Detectability of Anti-Platelet-Derived Growth Factor Receptor alpha Autoantibodies in Systemic Sclerosis. Arthritis Rheumatol. 2015, 67, 1891–1903. [Google Scholar] [CrossRef]
  8. Makino, K.; Makino, T.; Stawski, L.; Mantero, J.C.; Lafyatis, R.; Simms, R.; Trojanowska, M. Blockade of PDGF Receptors by Crenolanib Has Therapeutic Effect in Patient Fibroblasts and in Preclinical Models of Systemic Sclerosis. J. Invest. Dermatol. 2017, 137, 1671–1681. [Google Scholar] [CrossRef]
  9. Huang, J.; Beyer, C.; Palumbo-Zerr, K.; Zhang, Y.; Ramming, A.; Distler, A.; Gelse, K.; Distler, O.; Schett, G.; Wollin, L.; et al. Nintedanib inhibits fibroblast activation and ameliorates fibrosis in preclinical models of systemic sclerosis. Ann. Rheum. Dis. 2016, 75, 883–890. [Google Scholar] [CrossRef]
  10. Benfaremo, D.; Svegliati, S.; Paolini, C.; Agarbati, S.; Moroncini, G. Systemic Sclerosis: From Pathophysiology to Novel Therapeutic Approaches. Biomedicines 2022, 10, 163. [Google Scholar] [CrossRef]
  11. Grzesk, G.; Wozniak-Wisniewska, A.; Blazejewski, J.; Gorny, B.; Wolowiec, L.; Rogowicz, D.; Nowaczyk, A. The Interactions of Nintedanib and Oral Anticoagulants-Molecular Mechanisms and Clinical Implications. Int. J. Mol. Sci. 2020, 22, 282. [Google Scholar] [CrossRef]
  12. Liang, L.; Yan, X.E.; Yin, Y.; Yun, C.H. Structural and biochemical studies of the PDGFRA kinase domain. Biochem. Biophys. Res. Commun. 2016, 477, 667–672. [Google Scholar] [CrossRef] [PubMed]
  13. Tong, M.; Seeliger, M.A. Targeting conformational plasticity of protein kinases. ACS Chem. Biol. 2015, 10, 190–200. [Google Scholar] [CrossRef]
  14. Moroncini, G.; Maccaroni, E.; Fiordoliva, I.; Pellei, C.; Gabrielli, A.; Berardi, R. Developments in the management of advanced soft-tissue sarcoma—olaratumab in context. OncoTargets Ther. 2018, 11, 833–842. [Google Scholar] [CrossRef] [PubMed]
  15. Mendoza, F.A.; Piera-Velazquez, S.; Jimenez, S.A. Tyrosine kinases in the pathogenesis of tissue fibrosis in systemic sclerosis and potential therapeutic role of their inhibition. Transl. Res. 2021, 231, 139–158. [Google Scholar] [CrossRef] [PubMed]
  16. Moroncini, G.; Cuccioloni, M.; Mozzicafreddo, M.; Pozniak, K.N.; Grieco, A.; Paolini, C.; Tonnini, C.; Spadoni, T.; Svegliati, S.; Funaro, A.; et al. Characterization of binding and quantification of human autoantibodies to PDGFRalpha using a biosensor-based approach. Anal. Biochem. 2017, 528, 26–33. [Google Scholar] [CrossRef]
  17. Svegliati, S.; Amico, D.; Spadoni, T.; Fischetti, C.; Finke, D.; Moroncini, G.; Paolini, C.; Tonnini, C.; Grieco, A.; Rovinelli, M.; et al. Agonistic Anti-PDGF Receptor Autoantibodies from Patients with Systemic Sclerosis Impact Human Pulmonary Artery Smooth Muscle Cells Function In Vitro. Front. Immunol. 2017, 8, 75. [Google Scholar] [CrossRef]
  18. Chu, J.S.; Ge, F.J.; Zhang, B.; Wang, Y.; Silvestris, N.; Liu, L.J.; Zhao, C.H.; Lin, L.; Brunetti, A.E.; Fu, Y.L.; et al. Expression and prognostic value of VEGFR-2, PDGFR-beta, and c-Met in advanced hepatocellular carcinoma. J. Exp. Clin. Cancer Res. 2013, 32, 16. [Google Scholar] [CrossRef]
  19. Huang, W.; Fridman, Y.; Bonfil, R.D.; Ustach, C.V.; Conley-LaComb, M.K.; Wiesner, C.; Saliganan, A.; Cher, M.L.; Kim, H.R. A novel function for platelet-derived growth factor D: Induction of osteoclastic differentiation for intraosseous tumor growth. Oncogene 2012, 31, 4527–4535. [Google Scholar] [CrossRef]
  20. Jechlinger, M.; Sommer, A.; Moriggl, R.; Seither, P.; Kraut, N.; Capodiecci, P.; Donovan, M.; Cordon-Cardo, C.; Beug, H.; Grunert, S. Autocrine PDGFR signaling promotes mammary cancer metastasis. J. Clin. Invest. 2006, 116, 1561–1570. [Google Scholar] [CrossRef]
  21. Lin, C.L.; Tsai, M.L.; Chen, Y.H.; Liu, W.N.; Lin, C.Y.; Hsu, K.W.; Huang, C.Y.; Chang, Y.J.; Wei, P.L.; Chen, S.H.; et al. Platelet-Derived Growth Factor Receptor-alpha Subunit Targeting Suppresses Metastasis in Advanced Thyroid Cancer In Vitro and In Vivo. Biomol. Ther. 2021, 29, 551–561. [Google Scholar] [CrossRef]
  22. Lopez-Campistrous, A.; Adewuyi, E.E.; Benesch, M.G.K.; Ko, Y.M.; Lai, R.; Thiesen, A.; Dewald, J.; Wang, P.; Chu, K.; Ghosh, S.; et al. PDGFRalpha Regulates Follicular Cell Differentiation Driving Treatment Resistance and Disease Recurrence in Papillary Thyroid Cancer. EBioMedicine 2016, 12, 86–97. [Google Scholar] [CrossRef]
  23. Matei, D.; Kelich, S.; Cao, L.; Menning, N.; Emerson, R.E.; Rao, J.; Jeng, M.H.; Sledge, G.W. PDGF BB induces VEGF secretion in ovarian cancer. Cancer Biol. Ther. 2007, 6, 1951–1959. [Google Scholar] [CrossRef]
  24. Heldin, C.H.; Lennartsson, J. Structural and functional properties of platelet-derived growth factor and stem cell factor receptors. Cold Spring Harb. Perspect. Biol. 2013, 5, a009100. [Google Scholar] [CrossRef] [PubMed]
  25. Heldin, C.H.; Westermark, B. Mechanism of action and in vivo role of platelet-derived growth factor. Physiol. Rev. 1999, 79, 1283–1316. [Google Scholar] [CrossRef]
  26. Roskoski, R., Jr. The role of small molecule platelet-derived growth factor receptor (PDGFR) inhibitors in the treatment of neoplastic disorders. Pharmacol. Res. 2018, 129, 65–83. [Google Scholar] [CrossRef]
  27. Wynn, T.A. Cellular and molecular mechanisms of fibrosis. J. Pathol. 2008, 214, 199–210. [Google Scholar] [CrossRef]
  28. Alvarez, R.H.; Kantarjian, H.M.; Cortes, J.E. Biology of platelet-derived growth factor and its involvement in disease. Mayo Clin. Proc. 2006, 81, 1241–1257. [Google Scholar] [CrossRef]
  29. Bonner, J.C. Regulation of PDGF and its receptors in fibrotic diseases. Cytokine Growth Factor. Rev. 2004, 15, 255–273. [Google Scholar] [CrossRef] [PubMed]
  30. Trojanowska, M. Role of PDGF in fibrotic diseases and systemic sclerosis. Rheumatology 2008, 47 (Suppl. 5), v2–v4. [Google Scholar] [CrossRef] [PubMed]
  31. Aarthy, M.; Panwar, U.; Selvaraj, C.; Singh, S.K. Advantages of Structure-Based Drug Design Approaches in Neurological Disorders. Curr. Neuropharmacol. 2017, 15, 1136–1155. [Google Scholar] [CrossRef]
  32. Huckleby, A.E.; Saul, J.G.; Shin, H.; Desmarais, S.; Bokka, A.; Jeon, J.; Kim, S.K. Development of Hydroxamic Acid Compounds for Inhibition of Metallo-beta-Lactamase from Bacillus anthracis. Int. J. Mol. Sci. 2022, 23, 9163. [Google Scholar] [CrossRef]
  33. Santos, F.P.; Kantarjian, H.; Cortes, J.; Quintas-Cardama, A. Bafetinib, a dual Bcr-Abl/Lyn tyrosine kinase inhibitor for the potential treatment of leukemia. Curr. Opin. Investig. Drugs 2010, 11, 1450–1465. [Google Scholar]
  34. Kantarjian, H.; le Coutre, P.; Cortes, J.; Pinilla-Ibarz, J.; Nagler, A.; Hochhaus, A.; Kimura, S.; Ottmann, O. Phase 1 study of INNO-406, a dual Abl/Lyn kinase inhibitor, in Philadelphia chromosome-positive leukemias after imatinib resistance or intolerance. Cancer 2010, 116, 2665–2672. [Google Scholar] [CrossRef]
  35. Campochiaro, C.; De Luca, G.; Lazzaroni, M.G.; Armentaro, G.; Spinella, A.; Vigone, B.; Ruaro, B.; Stanziola, A.; Benfaremo, D.; De Lorenzis, E.; et al. Real-life efficacy and safety of nintedanib in systemic sclerosis-interstitial lung disease: Data from an Italian multicentre study. RMD Open. 2023, 9, e002850. [Google Scholar] [CrossRef]
  36. Labbe, C.M.; Rey, J.; Lagorce, D.; Vavrusa, M.; Becot, J.; Sperandio, O.; Villoutreix, B.O.; Tuffery, P.; Miteva, M.A. MTiOpenScreen: A web server for structure-based virtual screening. Nucleic Acids Res. 2015, 43, W448–W454. [Google Scholar] [CrossRef]
  37. Lagarde, N.; Rey, J.; Gyulkhandanyan, A.; Tuffery, P.; Miteva, M.A.; Villoutreix, B.O. Online structure-based screening of purchasable approved drugs and natural compounds: Retrospective examples of drug repositioning on cancer targets. Oncotarget 2018, 9, 32346–32361. [Google Scholar] [CrossRef]
  38. Ferraro, G.; Mozzicafreddo, M.; Ettari, R.; Corsi, L.; Monti, M.C. A Proteomic Platform Unveils the Brain Glycogen Phosphorylase as a Potential Therapeutic Target for Glioblastoma Multiforme. Int. J. Mol. Sci. 2022, 23, 8200. [Google Scholar] [CrossRef]
  39. Grosdidier, A.; Zoete, V.; Michielin, O. SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 2011, 39, W270–W277. [Google Scholar] [CrossRef]
  40. Grosdidier, A.; Zoete, V.; Michielin, O. Fast docking using the CHARMM force field with EADock DSS. J. Comput. Chem. 2011, 32, 2149–2159. [Google Scholar] [CrossRef]
  41. Adasme, M.F.; Linnemann, K.L.; Bolz, S.N.; Kaiser, F.; Salentin, S.; Haupt, V.J.; Schroeder, M. PLIP 2021: Expanding the scope of the protein-ligand interaction profiler to DNA and RNA. Nucleic Acids Res. 2021, 49, W530–W534. [Google Scholar] [CrossRef] [PubMed]
  42. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef] [PubMed]
  43. Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; et al. CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. J. Comput. Chem. 2010, 31, 671–690. [Google Scholar] [CrossRef] [PubMed]
  44. Cramer, R.D.; Patterson, D.E.; Bunce, J.D. Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. J. Am. Chem. Soc. 1988, 110, 5959–5967. [Google Scholar] [CrossRef]
  45. Ragno, R. www.3d-qsar.com: A web portal that brings 3-D QSAR to all electronic devices-the Py-CoMFA web application as tool to build models from pre-aligned datasets. J. Comput. Aided Mol. Des. 2019, 33, 855–864. [Google Scholar] [CrossRef]
  46. Xiong, G.; Wu, Z.; Yi, J.; Fu, L.; Yang, Z.; Hsieh, C.; Yin, M.; Zeng, X.; Wu, C.; Lu, A.; et al. ADMETlab 2.0: An integrated online platform for accurate and comprehensive predictions of ADMET properties. Nucleic Acids Res. 2021, 49, W5–W14. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the methodology and workflow implemented for the SBDD.
Figure 1. Schematic representation of the methodology and workflow implemented for the SBDD.
Ijms 24 09623 g001
Figure 2. Complex geometries prediction by molecular docking analysis: amino acids and the ligand of (A) Bafetinib/hiPDGFRα and (B) Flumatinib/hiPDGFRα complex involved in the interaction are shown as sticks. Weak bonds are also reported and emphasized in the legends.
Figure 2. Complex geometries prediction by molecular docking analysis: amino acids and the ligand of (A) Bafetinib/hiPDGFRα and (B) Flumatinib/hiPDGFRα complex involved in the interaction are shown as sticks. Weak bonds are also reported and emphasized in the legends.
Ijms 24 09623 g002
Figure 3. RMSDs of ligands and complexes trajectories: (A,C,E,G) RMSDs of Bafetinib, Radotinib, Imatinib, and Flumatinib versus protein backbone, respectively; (B,D,F,H) RMSDs of hiPDGFRα versus protein backbone complexed with Bafetinib, Radotinib, Imatinib, and Flumatinib, respectively.
Figure 3. RMSDs of ligands and complexes trajectories: (A,C,E,G) RMSDs of Bafetinib, Radotinib, Imatinib, and Flumatinib versus protein backbone, respectively; (B,D,F,H) RMSDs of hiPDGFRα versus protein backbone complexed with Bafetinib, Radotinib, Imatinib, and Flumatinib, respectively.
Ijms 24 09623 g003
Figure 4. CoMFA model prediction values of activities (Recalcd/Pred pAct) versus experimental activities obtained by molecular docking analysis considering the optimal PC values.
Figure 4. CoMFA model prediction values of activities (Recalcd/Pred pAct) versus experimental activities obtained by molecular docking analysis considering the optimal PC values.
Ijms 24 09623 g004
Table 1. The best compounds selected after the structure-based virtual screening and their affinities for the hiPDGFRα obtained after the molecular docking analysis by the SwissDock web server.
Table 1. The best compounds selected after the structure-based virtual screening and their affinities for the hiPDGFRα obtained after the molecular docking analysis by the SwissDock web server.
CompoundDatabaseΔG (kcal/mol)Kd,pred (M)
BafetinibDrugs-lib−12.261.01 × 10−9
RadotinibDrugs-lib−11.742.43 × 10−9
ImatinibDrugs-lib−11.672.74 × 10−9
FlumatinibDrugs-lib−11.642.88 × 10−9
NilotinibDrugs-lib−10.939.53 × 10−9
Ag-13958Drugs-lib−10.791.21 × 10−8
Tg100-801Drugs-lib−10.761.28 × 10−8
DitercaliniumDrugs-lib−10.153.60 × 10−8
DasatinibDrugs-lib−10.064.14 × 10−8
R428_BemcentinibDrugs-lib−9.984.78 × 10−8
GlisolamideDrugs-lib−9.915.35 × 10−8
MolPort-002-524-598NP-lib−9.836.19 × 10−8
Bms-833923Drugs-lib−9.531.01 × 10−7
BenfluorexDrugs-lib−9.401.28 × 10−7
MolPort-039-052-621NP-lib−9.231.69 × 10−7
MolPort-009-018-791_CurcuminNP-lib−9.102.12 × 10−7
47194043_iPPIiPPI-lib−9.082.18 × 10−7
CrenolanibDrugs-lib−9.082.18 × 10−7
MolPort-001-740-946_luteolin-7-O-glucuronideNP-lib−9.072.21 × 10−7
NintedanibDrugs-lib−8.803.48 × 10−7
24301892_iPPIiPPI-lib−8.783.60 × 10−7
24824231_divDiverse-lib−8.743.88 × 10−7
MolPort-001-741-358_EGCGNP-lib−8.704.18 × 10−7
29215783_iPPIiPPI-lib−8.486.00 × 10−7
MolPort-001-740-557_QuercetinNP-lib−8.406.88 × 10−7
Hispaglabridin_BFood-lib−8.367.38 × 10−7
MolPort-001-768-161_BinaphthaleneNP-lib−8.061.23 × 10−6
Table 2. Average short-range interaction energies calculated by energy module of Gromacs.
Table 2. Average short-range interaction energies calculated by energy module of Gromacs.
ComplexCoulombic Interaction Energy (kJ/mol)Lennard–Jones Energy (kJ/mol)
Bafetinib/hiPDGFRα−135.72 ± 0.96−257.94 ± 0.60
Radotinib/hiPDGFRα−142.66 ± 1.10−234.45 ± 1.50
Imatinib/hiPDGFRα−121.95 ± 1.20−234.49 ± 1.80
Flumatinib/hiPDGFRα−135.79 ± 2.10−244.41 ± 3.00
Table 3. Summary of the best CoMFA models.
Table 3. Summary of the best CoMFA models.
Fieldr2q2Optimal PC
STE1.0000.5255
ELE0.7610.5744
STE-ELE0.9990.6926
Table 4. ADMET results. Indexes and some descriptors are reported (part 1).
Table 4. ADMET results. Indexes and some descriptors are reported (part 1).
CompoundLipinskiPfizerGSKGolden TriangleMWLogPLogDnHAnHDTPSA
Recommended Range 100/6000/31/30/120/70/140
BenfluorexAcceptedRejectedRejectedAccepted351.144.2294.0393138.33
GlisolamideAcceptedAcceptedRejectedAccepted434.162.9531.32593130.4
ImatinibAcceptedAcceptedRejectedAccepted493.263.8053.1448289.51
NintedanibAcceptedAcceptedRejectedRejected539.253.4663.03692101.47
NilotinibAcceptedAcceptedRejectedRejected529.184.8943.88182100.85
RadotinibAcceptedAcceptedRejectedRejected530.184.5743.64392113.74
FlumatinibAcceptedAcceptedRejectedRejected562.244.0433.2792102.4
Ag-13958AcceptedAcceptedRejectedAccepted467.194.6253.47683100.52
R428_BemcentinibAcceptedAcceptedRejectedRejected506.294.544.01783101.74
BafetinibAcceptedAcceptedRejectedRejected576.264.2763.51892102.4
DitercaliniumRejectedRejectedRejectedRejected718.47.5625.0148264.28
Bms-833923AcceptedAcceptedRejectedAccepted473.225.1933.9026382.17
Tg100-801RejectedAcceptedRejectedRejected579.26.7454.8578192.7
Hispaglabridin_BAcceptedRejectedRejectedRejected390.186.7575.3244147.92
CrenolanibAcceptedAcceptedRejectedAccepted443.234.343.0937278.43
DasatinibAcceptedAcceptedRejectedAccepted487.162.8072.92293109.74
MolPort-001-740-557_QuercetinAcceptedAcceptedAcceptedAccepted302.042.1551.76775131.36
MolPort-001-768-161_BinaphthaleneAcceptedRejectedRejectedAccepted254.116.024.668000
MolPort-001-741-358_EGCGRejectedAcceptedRejectedAccepted458.081.8930.652118197.37
MolPort-009-018-791_CurcuminAcceptedAcceptedAcceptedAccepted368.132.7422.826293.06
MolPort-039-052-621RejectedAcceptedRejectedRejected542.125.0213.33105162.98
MolPort-002-524-598RejectedAcceptedRejectedRejected596.14−0.524−0.31610269.43
MolPort-001-740-946_Luteolin-7-glucuronideRejectedAcceptedRejectedAccepted462.080.8640.745127207.35
29215783_iPPIAcceptedRejectedRejectedAccepted410.085.884.1674048.67
24301892_iPPIAcceptedAcceptedRejectedAccepted454.135.4254.1047177.75
47194043_iPPIAcceptedRejectedRejectedAccepted354.144.8243.5254261.96
24824231_divAcceptedAcceptedAcceptedAccepted330.12.8632.8295275.27
Abbreviations: molecular weight (MW), logarithm of the n-octanol/water distribution coefficient (LogP), logarithm of the n-octanol/water distribution coefficients at pH = 7.4 (LogD), number of hydrogen bond acceptors (nHA), number of hydrogen bond donors (nHD), and topological polar surface area (TPSA).
Table 5. ADMET results. Indexes and some descriptors are reported (part 2).
Table 5. ADMET results. Indexes and some descriptors are reported (part 2).
CompoundQEDNP-LikenessCaco-2MDCKCarcinogenicitySkinSenECEIRespiratoryLD50 Oral
Recommended Range>0.67−5/5>−5.15>2.0 × 10−60/0.3 Excellent; 0.3/0.7 Medium; 0.7/1 Poor0
Benfluorex0.603−0.95−4.622.19 × 10−50.0460.0330.0030.0110.6860
Glisolamide0.613−1.57−5.4471.56 × 10−50.0260.1090.0030.0070.0070
Imatinib0.424−1.477−5.5767.27 × 10−60.0440.9460.0030.0070.9920
Nintedanib0.272−1.176−5.7162 × 10−50.0780.750.0030.0070.6940
Nilotinib0.302−1.657−5.1121.12 × 10−50.0270.9320.0030.0090.9280
Radotinib0.325−1.62−5.0651.06 × 10−50.0210.9440.0030.0090.9540
Flumatinib0.361−1.566−5.4927.93 × 10−60.0280.9480.0030.0060.9880
Ag-139580.314−1.989−5.466.03 × 10−60.5910.2450.0030.0130.9860
R428_Bemcentinib0.364−0.866−5.3426.99 × 10−60.950.5810.0030.0070.740
Bafetinib0.327−1.607−5.6585.79 × 10−60.0360.9490.0030.0060.950
Ditercalinium0.1490.003−5.9496.59 × 10−60.0860.9620.0030.0090.9990
Bms-8339230.294−0.978−5.5261.54 × 10−50.1690.3920.0030.010.8960
Tg100-8010.166−1.068−5.1931.92 × 10−50.2540.4770.0030.0090.7430
Hispaglabridin_B0.72.088−5.061.87 × 10−50.9060.8540.0030.1110.8840
Crenolanib0.504−0.993−5.4219.76 × 10−50.8250.5430.0030.0080.810
Dasatinib0.493−1.724−4.8711.25 × 10−50.9480.1390.0030.0080.8460
MolPort-001-740-557_Quercetin0.4341.701−5.2047.69 × 10−50.050.9190.0070.9360.0720
MolPort-001-768-161_Binaphthalene0.411−0.156−4.6911.45 × 10−50.8610.9450.0110.9830.0190
MolPort-001-741-358_EGCG0.2121.65−6.7175.05 × 10−60.0340.9690.0030.9360.0210
MolPort-009-018-791_Curcumin0.5480.722−4.8341.63 × 10−50.7060.9580.0070.7920.9510
MolPort-039-052-6210.2250.947−5.0592.03 × 10−50.390.940.0030.9120.8160
MolPort-002-524-5980.1371.984−6.4481.21 × 10−50.2430.6740.0030.2220.0180
MolPort-001-740-946_Luteolin-7-glucuronide0.2531.771−6.4713.08 × 10−50.2790.0790.0030.0130.0580
29215783_iPPI0.388−0.195−4.9163.44 × 10−50.5750.0330.0030.5130.7040
24301892_iPPI0.471−1.923−4.9421.95 × 10−50.5540.0420.0030.0540.8380
47194043_iPPI0.499−0.597−5.0347.22 × 10−60.4650.8820.0030.6230.6840
24824231_div0.725−1.188−5.2891.36 × 10−50.4270.8540.0030.0360.0170
Abbreviations: drug-likeness based on the concept of desirability (QED), Natural Product-likeness score (NP-likeness), human colon adenocarcinoma cell lines permeability (Caco-2), Madin−Darby Canine Kidney cells permeability (MDCK), skin sensitization (Skin Sen), eye irritation/corrosion (EI/EC) potential, respiratory toxicity (Respiratory), and acute toxicity during oral administration (LD50 oral).
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

Mozzicafreddo, M.; Benfaremo, D.; Paolini, C.; Agarbati, S.; Svegliati Baroni, S.; Moroncini, G. Screening and Analysis of Possible Drugs Binding to PDGFRα: A Molecular Modeling Study. Int. J. Mol. Sci. 2023, 24, 9623. https://doi.org/10.3390/ijms24119623

AMA Style

Mozzicafreddo M, Benfaremo D, Paolini C, Agarbati S, Svegliati Baroni S, Moroncini G. Screening and Analysis of Possible Drugs Binding to PDGFRα: A Molecular Modeling Study. International Journal of Molecular Sciences. 2023; 24(11):9623. https://doi.org/10.3390/ijms24119623

Chicago/Turabian Style

Mozzicafreddo, Matteo, Devis Benfaremo, Chiara Paolini, Silvia Agarbati, Silvia Svegliati Baroni, and Gianluca Moroncini. 2023. "Screening and Analysis of Possible Drugs Binding to PDGFRα: A Molecular Modeling Study" International Journal of Molecular Sciences 24, no. 11: 9623. https://doi.org/10.3390/ijms24119623

APA Style

Mozzicafreddo, M., Benfaremo, D., Paolini, C., Agarbati, S., Svegliati Baroni, S., & Moroncini, G. (2023). Screening and Analysis of Possible Drugs Binding to PDGFRα: A Molecular Modeling Study. International Journal of Molecular Sciences, 24(11), 9623. https://doi.org/10.3390/ijms24119623

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop