Molecular Docking Studies of a Cyclic Octapeptide-Cyclosaplin from Sandalwood
1
Plant Biotechnology Laboratory, Department of Biotechnology, Indian Institute of Technology Kharagpur, Kharagpur, West Bengal 721302, India
2
Department of Internal Medicine, The University of Texas Southwestern Medical Center, 5323 Harry Hines Blvd, Dallas, TX 75390, USA
*
Author to whom correspondence should be addressed.
Biomolecules 2019, 9(11), 740; https://doi.org/10.3390/biom9110740
Submission received: 19 October 2019
/
Accepted: 12 November 2019
/
Published: 15 November 2019
(This article belongs to the Section Bioinformatics and Systems Biology)
Abstract
:Natural products from plants, such as chemopreventive agents, attract huge attention because of their low toxicity and high specificity. The rational drug design in combination with structure-based modeling and rapid screening methods offer significant potential for identifying and developing lead anticancer molecules. Thus, the molecular docking method plays an important role in screening a large set of molecules based on their free binding energies and proposes structural hypotheses of how the molecules can inhibit the _target. Several peptide-based therapeutics have been developed to combat several health disorders, including cancers, metabolic disorders, heart-related diseases, and infectious diseases. Despite the discovery of hundreds of such therapeutic peptides however, only few peptide-based drugs have made it to the market. Moreover, the in silico activities of cyclic peptides towards molecular _targets, such as protein kinases, proteases, and apoptosis related proteins have not been extensively investigated. In this study, we explored the in silico kinase and protease inhibitor potentials of cyclosaplin, and studied the interactions of cyclosaplin with other apoptosis-related proteins. Previously, the structure of cyclosaplin was elucidated by molecular modeling associated with dynamics that were used in the current study as well. Docking studies showed strong affinity of cyclosaplin towards cancer-related proteins. The binding affinity closer to 10 kcal/mol indicated efficient binding. Cyclosaplin showed strong binding affinities towards protein kinases such as EGFR, VEGFR2, PKB, and p38, indicating its potential role in protein kinase inhibition. Moreover, it displayed strong binding affinity to apoptosis-related proteins and revealed the possible role of cyclosaplin in apoptotic cell death. The protein–ligand interactions using LigPlot displayed some similar interactions between cyclosaplin and peptide-based ligands, especially in case of protein kinases and a few apoptosis related proteins. Thus, the in silico analyses gave the insights of cyclosaplin being a potential apoptosis inducer and protein kinase inhibitor.
1. Introduction
Cancer is a well-recognized global health problem responsible for ∼7.6 million deaths (∼13% of all deaths) worldwide, which is expected to rise to 13.1 million by 2030 (WHO, 2012). Despite the advancements in the field of cancer research, there is still an urgency to discover and develop anti-cancer therapeutics. Natural products are of particular interest as chemopreventive agents because of their low toxicities and potential efficacies [1]. The conventional drug discovery techniques are time consuming and expensive processes [2]. Thus, rational drug design in combination with structure based modeling and rapid screening methods offer significant potential for identifying and developing lead anticancer molecules. The use of the molecular docking method addresses deducing the ligand binding sites with a protein of known three-dimensional structure. One of the computational approaches, docking, helps with screening a large set of molecules based on their free binding energies and proposes structural hypotheses of how the molecules could inhibit the _target. For example, docking studies were used to select and rationally design novel biguanides towards m/hTAAR1 (Guariento et al., 2018) [3]. Additionally, docking studies were used to identify selective and 5-HT1A receptor agonists (Franchin et al., 2017) [4]. Recently, several in silico-based studies have been performed on small molecules, including peptides, to identify their anti-cancerous properties [5]. Peptide based therapeutics have been effective at combating several health disorders, including cancers, metabolic disorders, cardiovascular diseases, neurological disorders, kidney diseases, and infectious diseases. Peptides are structurally diverse, have a wide spectrum of therapeutic action, have low absorption in body tissues, and are highly specific to _targets [6]. Several cyclic peptides with diverse biological activities, such as antibacterial activity, immunosuppressive activity, and anticancer activity, have been reported [7]. For example, tyrocidine and gramicidin S with antibacterial activity; cyclosporin A, displaying immunosuppressive activity; and Cyclo-RGDfV, with antiangiogenic activity [7,8,9]. Apart from their use as cytotoxic agents, peptides can also be used in drug formulations for enhancing biological activity, _targeted drug delivery, or transport across cellular membranes. Thus, revival of interest in therapeutic peptides and extensive research has seen peptides entering into clinical trials improve significantly over the decade [10]. Despite the discovery of hundreds of such therapeutic peptides, however, only few peptide-based drugs have made it to the market. Moreover, thus far, the activities of cyclic peptides towards molecular _targets such as protein kinases, proteases, and apoptosis related proteins, have never been explored. In this study we explorde the in silico kinase and protease inhibitor potentials of cyclosaplin and studied the interactions of cyclosaplin with other cancer-related proteins.
2. Materials and Methods
2.1. Softwares and Tools
ACD/ChemSketch 12.01, AutoDock Vina 1.1.2, Avogadro, CycloPsWeb, GROMACS, LigPlot, Modeller 9.2, MGL tools, Open Babel, Protein Data Bank (PDB), PubChem, PyMOL, and Swiss _target Prediction.
2.2. Ligand Preparation
The cyclic octapeptide (cyclosaplin) and various peptides (positive control) for specific proteins were used as ligands for docking studies (Table 1). The peptides (positive controls) were chosen based on their specific modes of action on cancer receptors. A molecular dynamics tool GROMACS and Modeller 9.2 program were used for the molecular modeling of cyclosaplin using appropriate energy minimization steps and simulations previously described [11,12]. Due to the unavailability of X-ray diffraction or NMR structure data for the chosen ligands, the ligand molecules were drawn in either ACD/Chem Basic freeware (ACD/ChemSketch 12.01) or using CycloPsWeb and saved as MDL mol file formats. The MDL files were converted to pdb format files using Open Babel. Further, all the ligand structures were energy-minimized using GROMACS and Modeller 9.2 program prior to docking studies [13]. The ligands used in the study are represented in the Table 1.
2.3. Lipinski Rule for Ligands
The peptide based-ligand molecules selected for docking experiments were screened for Lipinski’s rule of five. Lipinski’s rule of five [23] states that a drug molecule generally does not violate more than one of the following five rules
- o
- Molecular mass less than 500 Da;
- o
- High lipophilicity (expressed as LogP less than 5);
- o
- Less than 5 hydrogen bond donors;
- o
- Less than 10 hydrogen bond acceptors;
- o
- Molar refractivity between 40 and 130.
Lipinski’s rule of five was also checked in Supercomputing Facility for Bioinformatics and Computational Biology, IIT Delhi, wherein PDB structures of the molecules were uploaded to the online server (http://www.scfbio-iitd.res.in/utility/LipinskiFilters.jsp).
2.4. Protein Preparation
Swiss _target Prediction was used to predicting the potential _targets of cyclosaplin [24]. The protein structures were obtained from Protein Data Bank (PDB) [25]. The proteins selected for the study were epidermal growth factor’s receptor kinase domain (EGFR; PDB ID: 2GS2), vascular endothelial growth factor r2’s receptor kinase (VEGFR2; PDB ID: 1VR2), mitogen activated protein kinase (P38; PDB ID: 1P38), protein kinase B (PKB; PDB ID: 1GZN), phosphatase and tensin homolog tumor suppressor (PTEN; PDB ID:1D5R), matrix metalloproteases (MMP-2 (PDB ID: 1CK7), and MMP-9 (PDB ID: 1L6J)); and apoptosis related proteins (Procaspase 3 (PDB ID: 4JQZ), procaspase 7 (PDB ID: 1K88), caspase 9 (PDB ID: 2AR9), TRAIL (PDB ID: 1D2Q), and SURVIVIN (PDB ID: 1XOX)). EGFR kinase and procaspase 3 were previously used in our study and were used as controls in this study [11]. The files in pdb format for each receptor were converted to respective PDBQT format using MGL tools. The polar hydrogen atoms were added to the receptor molecules prior to docking studies. Three-dimensional affinity grids were created at the geometric centers of the _target proteins.
2.5. Docking Studies Using AutoDock Vina
The energy-minimized structure of cyclic octapeptide and ligands o (positive control) were docked with _target proteins using AutoDock Vina 1.1.2 [26]. The receptor and ligand files were represented in PDBQT file format, a modified pdb format containing atomic charges, atom type definitions for ligands, and topological information (rotatable bonds). For docking, the entire receptor was enclosed inside a grid box, with a grid spacing of 1 Å, keeping the receptor rigid and the ligand as a flexible molecule. The ligand’s backbone and side-chain were flexible and allowed to dock with the receptor to form all possible conformations. After defining the binding site and receptor–ligand preparation, docking runs were launched from the command prompt. The interaction energy between the ligand and the receptor was calculated for the entire binding site and expressed as affinity (kcal/mol).
2.6. Protein–Ligand Interactions
LigPlot was used to study protein–ligand interactions for a given pdb file encrypting the docking [27]. The LigPlot program self-generated schematic 2D representations of the interfaces of protein–ligand complexes from standard pdb file input. The output was in the form of informative representation of the intermolecular interactions and their strengths, including hydrogen bonds, hydrophobic contacts, and atom accessibilities. Hydrogen bonds are represented by dashed lines whereas hydrophobic contacts are depicted schematically. The amino acid residues of the protein involved with the above contacts are shown by an arc with spokes emerging towards the ligand atoms in contact and vice versa.
3. Results
3.1. Ligand Preparation
The ligand structures were drawn in CycloPsWeb or downloaded from PubChem and converted to pdb format using Open Babel and modeled using GROMACS and MODELLER 9.2 program along with cyclosaplin (Table 2). The structures were energy minimized and saved in PDBQT format by MGL tools (Figure 1). The energy minimized structure of cyclosaplin is represented in Figure S1.
3.2. Lipinski Rule
The ligands prepared for docking were screened for Lipinski’s rule of five. The commercially available or reported peptide inhibitors/inducers (positive control) were also tested against the _target proteins respectively (Table 3).
3.3. Protein Preparation
Swiss _target Prediction was used to predict the possible _targets of cyclosaplin (Figure 2a). A list of 50 potential _targets was displayed, out of which top _targets were selected from each class, such as kinases (33%), membrane receptors (40%), and proteases (27%). The cancer related proteins were downloaded from Protein Data Bank and (PDB) converted to PDBQT format using AutoDock tools (Figure 2b).
3.4. Docking Studies Using AutoDock Vina
The docking process was carried out using AutoDock Vina. The docking scores are graphically represented in Figure 3a and the binding affinities of ligands are represented as kcal/mol (Table 4). The affinity value of less than or closer to 5 kcal/mol depicts negligible binding, whereas values closer to 10 kcal/mol indicate efficient binding (Figure 3a). The protein–cyclosaplin docking studies are represented in Figure 3b, and protein–ligand interactions (positive control) are represented in Figure 4. Cyclosaplin exhibited stronger binding affinity (>5 kcal/mol for all the cancer-related proteins).
3.5. Protein–Ligand Interactions
The protein–ligand interaction study was performed using LigPlot. The interactions of the ligands cyclosaplin and various peptide-based ligands with amino acids residues of the _target proteins are shown in Table 5. The H-bonds and hydrophobic contacts between the docked complexes are shown in Figure 5 and Figure 6.
4. Discussion
Peptides are effective receptor-binding ligands; many other classes of ligands sharing this binding trait include small molecules, endogenous proteins, and antibodies [28]. Cyclic peptides have built-in, stable pharmacokinetic characteristics, including enzyme stability, conformational rigidity, improved receptor site selectivity, and pharmacological specificity. In addition, cyclic peptides are reported to be potent protein kinase inhibitors, protease inhibitors (MMP-2 and MMP-9), angiogenesis blockers, and apoptosis inducers [29,30,31,32]. In comparison to small molecules, cyclic peptides can be more selective, whereas the size of molecules can be less than protein molecules such as antibodies and growth factors. Therefore, in the present study, an attempt was made to investigate the potential of cyclosaplin and other reported peptide-based ligands (positive control) against specific cancer-related proteins.
In our previous study, cyclosaplin was isolated, purified, and characterized from Santalum album L. [11]. The cyclosaplin was molecularly modeled and the energy minimized structure was further used for docking studies (Figure S1). The ligands were energy minimized prior to docking studies (Table 1 and Table 2, Figure 1). All of the peptide-based ligands, along with cyclosaplin, were screened for Lipinski’s rule of five (Table 3). Some of these peptides violated the rules, yet displayed drug-like properties in the experimental studies in vitro. Cyclic peptides tend to have properties (e.g., MW, number of polar atoms, and total polar surface area) that put them outside conventional predictors of “drug-likeness,” such as Lipinski’s rule of five [23]. In spite of this, many compounds exhibited drug-like properties, including the potential to penetrate cellular membranes.
The potential _targets of cyclosaplin were predicted by Swiss _target Prediction [23] (Figure 2a) and the proteins used in docking studies were energy minimized, which is represented in Figure 2b. Relative binding affinities were scored for the cyclosaplin and peptide-based ligands, represented as kcal/mol (Table 4). The affinity value of less than five depicts negligible binding, whereas values closer to 10 kcal/mol indicate efficient binding. In addition, the docking scores for various cancer-related proteins was represented graphically, as shown in Figure 3. Docking studies revealed the strong binding affinities of cyclosaplin towards apoptosis-related proteins procaspase 3 (−7.8 kcal/mol; [11]), procaspase 7 (−8.7 kcal/mol), caspase 9 (−8.9 kcal/mol), TRAIL (−8.2 kcal/mol), SURVIVIN (−7.4 kcal/mol), and protease MMP-2 (−8.2 kcal/mol) (Figure 3a,b). Cyclosaplin also demonstrated effective binding affinities towards other cancer-related proteins, such as EGFR (−6.8 kcal/mol) [9], VEGFR2 (−7.8 kcal/mol), PKB (−8.1 kcal/mol), p38 (−8.3 kcal/mol), PTEN-tumor suppressor (−6.3 kcal/mol), and MMP-9 (−7.3 kcal/mol) (Table 4, Figure 3). The peptide-based ligands (positive control) reported in the literature or under clinical studies showed strong binding affinities with the specific proteins except for TRAIL (Figure 3). In case of TRAIL, the ligand remained unbound to the protein with a score of −6.4 kcal/mol. The result indicated the possible role of cyclosaplin in mediating apoptotic cell death. Cyclosaplin exhibited stronger binding affinity (>5 kcal/mol for all the protein _targets which is consistent with our previously shown experimental study were we have shown that the cyclosaplin exhibits significant anti-proliferative activity with an IC50 2.06 µg/mL in MDA-MB-231 cells (Mishra et al., 2014).
In contrast to most small molecule drugs, peptides have high affinity, strong specificity for _targets, and low toxicity, whereas, in contrast to chemotherapeutics antibodies, they have good penetration of tissues because of their small size [33,34,35,36]. Cyclization is also thought to minimize conformational entropy losses upon _target binding, although some studies have shown the impact of cyclization on binding entropy to be more complex [37].
The interaction of the cyclosaplin and other peptide-based ligands (positive control) with the amino acids of various cancer-related proteins were also determined (Table 5). We previously showed the structure–activity relationship for EGFR kinase with cyclosaplin [11], but in the present study, we demonstrated the possible interactions between protein and ligand with key amino acid residues involved in such interactions. In case of EGFR kinase, the peptide inhibitor CVRACGAD (cyclic) showed no similar interactions with cyclosaplin for amino acid residues of the protein (Table 5). The cyclosaplin interacted with Asp 960/Glu961 and Ser787/Tyr789, forming H-bonds and hydrophobic contacts respectively (Table 5, Figure 5A and Figure S2). Asp-960/Glu-961 facilitates the movement of the C-terminal tail of the EGF receptor to regulate asymmetrical dimer formation [27]. The side chain of Asp-960 interacts with that of Ser787, and mutation at this site enhanced protein kinase activity [38], whereas Tyr789 is the site for phosphorylation, the new potential binding site from the catalytic domain of EGFR [39]. The positive control CVRACGAD (Figure 6A and Figure S2B) forms a H-bond with Lys 721, whose side chains interact with the ATP forming salt-bridges in activated kinases [40]. In addition, it interacts with the glycine-rich nucleotide phosphate-binding loop (Gly695-Gly700) and DFG motif (Asp831-Gly833) within the A-loop [40]. The interaction between cyclosaplin and EGFR kinase occurs on Asp960, Glu961, Ser787, and Tyr789 with significant binding affinity. The residues mentioned above played a key role in dimer formation and the site for phosphorylation, highlighting that cyclosaplin could inhibit EGFR kinase by interacting with C-terminal region of EGFR (Table 5, Figure 5A). Previously, we have shown binding of cyclosaplin to EGFR through co-localization studies resulting in the sensitization of MDA-MB-231 cells and the induction of apoptosis [11]. Interestingly, certain common amino acid residues of most of the proteins shared trans-similarity; for example, residues involved in H-bond formation in cyclosaplin matched with residues forming hydrophobic contacts in peptide-based ligands. It is not necessary for hydrophobic interactions to occur only between the amino acids with hydrophobic side chains. It can occur between all the amino acid residues depending on their total degrees of hydrophobicity [41]. The architecture of VEGFR-2 involves several important loop domains, including glycine-rich loop (also refers to nucleotide binding loop) at residues 841–846, the catalytic loop at residues 1026–1033, and the activation loop at residues 1046–1075 [42]. The active sites around the ATP-binding domain of VEGFR-2 consist of three hydrophobic regions (regions 1–3) and one polar region (region 4). Between region 1 and the region 2, Lys866, Glu883, and Asp1044 are crucial for receptor activation [42]. Region 3 contains only a few residues, including Leu838 and Phe916. The unique polar region involves several residues, such as Asn921, Cys1043, Arg1030, and Asn1031 [43]. The interaction between cyclosaplin and VEGFR2 occurs on Glu 885, Asn923, Asp1046, Cys919, and Lys868 with strong binding affinity indicating that cyclosaplin could inhibit VEGFR-2 activity by interacting with the ATP-binding site of VEGFR-2 (Table 5, Figure 5B, Figure 6B, and Figure S2). A similar residual interaction occurred in the case of antiangiogenic peptide, cilengitide, and VEGFR2 kinase. We envisage that a fixed geometry ascertained due to cyclization in peptides could help it bind to receptors more effectively. The RGD peptide or RGD-like peptides are good examples of cyclic peptides as receptor binding molecules (Figure S4). The binding affinity of cyclosaplin towards α5β3 was closer to 10 kcal/mol (−9.5 kcal/mol), indicating strong binding (Table S1). Some common amino acid residues, such as Arg274, Asp275, Lys181, Phe163, and Thr199 of PKB interacted with both cyclosaplin (RLGDGCTR) and RPRTSSF (Table 5, Figure 5C, Figure 6C, Figures S2 and S3). Mutational analysis of Arg274 in Akt2 is essential for shielding Ther308 in the activation loop against dephosphorylation [44]. The α-helix at C-terminal (αC helix) of the N lobe plays a vital role in regulating the catalytic functions in all the protein kinases [45]. In the inactive state of PKB, His 196, and Glu 200 of the αC helix are disordered, and contacts between Glu 200 and Lys 181, and those between His 196 and pThr 309, are not formed [45]. The interaction between cyclosaplin and PKB occurs on Arg274, Lys181, Phe163, Thr199, Tyr273, and Leu183 with strong binding, indicating its possible role as PKB inhibitor (Table 5, Figure 5C). Moreover, the above interacting amino acid residues are also common to RPRTSSF, the positive control used in this study (Figure 6C). In p38, both the peptide-based ligands (FWCS and cyclosaplin) had interactions with common amino acid residues involved in phosphate and ATP binding sites (Table 5, Figure 5D, Figure 6D, Figures S2 and S3). Among all of the MAP kinases, the phosphorylation sites (Thr-180 and Tyr-182), and the putative phosphate binding ligands (Arg67, Arg70, Arg149, Arg173, Arg186, and Arg189) are conserved in homologous positions, and thus, may interact similarly in different active MAP kinases [46]. The available structural data revealed that most of the small molecule inhibitors of protein kinases bind in the ATP binding pocket [46,47]. ATP binding sites of p38 are the residues corresponding to Glu71, Lys53, and Asp168 [48]. Several of these kinases play a crucial role in cellular proliferation and migration. A scratch assay or cell migration assay was performed in MDA-MB-231 cells with and without cyclosaplin treatment (10 µg/mL) for 0 and 12 h (Figure S5). Cyclosaplin abrogated cellular migration of MDA-MB-231 cells effectively, suggesting the kinases being down-regulated.
Some of the amino acids of MMP-2 interacted with both the peptide ligands (CTTHWGFTLC and cyclosaplin), forming H-bonds and hydrophobic contacts (Figure 5F and Figure 6F). Similarly, in the case of procaspase 7 (Figure 5H and Figure 6H), caspase 9 and survivin (Figure 5 and Figure 6) a few amino acids shared similar interactions with both the peptide-based ligands (positive control; cyclosaplin).
In caspase 9 (Figure 5I and Figure 6I), Gly276 interacted with both RGDS and cyclosaplin, forming hydrophobic contacts, whereas in surviving, Gln92, Phe13, Phe93, Leu14, Leu96, and Lys15 formed interactions with both the peptide based ligands (Figure 5J and Figure 6J). No common interactions were observed in cases of PTEN, MMP-9, and Procaspase 3 (Figure 5 and Figure 6), whereas the positive control failed to interact with TRAIL. Additionally, the apoptotic mechanism was further validated by analyzing the alterations in the mRNA expression of Bax and caspase 3, using GAPDH as control. An upregulation in the mRNA expression of caspase 3 occurred after cyclosaplin treatment in MDA-MB-231 cells (Figure S6a). Similarly, the level of Bax, proapoptotic gene increased in cyclosaplin induced MDA-MB-231 cells after 24 h of treatment (Figure S6a). There was an approximately two-fold increase in mRNA levels of caspase 3 and Bax (Figure S6b). In addition to this, we showed in an earlier study, the induction of apoptosis in MDA-MB-231 cells through caspase-3 activity [11]. Arg1 and Arg8 of cyclosaplin interacted well with the amino acid residues of cancer-related proteins. This could have been possible because Arg side chains provide positive charges as well as hydrogen bonding capabilities to attract the peptide to the negative surface charges of the protein. The binding affinity of cyclic peptide TYY, along with its interaction with EphA4 receptor tyrosine kinase by using AutoDock 4 and LigPlot have been reported elsewhere [49].
Previously, in a 3D cell culture, the efficacy of cyclosaplin was shown as an anticancer agent [50]. Apart from anticancer activity, the other biological activities, such as antimicrobial activity, antiviral activity, and immunomodulatory function, need to be investigated for cyclosaplin. In this context, several analogs of cyclosaplin should be designed and screened in silico for the above-mentioned biological activities prior to in vitro studies. Thus, the in silico experiments gave a clear insight of cyclosaplin potential as an apoptosis inducer and a potential protein kinase inhibitor.
5. Conclusions
The structure of the cyclic octapeptide was elucidated previously by molecular modeling associated with dynamics and was used in the docking studies. Docking studies showed the strong affinity of cyclosaplin towards cancer-related proteins, especially protein kinases and apoptosis-related proteins. Thus, the in silico analyses revealed the potential of cyclosaplin as an apoptosis inducer and a protein kinase inhibitor. Based on these studies, appropriate in vitro and in vivo experiments should be designed rationally to validate its biological activity.
Supplementary Materials
The supplementary materials are available online at https://www.mdpi.com/2218-273X/9/11/740/s1.
Author Contributions
Conceptualization, A.M. and S.D.; Methodology, A.M..; Software, S.D.; Validation, A.M., S.D. and Formal Analysis, A.M.; Investigation, A.M. and S.D.; Resources, S.D.; Writing-Original Draft Preparation, A.M.; Writing-Review and Editing, S.D.; Visualization, A.M.; Supervision, S.D.; Project Administration, S.D.
Funding
This research received no external funding.
Acknowledgments
The Indian Institute of Technology Kharagpur supported the work and the authors would like to acknowledge the Central Research Facility, IIT Kharagpur for providing instrumental support. I am grateful to Saswat Dash for editing and proof-reading the manuscript.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Crowell, J.A. The chemopreventive agent development research program in the Division of Cancer Prevention of the US National Cancer Institute: An overview. Eur. J. Cancer 2005, 41, 1889–1910. [Google Scholar] [CrossRef] [PubMed]
- DiMasi, J.A.; Hansen, R.W.; Grabowski, H.G. The price of innovation: New estimates of drug development costs. J. Health Econ. 2002, 22, 151–185. [Google Scholar] [CrossRef]
- Guariento, S.; Tonelli, M.; Espinoza, S.; Gerasimov, A.S.; Gainetdinov, R.R.; Cichero, E. Rational design, chemical synthesis and biological evaluation of novel biguanides exploring species-specificity responsiveness of TAAR1 agonists. Eur. J. Med. Chem. 2018, 146, 171–184. [Google Scholar] [CrossRef] [PubMed]
- Franchini, S.; Manasieva, L.I.; Sorbi, C.; Battisti, U.M.; Fossa, P.; Cichero, E.; Denora, N.; Iacobazzi, R.M.; Cilia, A.; Pirona, L.; et al. Synthesis, biological evaluation and molecular modeling of 1-oxa-4-thiaspiro-and 1,4-dithiaspiro[4.5]decane derivatives as potent and selective 5-HT1A receptor agonists. Eur. J. Med. Chem. 2017, 125, 435–452. [Google Scholar] [CrossRef]
- Dileep, K.V.; Kelly, M.; Hardin, E. Approaches in the Chemoprevention of Breast Cancer. Cancer Sci. 2013, 5, 1948–5956. [Google Scholar]
- Agyei, D.; Danquah, M.K. Industrial-scale manufacturing of pharmaceuticalgrade bioactive peptides. Biotechnol. Adv. 2011, 29, 272–277. [Google Scholar] [CrossRef]
- Joo, S.H. Cyclic peptides as therapeutic agents and biochemical tools. Biomolecules 2012, 20, 19–26. [Google Scholar] [CrossRef]
- Gurrath, M.; Müller, G.; Kessler, H.; Aumailley, M.; Timpl, R. Conformation/activity studies of rationally designed potent anti-adhesive RGD peptides. Eur. J. Biochem. 1992, 210, 911–921. [Google Scholar] [CrossRef]
- Hotchkiss, R.D. The chemical nature of gramicidin and tyrocidine. J. Biol. Chem. 1941, 141, 171–185. [Google Scholar]
- Vlieghe, P.; Lisowski, V.; Martinez, J.; Khrestchatisky, M. Synthetic therapeutic peptides: Science and market. Drug Discov. Today 2010, 15, 40–56. [Google Scholar] [CrossRef]
- Mishra, A.; Gauri, S.S.; Mukhopadhyay, S.K.; Chatterjee, S.; Das, S.S.; Mandal, S.M.; Dey, S. Identification and structural characterization of a new pro-apoptotic cyclic octapeptide cyclosaplin from somatic seedling of Santalum album L. Peptides 2014, 54, 148–158. [Google Scholar] [CrossRef] [PubMed]
- Gilanski, K.; Rychewlski, L. Protein structure prediction of CAPS5 comparative modelling and fold recognition _targets using consensus alignment and approach and 3D assessment. Proteins 2003, 53, 410–417. [Google Scholar]
- Bikadi, Z.; Hazai, E. Application of the PM6 semi-empirical method to modeling proteins enhances docking accuracy of AutoDock. J. Cheminform. 2009, 1, 1–16. [Google Scholar] [CrossRef] [PubMed]
- Vila, M.C.; Giodarno, R.J.; Sidman, R.L.; Bronk, L.F. From combinatorial peptide selection to drug prototype (II): _targeting the epidermal growth factor receptor pathway. Proc. Nat. Acad. Sci. USA 2010, 107, 5118–5123. [Google Scholar] [CrossRef]
- Alghisi, G.; Ponsonnet, L.; Rüegg, C. The integrin antagonist cilengitide activates αVβ3, disrupts VE-cadherin localization at cell junctions and enhances permeability in endothelial cells. Public Libr. Sci. One 2009, 4, e4449. [Google Scholar] [CrossRef] [PubMed]
- Tal-Gan, Y.; Hurevich, M.; Klein, S.; Ben-Shimon, A.; Rosenthal, D.; Hazan, C.; Shalev, D.E.; Niv, M.Y.; Levitzki, A.; Gilon, C. Backbone cyclic peptide inhibitors of protein kinase B (PKB/Akt). J. Med. Chem. 2011, 54, 5154–5164. [Google Scholar] [CrossRef] [PubMed]
- Gill, K.; Nigam, L.; Singh, R.; Kumar, S.; Dey, S. The Rational Design of Specific Peptide Inhibitor against p38α MAPK at Allosteric-Site: A Therapeutic Modality for HNSCC. PLoS ONE 2014, 9, e101525. [Google Scholar] [CrossRef]
- Zhu, Z.; Jia, J.; Lu, R.; Lu, Y.; Fu, Z.; Zhao, L.; Wang, L.; Jin, M.; Zhao, L.; Gao, W.; et al. Expression of PTEN, p27, p21 and AKT mRNA and protein in human BEL-7402 hepatocarcinoma cells in transplanted tumors of nude mice treated with the tripeptide tyroservatide (YSV). Int. J. Cancer 2006, 118, 1539–1544. [Google Scholar] [CrossRef]
- Koivunen, E.; Arap, W.; Valtanen, H.; Rainisalo, A.; Medina, O.P.; Heikkilä, P.; Sorsa, T. Tumor _targeting with a selective gelatinase inhibitor. Nat. Biotechnol. 1999, 17, 768–774. [Google Scholar] [CrossRef]
- Aguzzi, M.S.; Giampeitri, C.; Marchis, F.D. RGDS peptide induces caspase 8 and caspase 9 activation in human endothelial cells. Blood 2004, 103, 4180–4187. [Google Scholar] [CrossRef]
- Heins, M.S.; Quax, W.J. Implications of a Newly Discovered DR5 Specific Antagonistic Peptide for Neurodegenerative Disorders. Mol. Cell. Pharmacol. 2010, 277, 1653–1665. [Google Scholar]
- Hirohashi, Y.; Torigoe, T.; Maeda, A.; Nabeta, Y. An HLA-A24-restricted Cytotoxic T Lymphocyte Epitope of a Tumor-associated Protein, Survivin. Clin. Cancer Res. 2002, 8, 1731. [Google Scholar] [PubMed]
- Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 2001, 46, 3–26. [Google Scholar] [CrossRef]
- Gfeller, D.; Grosdidier, A.; Wirth, M. Swiss _target Prediction: A web server for _target prediction of bioactive small molecules. Nucleic Acid Res. 2013, 42, 32–38. [Google Scholar] [CrossRef]
- Berman, H.M.; Henrick, K.; Nakamura, H. Announcing the world wide Protein Data Bank. Nat. Struct. Biol. 2000, 10, 980. [Google Scholar]
- Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef]
- Wallace, A.C.; Laskowski, R.A.; Thornton, J.M. LigPlot: A program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995, 8, 127–134. [Google Scholar] [CrossRef]
- Roxin, A.; Zheng, G. Flexible or fixed: A comparative review of linear and cyclic cancer-_targeting peptides. Future Med. Chem. 2012, 4, 1601–1618. [Google Scholar] [CrossRef]
- KCathcart, J.; Pulkoski-Gross, A.; Cao, J. _targeting matrix metalloproteinases in cancer: Bringing new life to old ideas. Genes Dis. 2015, 2, 26–34. [Google Scholar] [CrossRef]
- Koskimaki, J.E.; Pandey, N.B.; Tamiz, A.P.; Popel, A.S. Structure-activity relationship study of collagen-derived anti-angiogenic biomimetic peptides. Chem. Biol. Drug Des. 2012, 1, 27–37. [Google Scholar]
- Pero, S.C.; Shukla, G.S.; Armstrong, A.L.; Peterson, D.; Fuller, S.P.; Godin, K.; Kingsley-Richards, S.L.; Weaver, D.L.; Bond, J.; Krag, D.N. Identification of a small peptide that inhibits the phosphorylation of ErbB2 and proliferation of ErbB2 overexpressing breast cancer cells. Int. J. Cancer 2004, 111, 951–960. [Google Scholar] [CrossRef] [PubMed]
- Ubol, S.; Kramyu, J.; Masrinoul, P.; Kachangchaeng, C.; Pittayanurak, P.; Sophasan, S.; Reutrakul, V. A novel cycloheptapeptide exerts strong anticancer activity via stimulation of multiple apoptotic pathways in Caspase-3 deficient cancer cells. Anticancer Res. 2007, 27, 2473–2480. [Google Scholar] [PubMed]
- Janin, J.; Henrick, K.; Moult, J. A critical assessment of predicted interaction. Proteins Struct. Funct. Genet. 2003, 52, 2–9. [Google Scholar] [CrossRef] [PubMed]
- Kaplan, N.; Morpurgo, N.; Linial, M. Novel families of toxin-like peptides in insects and mammals: A computational approach. J. Mol. Biol. 2007, 369, 553–566. [Google Scholar] [CrossRef] [PubMed]
- Mader, J.S.; Salsman, J.; Conrad, D.M.; Hoskin, D.W. Bovine lactoferricin inhibits basic fibroblast growth factor and vascular endothelial growth factor 165-induced angiogenesis by competing for heparin-like binding sites on endothelial sites. Am. J. Pathol. 2006, 169, 1753–1766. [Google Scholar] [CrossRef]
- Ruegg, C. Vascular integrins in tumor angiogenesis: Mediators and therapeutic _targets. Endothelium 2006, 13, 113–135. [Google Scholar]
- DeLorbe, J.E.; Clements, J.H.; Teresk, M.G.; Befield, A.P.; Plake, H.R.; Millspaugh, L.E.; Martin, S.F. Thermodynamic and structural effects of conformational constraints in protein-ligand interactions: Entropic paradoxy associated with ligand pre-organization. J. Am. Chem. Soc. 2009, 131, 16758–16770. [Google Scholar] [CrossRef]
- Yang, K.S.; Macdonald-Obermann, J.L.; Linda, J.P. Asp-960/Glu-961 controls the movement of the c-terminal tail of the epidermal growth factor receptor to regulate asymmetric dimer formation. J. Biol. Chem. 2010, 285, 24014–24022. [Google Scholar] [CrossRef] [Green Version]
- Ward, C.W.; Gough, K.H.; Rashke, M.; Wan, S.S.; Tribbick, G.; Wang, J. Systematic mapping of potential binding sites for Shc and Grb2 SH2 domains on insulin receptor substrate-1 and the receptors for insulin, epidermal growth factor, platelet-derived growth factor, and fibroblast growth factor. J. Biol. Chem. 1996, 271, 5603–5609. [Google Scholar] [CrossRef] [Green Version]
- Stamos, J.; Swilkowski, M.X.; Eigenbrot, C. Structure of the epidermal growth factor receptor kinase domain alone and in complex with a 4 anilinoquinazoline inhibitor. J. Biol. Chem. 2002, 277, 46265–46272. [Google Scholar] [CrossRef] [Green Version]
- Kyte, J.; Doolittle, R.F. A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 1982, 157, 105–132. [Google Scholar] [CrossRef] [Green Version]
- Wang, Z.; Wang, N.; Han, S.; Wang, D.; Mo, S.; Yu, L.; Huang, H.; Tsui, K.; Shen, J.; Chen, J. Dietary compound Isoliquiritigenin inhibits breast cancer neoangiogenesis via VEGF/VEGFR-2 Signaling Pathway. PLoS ONE 2013, 8, e68566. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Olsson, A.K.; Dimberg, A.; Kreuger, J.; Claesson-Welsh, L. VEGF receptor signalling-in control of vascular function. Nat. Rev. Mol. Cell Biol. 2006, 7, 359–371. [Google Scholar] [CrossRef] [PubMed]
- Chan, T.O.; Zhang, J.; Rodeck, U.; Pascal, J.M.; Armen, R.S.; Spring, M.; Dumitru, C.D.; Myers, V.; Li, X.; Cheung, J.Y.; et al. Resistance of Akt-kinases to dephosphorylation through ATP-dependent conformational plasticity. Proc. Natl. Acad. Sci. USA 2011, 108, e1120–e1127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Yang, J.; Cron, P.; Good, V.M.; Thompson, V.; Hemmings, B.A.; Barford, D. Crystal structure of an activated Akt/protein kinase B ternary complex with GSK3-peptide, A.M.P.P.N.P. Nat. J. Struct. Biol. 2002, 12, 940–944. [Google Scholar] [CrossRef] [PubMed]
- Wang, Z.; Harkins, P.C.; Ulevitch, R.J.; Han, J.; Cobb, M.H.; Goldsmith, E.J. The structure of mitogen-activated protein kinase p38 at 2.1-Å resolution. Proc. Natl. Acad. Sci. USA 1997, 94, 2327–2332. [Google Scholar]
- De Azevedo, W.F.J.; Dieckmann, H.J.M.; Gahmen, U.S.; Worland, P.J.; Sausville, E.; Kim, S.H. Structural basis for specificity and potency of a flavonoid inhibitor of human CDK2, a cell cycle kinase. Proc. Natl. Acad. Sci. USA 1996, 83, 2735–2740. [Google Scholar] [CrossRef] [Green Version]
- Wilson, K.P.; Fitzgibbon, M.J.; Caron, P.N.; Griffith, J.P.; Chen, W.Y.; McCaffrey, P.G.; Chambers, S.P.; Su, M.S.S. Crystal structure of p38 mitogen-activated protein kinase. J. Biol. Chem. 1996, 271, 27696–27700. [Google Scholar] [CrossRef] [Green Version]
- Han, X.; Xu, Y.; Yang, Y.; Xi, J.; Tian, W. Discovery and characterization of a novel cyclic peptide that effectively inhibits ephrin binding to the EphA4 receptor and displays anti-angiogenesis activity. PLoS ONE 2013, 8, e80183. [Google Scholar] [CrossRef] [Green Version]
- Mishra, A.; Mukhopadhyay, S.K.; Dey, S. Evaluation of Cyclosaplin Efficacy Using a Silk Based 3D Tumor Model. Biomolecules 2019, 9, 123. [Google Scholar] [CrossRef] [Green Version]
Figure 1.
The different peptide-based ligands for _targeting cancer-related proteins used in docking studies. Cyan blue = carbon, grey = hydrogen, deep blue = nitrogen, red = oxygen, and yellow = sulfur.
Figure 1.
The different peptide-based ligands for _targeting cancer-related proteins used in docking studies. Cyan blue = carbon, grey = hydrogen, deep blue = nitrogen, red = oxygen, and yellow = sulfur.
Figure 2.
(a) Possible _targets of cyclosaplin as predicted by Swiss _target Prediction. (b) Different energy minimized proteins (rainbow spectrum) used in docking studies. (A) EGFR kinase. (B) VEGFR2 kinase. (C) PKB. (D) p38. (E) PTEN. (F) MMP-2. (G) MMP-9. (H) Procaspase 3. (I) Procaspase 7. (J) Caspase 9. (K) TRAIL. (L) SURVIVIN.
Figure 2.
(a) Possible _targets of cyclosaplin as predicted by Swiss _target Prediction. (b) Different energy minimized proteins (rainbow spectrum) used in docking studies. (A) EGFR kinase. (B) VEGFR2 kinase. (C) PKB. (D) p38. (E) PTEN. (F) MMP-2. (G) MMP-9. (H) Procaspase 3. (I) Procaspase 7. (J) Caspase 9. (K) TRAIL. (L) SURVIVIN.
Figure 3.
(a) Docking scores in kcal/mol for various cancer-related proteins. The binding affinities closer to 10 indicate efficient binding. (b) Cyclosaplin bound to different receptors. Cyclosaplin is shown in white, indicated by the arrows, and proteins are depicted with rainbow’s spectrum. (A) EGFR kinase. (B) VEGFR2 kinase. (C) PKB. (D) p38. (E) PTEN. (F) MMP-2. (G) MMP-9. (H) Procaspase 3 (previous study [11]). (I) Procaspase 7. (J) Caspase 9. (K) TRAIL. (L) SURVIVIN.
Figure 3.
(a) Docking scores in kcal/mol for various cancer-related proteins. The binding affinities closer to 10 indicate efficient binding. (b) Cyclosaplin bound to different receptors. Cyclosaplin is shown in white, indicated by the arrows, and proteins are depicted with rainbow’s spectrum. (A) EGFR kinase. (B) VEGFR2 kinase. (C) PKB. (D) p38. (E) PTEN. (F) MMP-2. (G) MMP-9. (H) Procaspase 3 (previous study [11]). (I) Procaspase 7. (J) Caspase 9. (K) TRAIL. (L) SURVIVIN.
Figure 4.
Peptide based ligands bound to specific proteins. Ligands are shown in white, indicated by the arrows, and proteins are depicted in rainbow spectrum. (A) CVRACGAD bound to EGFR kinase, (B) Cilengitide bound to VEGFR2 kinase, (C) RPRTSSF bound to PKB, (D) FWCS bound to p38, (E) YSV bound to PTEN, (F) CTTHWGFTLC bound to MMP-2, (G) CRRHWGFEFC bound to MMP-9, (H) Cilengitide bound to Procaspase 3, (I) RGDS bound to Procaspase 7, (J) RGDS bound to Caspase 9, (K) CKVILTHRC unbound to TRAIL, and (L) AYACNTSTL bound to SURVIVIN.
Figure 4.
Peptide based ligands bound to specific proteins. Ligands are shown in white, indicated by the arrows, and proteins are depicted in rainbow spectrum. (A) CVRACGAD bound to EGFR kinase, (B) Cilengitide bound to VEGFR2 kinase, (C) RPRTSSF bound to PKB, (D) FWCS bound to p38, (E) YSV bound to PTEN, (F) CTTHWGFTLC bound to MMP-2, (G) CRRHWGFEFC bound to MMP-9, (H) Cilengitide bound to Procaspase 3, (I) RGDS bound to Procaspase 7, (J) RGDS bound to Caspase 9, (K) CKVILTHRC unbound to TRAIL, and (L) AYACNTSTL bound to SURVIVIN.
Figure 5.
Interaction of cyclosaplin with various cancer-related proteins. (A) EGFR Kinase. (B) VEGFR2 Kinase. (C) PKB. (D) p38. (E) PTEN. (F) MMP-2. (G) Procaspase 3. (H) Procaspase 7. (I) Caspase 9. (J) SURVIVIN.
Figure 5.
Interaction of cyclosaplin with various cancer-related proteins. (A) EGFR Kinase. (B) VEGFR2 Kinase. (C) PKB. (D) p38. (E) PTEN. (F) MMP-2. (G) Procaspase 3. (H) Procaspase 7. (I) Caspase 9. (J) SURVIVIN.
Figure 6.
Protein–ligand interactions using LigPlot. (A) CVRACGAD and EGFR kinase. (B) Cilengitide and VEGFR2 kinase. (C) RPRTSSF and PKB. (D) FWCS and p38. (E) YSV and PTEN. (F) CTTHWGFTLC and MMP-2. (G) Cilengitide and Procaspase 3. (H) RGDS and Procaspase 7. (I) RGDS and Caspase 9. (J) AYACNTSTL and SURVIVIN.
Figure 6.
Protein–ligand interactions using LigPlot. (A) CVRACGAD and EGFR kinase. (B) Cilengitide and VEGFR2 kinase. (C) RPRTSSF and PKB. (D) FWCS and p38. (E) YSV and PTEN. (F) CTTHWGFTLC and MMP-2. (G) Cilengitide and Procaspase 3. (H) RGDS and Procaspase 7. (I) RGDS and Caspase 9. (J) AYACNTSTL and SURVIVIN.
Table 1.
Ligands used in the study.
| S.No. | Ligand | References |
|---|---|---|
| 1 | CVRACGAD (Cyclic) | [14] |
| 2 | Cilengitide (Cyclic) | [15] |
| 3 | RPRTSSF (Cyclic) | [16] |
| 4 | FWCS (Linear) | [17] |
| 5 | YSV (Linear) | [18] |
| 6 | CTTHWGFTLC (Cyclic) | [19] |
| 7 | CRRHWGFEFC (Cyclic) | [19] |
| 8 | RGDS (Linear) | [20] |
| 9 | CKVILTHRC (Cyclic) | [21] |
| 10 | AYACNTSTL (Linear) | [22] |
| 11 | Cyclosaplin (Cyclic) | [11] |
Table 2.
Molecular weight and molecular formulae of the ligands.
| S.No. | Ligand | Molecular Weight (Da) | Molecular Formula |
|---|---|---|---|
| 1 | CVRACGAD (Cyclic) | 791.9 | C29H49N11O11S2 |
| 2 | Cilengitide (Cyclic) | 588.6 | C27H40N8O7 |
| 3 | RPRTSSF (Cyclic) | 875.0 | C39H66N14O9 |
| 4 | FWCS (Linear) | 541.6 | C26H31N5O6S1 |
| 5 | YSV (Linear) | 367.4 | C17H25N3O6 |
| 6 | CTTHWGFTLC (Cyclic) | 1166.3 | C52H71N13O14S2 |
| 7 | CRRHWGFEFC (Cyclic) | 1338.5 | C60H79N19O13S2 |
| 8 | RGDS (Linear) | 433.4 | C15H27N7O8 |
| 9 | CKVILTHRC (Cyclic) | 1070.3 | C45H79N15O11S2 |
| 10 | AYACNTSTL (Linear) | 943.0 | C39H62N10O15S1 |
| 11 | Cyclosaplin (Cyclic) | 858.9 | C33H60N14O12S1 |
Table 3.
Physiochemical parameters of ligand molecules screened for Lipinski’s rule.
| Ligand | Molecular Weight (Da) | Hydrogen Bond Donor | Hydrogen Bond Acceptor | LogP | Molar Refractivity | Rules Satisfied |
|---|---|---|---|---|---|---|
| CVRACGAD | 791.9 | 13 | 13 | −4.7 | 206.6 | 1/5 |
| Cilengitide | 588.6 | 7 | 8 | −1.4 | 170.9 | 2/5 |
| RPRTSSF | 875.0 | 15 | 12 | −5.9 | 236.7 | 0/5 |
| FWCS | 541.6 | 8 | 7 | −0.7 | 143.2 | 2/5 |
| YSV | 367.4 | 6 | 6 | −1.0 | 93.14 | 4/5 |
| CTTHWGFTLC | 1166.3 | 16 | 17 | −4.0 | 331.1 | 1/5 |
| CRRHWGFEFC | 1338.5 | 20 | 17 | −3.6 | 381.6 | 1/5 |
| RGDS | 433.4 | 10 | 8 | −4.7 | 99.7 | 4/5 |
| CKVILTHRC | 1070.3 | 16 | 16 | −3.3 | 306.4 | 1/5 |
| AYACNTSTL | 943.0 | 16 | 16 | −6.1 | 230.0 | 0/5 |
| Cyclosaplin | 858.9 | 17 | 13 | −6.5 | 243.0 | 0/5 |
Table 4.
Comparative binding affinity of different ligands with receptors.
| S.No. | Receptor | Ligand | Binding Affinity (kcal/mol) |
|---|---|---|---|
| 1 | Epidermal Growth Factor Receptor Kinase | CVRACGAD Cyclosaplin | −7.7 −6.8 |
| 2 | Vascular Endothelial Growth Factor r 2 Receptor Kinase | Cilengitide Cyclosaplin | −8.1 −7.8 |
| 3 | Protein Kinase B | RPRTSSF Cyclosaplin | −7.5 −8.1 |
| 4 | p38 (Mitogen Activated Protein Kinase) | FWCS Cyclosaplin | −8.9 −8.3 |
| 5 | PTEN | YSV Cyclosaplin | −7.8 6.3 |
| 6 | Matrix metalloproteinase-2 (MMP-2) | CTTHWGFTLC Cyclosaplin | −7.8 −8.2 |
| 7 | Matrix metalloproteinase-9 (MMP-9) | CRRHWGFEFC Cyclosaplin | −8.4 −7.3 |
| 8 | Procaspase 3 | Cilengitide Cyclosaplin | −8.1 −7.8 |
| 9 | Procaspase 7 | RGDS Cyclosaplin | −6.8 −8.7 |
| 10 | Caspase 9 | RGDS Cyclosaplin | −6.7 −8.9 |
| 11 | TRAIL | CKVILTHRC Cyclosaplin | −6.4 −8.2 |
| 12 | SURVIVIN | AYACNTSTL Cyclosaplin | −7.2 −7.4 |
Table 5.
Molecular interactions of ligands with amino acids of proteins (amino acids showing similar interactions are marked in bold; black = cis similarity; red = trans similarity).
Table 5.
Molecular interactions of ligands with amino acids of proteins (amino acids showing similar interactions are marked in bold; black = cis similarity; red = trans similarity).
| S.No. | Protein | Ligand | Hydrophilic Interactions | Hydrophobic Contacts | No. of H-Bonds |
|---|---|---|---|---|---|
| 1 | EGFR Kinase | CVRACGAD | Lys721, Thr766 | Ala719, Asp831, Gly695, Gly772, Leu694, Leu768 | 2 |
| Cyclosaplin | Glu961, Asp960 (4) | Arg962, Asp950, Gln788, Gln952, Gly786, Met963, Pro951, Ser787, Tyr789 | 5 | ||
| 2 | VEGFR 2 Kinase | Cilengitide | Asn923, Asp 1046, Thr926 | Ala866, Arg1032, Arg1066, Asp1064, Cys919, Cys1045, Gly841, Leu840, Phe918, Ser925, Val848 | 3 |
| Cyclosaplin | Arg842, Arg842, Glu885 | Ala1065, Asn923, Asp1046, Asp1064, Gly843, Gly846, Gln847, Leu840, Leu1035, Lys868, Thr926, Val848 | 3 | ||
| 3 | Protein Kinase B | RPRTSSF | Arg274, Glu315, Lys181, Lys277, Tyr327 | Asp275, Glu200, Gly295, Leu317, Lys160, Phe163, Thr162, Thr199, Val198 | 5 |
| Cyclosaplin | Arg274, Arg274, Asp275, Asp293, Val272, Val272 | Ile188, Leu183, Lys181, Phe163, Thr199, Tyr273 | 6 | ||
| 4 | p38 (Mitogen Activated Protein Kinase) | FWCS | Arg173, Arg 67 | Asp168, Glu71, Glu178, His64, Leu74, Leu75, Leu171, Ly53, Phe169, Thr68, Tyr35 | 2 |
| Cyclosaplin | Asp112, Asp112, Asp150, Asn115, Lys53, Phe169, Phe169, Ser154 | Arg173, Asn114, Asn155, Asp168, Gly170, Leu167, Met109, Tyr35, Val38 | 8 | ||
| 5 | PTEN | YSV | Arg172, Arg173 | Asp324, Leu318, Leu320, Phe279, Tyr176, Tyr177, Tyr180, Val275 | 2 |
| Cyclosaplin | Ala72, Ala72, Gln87, Gln97 | Glu91, Glu99, Leu100, Pro89, Tyr88 | 4 | ||
| 6 | MMP-2 | CTTHWGFTLC | Arg98, Gly371, Gly394, Thr 511, Thr547 | Gln393, Gly216, Lys99, Lys372, Met373, Ser365, Ser546, Tyr395, Tyr425, Tyr427 | 6 |
| Cyclosaplin | Gly394, Tyr425 | Asp392, Glu515, Gln393, Gly216, Gly371, Phe512, Pro100, Pro514, Leu 548, Ser546, Thr426, Tyr427, Thr511, Tyr277, Tyr395 | 2 | ||
| 7 | MMP-9 | CRRHWGFEFC | Leu371, Arg2, Cys1 | Arg370, Arg424, Glu427, Gln391, Gly392, Lys92, Phe425, Pro97, Pro233, Ser240, Ser242, Thr426 Tyr393, Tyr423 | 3 |
| Cyclosaplin | Arg221, Thr331 | Arg279, Asp226, Asp284, Gly227, Gly285, Pro219, Pro272, Thr220, | 4 | ||
| 8 | Procaspase 3 | Cilengitide | Ala33, Arg238, Asn32, Ser36 (2), Tyr37, Tyr274, Tyr276 | Asn35, Glu272, Leu230, Lys38, His234 | 8 |
| Cyclosaplin | Gln261, Glu124, Lys186 | Arg164, Gly125, Ile126, Ile187, Leu136, Lys186, Pro188, Tyr197, Val189 | 3 | ||
| 9 | Procaspase 7 | RGDS | Arg87, Asn88, Gly228, Gln184 Thr189, Tyr229 | Arg187, Gly188, Gln287, His144, Lys285, Pro227, Ser239, Thr90 | 6 |
| Cyclosaplin | Arg170, Glu176, Glu284, Gly168, Phe282, Phe282, Ser277 | Ala169, Ala217, Arg167, Asp204, Gln276, Lys286, Leu175, Gln287, Ile288, Val215, Glu216, His283 | 7 | ||
| 10 | Caspase 9 | RGDS | Gly269, Ser339 | Ala149, Arg408, Asp150, Asp340, Gly276, Gly277, Gln399, Ile154, Ile396, Lys398, Met400, Thr337 | 2 |
| Cyclosaplin | Pro273, Ser144, Ser144 | Arg146, Asp228, Glu143, Gly147, Gly225, Gly276, Ile154, Leu155, Lys278, Lys414, Ser274 | 3 | ||
| 11 | SURVIVIN | AYACNTSTL | Arg18, Arg37, Cys31, Gly30 | Glu29, Glu36, Glu40, Gln92, Ile74, Leu14, Leu96, Leu104, Lys15, Lys90, Met1, Phe13, Phe93, Thr34 | 4 |
| Cyclosaplin | Gln92, Glu94, Lys91, Phe13 | Asp16, Gly2, Leu14, Leu96, Leu102, Lys15, Phe 93, Pro4 | 4 |
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Mishra, A.; Dey, S. Molecular Docking Studies of a Cyclic Octapeptide-Cyclosaplin from Sandalwood. Biomolecules 2019, 9, 740. https://doi.org/10.3390/biom9110740
AMA Style
Mishra A, Dey S. Molecular Docking Studies of a Cyclic Octapeptide-Cyclosaplin from Sandalwood. Biomolecules. 2019; 9(11):740. https://doi.org/10.3390/biom9110740
Chicago/Turabian StyleMishra, Abheepsa, and Satyahari Dey. 2019. "Molecular Docking Studies of a Cyclic Octapeptide-Cyclosaplin from Sandalwood" Biomolecules 9, no. 11: 740. https://doi.org/10.3390/biom9110740
APA StyleMishra, A., & Dey, S. (2019). Molecular Docking Studies of a Cyclic Octapeptide-Cyclosaplin from Sandalwood. Biomolecules, 9(11), 740. https://doi.org/10.3390/biom9110740
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
