Next Article in Journal
Variation of the Main Alkaloid Content in Equisetum palustre L. in the Light of Its Ontogeny
Previous Article in Journal
Effects of Cassava Juice (Manihot esculenta Crantz) on Renal and Hepatic Function and Motor Impairments in Male Rats
Previous Article in Special Issue
Isolation of an Anti–tumour Disintegrin: Dabmaurin–1, a Peptide Lebein–1–like, from Daboia mauritanica Venom
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 12
Issue 11
10.3390/toxins12110709
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Structural Insight into Integrin Recognition and Anticancer Activity of Echistatin

by
Yi-Chun Chen
†,
Yao-Tsung Chang
†,
Chiu-Yueh Chen
,
Jia-Hau Shiu
,
Chun-Ho Cheng
,
Chun-Hao Huang
,
Ju-Fei Chen
and
Woei-Jer Chuang
*
Department of Biochemistry and Molecular Biology, National Cheng Kung University College of Medicine, Tainan 70101, Taiwan
*
Author to whom correspondence should be addressed.
†
These authors contributed to this work equally.
Toxins 2020, 12(11), 709; https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12110709
Submission received: 20 October 2020 / Revised: 1 November 2020 / Accepted: 6 November 2020 / Published: 9 November 2020
(This article belongs to the Special Issue Medically Relevant Snake Toxins, Current Antivenoms and the Development of Next-Generation Antivenoms)
Browse Figures
Review Reports Versions Notes

Abstract

:
Echistatin (Ech) is a short disintegrin with a long 42NPHKGPAT C-terminal tail. We determined the 3-D structure of Ech by X-ray crystallography. Superimposition of the structures of chains A and B showed conformational differences in their RGD loops and C-termini. The chain A structure is consistent with our NMR analysis that the GPAT residues of the C-terminus cannot be observed due to high flexibility. The hydrogen bond patterns of the RGD loop and between the RGD loop and C-terminus in Ech were the same as those of the corresponding residues in medium disintegrins. The mutant with C-terminal HKGPAT truncation caused 6.4-, 7.0-, 11.7-, and 18.6-fold decreases in inhibiting integrins αvβ3, αIIbβ3, αvβ5, and α5β1. Mutagenesis of the C-terminus showed that the H44A mutant caused 2.5- and 4.4-fold increases in inhibiting αIIbβ3 and α5β1, and the K45A mutant caused a 2.6-fold decrease in inhibiting αIIbβ3. We found that Ech inhibited VEGF-induced HUVEC proliferation with an IC50 value of 103.2 nM and inhibited the migration of A375, U373MG, and Panc-1 tumor cells with IC50 values of 1.5, 5.7, and 154.5 nM. These findings suggest that Ech is a potential anticancer agent, and its C-terminal region can be optimized to improve its anticancer activity.
Keywords:
echistatin; disintegrins; crystal structure; cancer cell migration; HUVEC proliferation; anticancer drug design
Key Contribution: Disintegrins are cysteine-rich venom proteins that have received attention due to their potential use for cancer treatment. This is the first study to report the crystal structure of Ech, a short disintegrin. We found that Ech has a unique C-terminal tail structure that can be used to improve its anticancer activity.

1. Introduction

Echistatin (Ech) is a short disintegrin isolated from Echis carinatus that contains 49 residues with 4 disulfide bonds. Ech is the most extensively characterized disintegrin and has a conserved RGD motif and a long C-terminal NPHKGPAT tail [1,2,3,4,5,6]. Several studies support a functional role for the C-terminal residues of Ech in modulating its binding affinity and selectivity towards integrin αIIbβ3 [6,7,8]. For example, a synthetic C-terminal peptide of Ech inhibited the binding of integrin αIIbβ3 to Ech, activated the binding of integrin to immobilized ligands, induced the expression of activation-dependent conformational epitopes of integrin, and increased the binding of fibrinogen to peptide-treated platelets [6,7,9,10]. An RGD peptide fused with the C-terminal KGPAT residues of Ech selectively binds αvβ3 integrin but not αvβ5 integrin [11]. Replacement of contortrostatin’s C-terminal FHA with Ech’s C-terminal HKGPAT tail results in a higher binding affinity to integrin α5β1 [12]. In our previous study, we also demonstrated that the length of the C-terminus and the K45 residue of Ech are important for its interactions with integrin αIIbβ3 [13]. These results show that the C-terminal tail of Ech may act in synergy with the RGD loop to interact with integrins. However, little is known about the role of the C-terminus in recognizing RGD-dependent integrins.
Many 3-D structures of disintegrins have been determined by X-ray crystallography and NMR spectroscopy [6,14,15,16,17]. Comparison of the 3-D structures of disintegrins reveals a very similar fold with an elongated shape [14,15,16,17,18]. Their secondary structures are composed of a series of short antiparallel β-sheets. The RGD motif is located at the apex of the integrin-binding loop, which is structurally close to the C-terminus. The 3-D structure of Ech has also been determined by NMR spectroscopy [3,6,13,19]. Similar to the reported structures of other disintegrins, Ech has a well-defined core composed of a series of nonclassical turns crosslinked by disulfide bonds (C2-C11, C7-C32, C8-C37, and C20-C39). The RGD motif is located in a protruding loop that is formed by two antiparallel β-strands (residues 16–20 and 29–31). A dynamics study showed that the RGD loop and the C-terminal region exhibited concerted motion [6]. This is consistent with functional studies showing that the C-terminus of Ech plays an important role in integrin recognition [6,7,8,9,10,11,12,13]. However, the C-terminal region of Ech is flexible based on the reported NMR structure [3,6,19]. It is essential to determine the structural conformation of Ech’s C-terminus in a high-resolution crystal structure to elucidate its structure–function relationships.
Integrins have been found to be associated with tumor angiogenesis, tumor cell migration, proliferation, and metastasis [20,21,22,23]. Disintegrins are potent integrin inhibitors [24,25,26,27,28,29,30,31,32,33]. Many in vitro and in vivo studies have demonstrated that RGD-containing disintegrins have antiangiogenic and antimetastatic activities [24,25,26,27,28,29,30,31,32,33]. They are effective in inhibiting angiogenesis, tumor cell adhesion, and tumor cell migration with micromolar to high nanomolar potency [24,25]. Ech also significantly inhibited angiogenesis, lung metastasis of B16F10 melanoma, metastasis of MV3 melanoma, and the hyaluronan-mediated migration of glioblastoma cells in fibrin gels [26]. To evaluate the inhibitory activities of recombinant Ech and its mutants, we determined the IC50 values of VEGF-induced HUVEC proliferation and tumor cell migration, as well as their activities in inhibiting integrins αvβ3, αvβ5, and α5β1.
To analyze how the C-terminus of Ech interacts with the RGD loop to recognize integrins, we determined its 3-D structure using X-ray crystallography. Site-directed mutagenesis and deletion studies on Ech were used to identify the essential residues involved in the recognition of integrins, as well as its antiangiogenic and antimetastatic activities. We found that the major differences in disintegrins were the amino acid composition, length, and conformation of the C-terminus. These structural and functional studies of Ech serve as the basis for insight into their structure–function relationships and for the design of anticancer drugs.

2. Results

2.1. The Overall Structure of Ech

The diffraction data of Ech were obtained to a resolution of 1.80 Å from a single frozen crystal, which contained two molecules per asymmetric unit. Figure 1A–D show a representative electron-density fit including the RGD loop from residues K21 to D30 and the C-terminus from residues P40 to T49. Ramachandran plot analysis of the final structure, obtained from the program MolProbity [34], showed 95.7% of the residues in the favored region, 4.3% in the additionally allowed region, and no residues in the outlier region (Table 1). The data collection and structure refinement statistics are summarized in Table 1. The overall structure of chains A and B contained four disulfide bonds (C2-C11, C7-C32, C8-C37, and C20-C39) and one pair of antiparallel β-strands (residues 16–20 and 29–31) (Figure 2A,B).
Superimposition of the backbones of chains A and B with an RMSD value of 0.226 Å demonstrated that their overall structures were very similar, except for the RGD loop and C-terminal loop (Figure 2C–E). Superimposition of the RGD loop (residues K21-D30) and C-terminus (residues P40-K45) resulted in RMSD values of 0.162 Å and 0.581 Å, respectively (Figure 2D,E). Structural analysis showed that the RGD loop of chain B was squeezed by residues K21, R41, and K45 in contact with sulfate ions and residue S4 of chain A through crystal packing (Figure S1). Therefore, the R and D side chains of the chain B structure interacted with each other and formed hydrogen bonds (Figure 2D). In contrast to chain B, the R and D side chains in the RGD loop of chain A are far from each other (Figure 2D). Furthermore, the backbones of the C-terminal GPAT residues of chain A showed a weaker electron density map (Figure 1C). This was consistent with our unpublished NMR data and NMR nuclear Overhauser effect (NOE) analysis, which showed that the GPAT residues of the C-terminus exhibited low or negative NOE values (Figure S2). The average NOE value of Ech is 0.605. In addition, residues R24-D27 in the RGD loop, E1-E3 in the N-terminus, and K45-T49 in the C-terminus had NOE values less than the average NOE value of Ech, suggesting that these regions of Ech were flexible. In contrast, the conformation of the C-terminal tail of Ech chain B was well defined (Figure 1D), possibly because the C-terminal residue K45 can interact with the O atoms of SO42- (Figure S1, Table S1). Moreover, the C-terminal residue T49 can interact with water molecules (Figure S3).
Although the tertiary fold of Ech was similar to the reported NMR structure (2Ech), superimposition of the backbones of Ech chain A and the NMR structure (2Ech) resulted in RMSD values of 1.688 Å (Figure 2F,G) [13]. Nevertheless, the analyses of their structures showed that the RGD motif is located at the apex of the integrin-binding loop, which is structurally close to the C-terminus. These findings were consistent with the reported 3-D structures of medium disintegrins [14,15,16,17,18]. The resulting high-resolution crystal structure can provide an improved molecular basis to design anticancer agents using Ech as a scaffold.

2.2. Structural Comparison of the RGD Loop and the Interactions between the RGD Loop and C-Termini

Hydrogen bond connectivity and atom distances of the RGD loop were compared among Ech chains A and B and trimestatin, a medium disintegrin with an ARGDFP motif, to identify their structural differences in detail (Figure 3). We found that the Cα to Cα distance between residues R and D in chain A, chain B, and trimestatin of the RGD motif was 6.5, 5.8, and 6.7 Å, and the Cζ to Cγ distance between residues R and D in chain A, chain B, and trimestatin was 12.4, 4.8, and 11.1 Å (Figure 3A–C, Table 2). These results showed that Ech chain A and trimestatin exhibited similarity in the conformation of the RGD motif and in the sidechain orientation of the R and D residues, which were pointed away from each other. Moreover, comparison of the Cα(Ri)-Cα(Di+2), Cβ(Ri)-Cβ(Di+2), Cζ(Ri)-Cγ(Di+2), and Cα(Ri)-Cα(Xi+3) distances (Å) of Ech chain A and Ech’s NMR structure (2Ech) revealed that the RGD motif in our X-ray structure (Ech chain A) was similar to that in the solution NMR structure (2Ech) (Table 2). Interestingly, we found a conserved and extensively connected hydrogen bond network in the RGD loop. The residues located on both sides of the RGD motif in Ech formed hydrogen bonds, including K21 with D30, A23 with D29, and R24 with D27 (Figure 3D,E, Tables S2 and S3). The corresponding residues of trimestatin also formed hydrogen bonds, including R46 with D55 and A48 with D54 (Figure 3F, Table S4). These results showed that short and medium disintegrins have similar RGD motif and loop conformations, suggesting that these conserved interactions may be important to integrin recognition because they maintain the conformation of the RGD loop.
Similar to the conformation of the RGD loop, the interactions between the RGD loop and C-terminus were conserved in short and medium disintegrins (Figure 4A–C, Tables S2–S4). The structural analysis of Ech chains A and B showed that the side chain of the D30 residue in the RGD loop and the amide proton of the N42 residue in the C-terminus formed hydrogen bonds (Figure 4A,B). The corresponding residues of trimestatin also formed hydrogen bonds between the side chain of residue D55 in the RGD loop and the amide proton of residue W67 (Figure 4C). In contrast to Ech chain A and trimestatin, Ech chain B formed additional hydrogen bonds, including residue M28 with residue K45 and residue D30 with residue K45 (Figure 4B). These results implied that the longer C-terminal tail, which forms the loop-like structure, may play an important role in the interaction with integrins.

2.3. Structural Comparison of the C-Terminal Regions

In addition to the similar conformation of the RGD loop and the interactions between the RGD loop and C-terminus, we also found that short and medium disintegrins shared a similar interaction in their C-termini (Figure 4D–F, Tables S2–S4). The N42 residue can interact with the mainchain amide proton of the P43 residue in the Ech chain A (Figure 4D) and chain B (Figure 4E). The corresponding residues of trimestatin also formed hydrogen bonds between the sidechain amide proton of residue W67 and the carboxyl group of residue N68 (Figure 4F). Thus, we know that these conserved regions of the C-terminus are essential for its structure. Furthermore, in contrast to trimestatin, the extended C-terminal tail residue H44 can interact with N42 in Ech chain A (Figure 4D) and chain B (Figure 4E). These results showed that the conformation of the C-terminal regions of Ech may be important in integrin recognition.

2.4. Structural Differences between the Integrin Complexes of Ech Chain A and Chain B

The docking of Ech chain A and chain B into integrin αvβ3 was used to identify the integrin-interacting residues. The analysis showed that the docking of Ech chain A and chain B into integrin αvβ3 resulted in different binding modes, especially in the C-terminus. (Figure 5A,B). The docking of Ech chain A into integrin αvβ3 showed that the side chain of the R22 residue of Ech chain A may interact with the main chain of N313 in the β3 subunit and that the main chain of the G25 residue of Ech chain A can also interact with the side chain of D218 in the αv subunit (Table 3 and Table S5). In contrast, Ech chain B did not show these kinds of interactions but instead interaction between the side chain of M28 and the main chain of D126 in the β3 subunit (Table 4 and Table S6). On the other hand, the main chain of C-terminal residue H44 of Ech chain A can interact with Y122 of the β3 subunit, and this contact was not found in Ech chain B. However, the side chain of C-terminal residue K45 of Ech chain B can interact with the side chain of M180 of the β3 subunit. These results showed that the C-terminal tails of Ech chains A and B were involved in Ech and β3 subunit interactions. These results suggested that Ech with a longer C-terminal tail may increase its recognition site and binding affinity to αvβ3 integrin.

2.5. The C-Terminal Residues of Ech are Involved in Mediating Integrin Recognition

The lengths of the C-terminus and the K45 residue of Ech were found to be important for interacting with integrin αIIbβ3 [13]. In this study, the inhibition of the binding of cell-expressed integrins αvβ3, αIIbβ3, αvβ5, and α5β1 to their ligands by Ech and its mutants was used to determine their activity and selectivity (Table 5 and Table 6). Compared to wild-type Ech, the mutant with the truncation of the C-terminal HKGPAT tail caused 6.4-, 7.0-, 11.7-, and 18.6-fold decreases in the inhibition of integrins αvβ3, αIIbβ3, αvβ5, and α5β1, respectively (Table 5). The effect of the C-terminal region on the change in the relative binding affinity to integrins was α5β1 > αvβ5 > αIIbβ3 > αvβ3. These results indicated that the C-terminal tail of Ech played an important role in recognizing not only integrin αIIbβ3 but also integrins αvβ3, αvβ5, and α5β1 [7,10,13].
Alanine scanning mutagenesis of the C-terminal tail of Ech showed that wild-type, P43A, H44A, and K45A mutants exhibited similar activity in inhibiting the adhesion of CHO cells that expressed integrin αvβ3 to immobilized fibrinogen with IC50 values of 20.7, 13.8, 16.4, and 14.7 nM (Table 6). In contrast to integrin αvβ3, some differences were found in their activities in inhibiting integrins αIIbβ3 and α5β1. They inhibited the adhesion of CHO cells that expressed integrin αIIbβ3 to immobilized fibrinogen with IC50 values of 51.5, 67.7, 20.7, and 132.5 nM and inhibited K562 cell adhesion to immobilized fibronectin with IC50 values of 132.6, 58.5. 30.0, and 84.9 nM. Interestingly, the H44A mutant caused 2.5- and 4.4-fold increases in the inhibition of αIIbβ3 and α5β1, and the K45A mutant caused a 2.6-fold decrease in the inhibition of integrin αIIbβ3. These results suggest that residues H44 and K45 may play different roles in integrin recognition. Indeed, in our docking structure of integrin αvβ3 with Ech chain A, we found that the main chain of C-terminal residue H44 in Ech chain A formed a hydrogen bond with the side chain of residue Y122 in the β3 subunit (Figure 5A, Table 3 and Table S5). Moreover, in the docking structure of integrin αvβ3 with Ech chain B, we found that the sidechain of C-terminal residue K45 in Ech chain B formed a hydrogen bond with the side chain of residue M180 in the β3 subunit (Figure 5B, Table 4 and Table S6). Our findings suggest that a longer C-terminal tail may improve the recognition site of Ech and increase its binding affinity to αvβ3 integrin, endowing Ech with inhibitory activity against αvβ3 integrin-related cancer. In the future, we can apply this property in Ech-based drug design for cancer therapy.

2.6. Inhibition of VEGF-Induced HUVEC Proliferation by Ech

Angiogenesis is a critical process for cancer cell growth and metastasis [35]. Endothelial cell proliferation is the key step in angiogenesis progression. Many studies have shown that disintegrins are potent angiogenic inhibitors [25,26,27,28,29,30,31,32]. The results of our analysis showed that Ech could inhibit HUVEC proliferation with an IC50 value of 103.2 nM (Figure 6, Table 7). In comparison with the reported disintegrin with IC50 values of 100–400 nM [17,30,31,32], Ech exhibited similar activity in inhibiting HUVEC proliferation (Table 7). These results suggest that Ech could be a good scaffold to design potent antiangiogenic agents.

2.7. C-Terminal Residue K45 of Ech Is Responsible for Its Anti-Metastatic Ability

Many studies have shown that disintegrins inhibit the adhesion and migration of tumor cells [24,25,26,29]. In our previous study, we also demonstrated that recombinant Ech effectively inhibited human A375 melanoma with an IC50 value of 1.5 nM (Table 8) [13]. To study the inhibitory effect of Ech on the migration of other tumor cells, we examined its antimigration activities in human glioblastoma U373MG and pancreatic Panc-1 cells. Similar to the results of melanoma cell inhibition by Ech, it inhibited the migration of U373MG and Panc-1 tumor cells with IC50 values of 5.7 and 154.5 nM, respectively (Figure 7 and Figure 8, Table 8). We found that the K45E mutant exhibited 8.9-, 4.8-, and 10.4-fold decreases in inhibiting the migration of A375, U373MG, and Panc-1 tumor cells (Figure 7 and Figure 8, Table 8), showing the importance of the K45 residue in inhibiting tumor cell migration. According to our X-ray structure, the conformation of Ech’s C-terminus was stabilized by the carboxyl group of residue M28 and the side-chain Nζ proton of residue K45 and between the side chain of residue D30 and the side-chain Nζ proton of residue K45 (Figure 4B). Furthermore, our docking structure of integrin αvβ3 with Ech chain B also showed that the side chain of C-terminal residue K45 in Ech chain B can interact with the side chain of residue M180 in the β3 subunit (Figure 5B, Table 4 and Table S6). These findings suggested that the conformation of the C-terminus plays an important role in αvβ3 integrin recognition, resulting in the inhibition of tumor cell migration.

3. Discussion

RGD-containing disintegrins have potent antiangiogenic and antimetastatic activities [23,24,25,26,27,28,29,30,31,32]. To elucidate the structure–function relationships of Ech, we determined the crystal structure of Ech and its anti-integrin, antiangiogenic, and antimigration activities. We found that the conformation of the RGD loop and the interactions between the RGD loop and C-terminus were similar to those in the reported medium disintegrin trimestatin [15]. In contrast, significant differences between short and medium disintegrins were found in their C-terminal conformations. We also showed that the C-terminus of Ech plays important roles in integrin recognition and the inhibition of tumor cell migration. This structural and functional study of Ech can serve as a basis for insight into its structure–function relationships and for the design of anticancer drugs.
Ech is the most studied disintegrin, and Aggrastat (Tirofilban), an antiplatelet drug, was designed based on its RGD loop structure [33]. It is known that the RGD loop, linker, and C-terminal regions of disintegrins are involved in integrin recognition [18,25,36]. The incorporation of C-terminal structural information into current drugs may improve its potency and selectivity [37]. Previous studies also showed that the C-terminal residues of Ech possess a functional role in modulating the binding affinity to integrin αIIbβ3 [6,7,8,9,10,13]. The analysis of Ech’s crystal structure showed that the sidechain of the D30 residue in the RGD loop and the amide proton of the N42 residue in the C-terminus formed hydrogen bonds. This is consistent with the report that the C-terminal tail of Ech acts in synergy with the RGD loop to interact with integrins [6,7,10,13]. According to our X-ray crystal structure and NMR NOE analysis, the residues from G46 to T49 of the C-terminus are flexible and cannot be observed in the X-ray structure. Disordered regions of proteins frequently have important biological functions. Interestingly, the carboxyl group of residue M28 interacts with the side-chain Nζ proton of residue K45, and the side chain of residue D30 interacts with the side chain Nζ proton of residue K45, resulting in a unique turn structure for Ech. These results showed that the K45 residue stabilizes the conformation of Ech’s C-terminus and is important for its integrin recognition. The unique conformation of Ech’s C-terminus can be used to design integrin-specific drugs.
Alanine scanning mutagenesis of the C-terminal tail of Ech showed that the H44A mutant caused 2.5- and 4.4-fold increases in the inhibition of αIIbβ3 and α5β1, and the K45A mutant caused a 2.6-fold decrease in the inhibition of integrin αIIbβ3. These results suggested that residues H44 and K45 may play different roles in integrin recognition. Moreover, we found a correlation between our docking results and these functional studies. In Ech chain A, the main chain of C-terminal residue H44 formed a hydrogen bond with the side chain of residue Y122 in the β3 subunit (Figure 5A, Table 3 and Table S5). On the other hand, in Ech chain B, the sidechain of C-terminal residue K45 interacts with the sidechain of residue M180 in the β3 subunit (Figure 5B, Table 4 and Table S6). The different interaction modes between H44/K45 and αvβ3 integrin may provide information on how H44/K45 binds to other integrins.
Integrins are highly expressed in many cancers, including melanoma, pancreatic cancer, glioblastoma, breast cancer, prostate cancer, ovarian cancer, and cervical cancer, and are involved in cancer cell progression [20,21,22,23]. A number of studies have reported that Ech is a potent antagonist of integrins αvβ3, αIIbβ3, and α5β1 for cancer therapy and inhibits HUVEC adhesion to immobilized fibronectin and vitronectin, as well as cell proliferation, migration, invasion, and adhesion of metastatic human osteosarcoma [6,10,24,25,26,38]. T cells genetically engrafted with Ech-containing chimeric antigen receptor could also efficiently inhibit integrin-expressing solid tumors [39]. We showed that Ech could inhibit VEGF-induced HUVEC proliferation with an IC50 value of 103.2 nM and inhibit the migration of human melanoma A375, glioblastoma U373MG, and pancreatic Panc-1 tumor cells with IC50 values of 1.5, 5.7, and 154.5 nM. These findings indicate that Ech is a potent antiangiogenic and antimetastatic agent. The resulting structure–function relationships found in this study can provide a molecular basis for the use of Ech as a molecular scaffold to design anticancer drugs in the future.

4. Conclusions

In conclusion, this is the first report to determine a high-resolution crystal structure of a short disintegrin with a resolution of 1.80 Å. In comparison with the crystal structure of medium disintegrin, they have similar conformations for the RGD loop and the same hydrogen bond connectivity between the RGD loop and C-terminus. The major difference was found from the conformation of their C-termini. According to the results of the functional analysis, Ech may be a potential candidate for anticancer therapeutics with antiangiogenic and antimigration properties, and the C-terminal region of Ech can be optimized to improve its anticancer activity. The structure–function relationships of Ech can provide a molecular basis to design anticancer drugs using Ech as a molecular scaffold.

5. Materials and Methods

5.1. Expression, Purification, and Characterization of Ech and Its Mutants

The expression and purification of Ech and its mutants in P. pastoris were performed following previously described protocols [13,14,18,40,41,42]. In brief, Ech and its mutants were expressed in the P. pastoris X33 strain with a 6xHis-tag, and purified to apparent homogeneity by Ni2+-chelating chromatography and HPLC. Wild-type Ech without a 6xHis-tag was expressed for crystallization and purified to homogeneity by CaptoMMC chromatography and HPLC. The yields of Ech and its mutants produced in P. pastoris were 2–7 mg/L. Mass spectrometry analysis showed that the experimental molecular weights of all recombinant proteins were in excellent agreement with theoretical values (less than 1 Da), which were calculated by assuming that all cysteines form disulfide bonds. Thus, our results indicated the formation of four disulfide bonds in Ech and its mutants.

5.2. Mass Spectrometric Measurements

The molecular weights of Ech and its mutants were confirmed using an LTQ Orbitrap XL mass spectrometer equipped with an electrospray ionization source (Thermo Fisher Scientific, Waltham, MA, USA) as previously described [18].

5.3. Cell Adhesion Assay

The cell adhesion assay was performed using previously described protocols [14]. In brief, Chinese hamster ovary (CHO) cells expressing integrins αvβ3 (CHO-αvβ3) and αIIbβ3 (CHO-αIIbβ3) were kindly provided by Dr Y. Takada (Scripps Research Institute) and maintained in Dulbecco’s Modified Eagle’s medium (DMEM). Human erythroleukemia K562 cells and the human colon adenocarcinoma cell line HT-29 were purchased from ATCC. K562 cells were cultured in RPMI-1640 medium containing 5% fetal calf serum, and HT29 cells were maintained in DMEM. Ninety-six-well microtiter plates (Costar, Corning, NY, USA) were coated with 100 μL of PBS buffer containing 200 μg/mL fibrinogen (for integrins αvβ3 and αIIbβ3), 50 μg/mL fibronectin (for integrins α5β1), and 5 μg/mL vitronectin (for integrins αvβ5) at 4 °C overnight. Wells were blocked by incubating with 200 μL of heat-denatured 1% bovine serum albumin (BSA) (Calbiochem, Darmstadt, Germany) at room temperature for 40 min. The heat-denatured BSA was discarded and was washed twice with 200 μL of PBS. CHO-αvβ3, CHO-αIIbβ3, and HT29 cells were detached by trypsinization and diluted to 3 × 105 cells/mL; K562 cells were diluted to 2.5 × 105 cells/mL; and 105 μL of cells was used for the assay. Cells were premixed with or without Ech and its mutants for 15–30 min at 37 °C in a 5% CO2 atmosphere prior to dispensation into microtiter wells. The cells were then seeded to the coated plate and incubated at 37 °C (5% CO2) for 1 h. The non-adhered cells were removed by washing twice with 200 μL of PBS. The adhered cells were fixed with 100 μL of 3.7% paraformaldehyde and quantified using a crystal violet assay. Colorization was performed by adding 150 μL of colorizing solution containing 50% alcohol and 0.1% acetic acid. The absorbance was read at 600 nm by an automatic microplate reader. Inhibition was defined as inhibition = 100 – [OD600 (Rho protein-treated sample)/OD600 (untreated sample)] × 100%. The reported IC50 values are the average of at least three separate experiments.

5.4. Crystallization of Ech

Ech was dialyzed twice against 1 L ddH2O at 4 °C, lyophilized, and then dissolved in 50 mM Tris–HCl buffer (pH 8.0) at a concentration of ~10 mg/mL. Ech was initially screened using commercial kits (The JCSG Core Suites, QIAGEN, Hilden, Germany) under 384 conditions. The CombiClover Junior sitting plate (24-well, Emerald BioSystems, Bainbridge Island, WA, USA) was used for crystallization. The crystals of Ech were obtained by mixing 0.8 μL protein solution and the reservoir solution containing 2.4 M (NH4)2SO4 (pH 9.0) and 0.1 M Bicine at pH 9.0 by the sitting-drop vapor diffusion method.

5.5. Diffraction Data Collection and Processing

All crystals were flash cooled to 100 K in a stream of cold nitrogen prior to data collection. X-ray diffraction data were collected on beamline TPS 05A at the National Synchrotron Radiation Research Center (NSRRC, Taiwan). The diffraction images were processed using the XDS program package [43]. The crystal data and data collection statistics of Ech are summarized in Table 1.

5.6. Structure Determination and Refinement

The initial phase of Ech was determined by molecular replacement, and the known structure of Rhodostomin (PDB ID 4RQG) in the PDB was used as our model protein. All structure determinations were carried out using alternate cycles of model building in Coot [44] and refinement in REFMAC5 [45]. The stereochemistry and structures of the final models were analyzed by MolProbity [34]. Ech belonged to the P3221 space group with parameters a = 33.79, b = 33.79, and c = 120.90 Å. This structure was subsequently refined at 1.80 Å resolution using REFMAC5 to an R-factor of 0.212 and R-free of 0.255. The crystal structure of Ech was deposited in PDB ID 6LSQ.

5.7. NOE Analysis

The NOE experiment was conducted according to previously described protocols [14]. In brief, the backbone dynamics of Ech were determined by two-dimensional proton-detected heteronuclear NMR spectroscopy. NMR experiments were performed at 27 °C on a Bruker Avance 600 MHz spectrometer. The steady-state 1H-15N NOEs were collected from 1H-detected 1H-15N correlation spectra recorded with sensitivity-enhanced pulse sequences. In the NOE experiment, two spectra, one with the saturated NOE and one without, were collected. The steady-state NOEs were the ratios that were calculated between the crosspeak intensities with and without saturation of the proton resonances (i.e., NOE = Isat/Inonsat). The reported NOE value was the average value of three experiments.

5.8. Molecular Docking

The HADDOCK webserver was used for the docking of Ech into integrin αvβ3 [46]. The starting structure for docking was the X-ray structure of Echistatin. The interaction restraints were derived from the X-ray structure of integrin αvβ3 in complex with the 10th type III RGD domain of wild-type fibronectin (PDB ID 4MMX) [47] by using CCP4i software [48]. The defined distance threshold was 4 Å, and the interaction restraints between the RGD motif and integrin were used for the calculation. The input restraints between the R49, G50, and D51 residues of Ech and integrin αvβ3 were 27, 11, and 40, respectively. They were the contacts between the R49 residue and residues D150, Y178, Q180, and D218 of integrin αv; between residues G50 and Y178 of αv and residues R216 and A218 of β3; and between residue D51 and residues S121, Y122, S123, R214, N215, R216, D217, and A218 and Mn2+ of the MIDAS of β3. An additional 0.5 Å distance was defined as the minimum and maximum limits of the interaction restraints. These restraints were used to perform the standard HADDOCK protocol for protein docking with minor modifications. This protocol combines three stages of molecular dynamics calculations, including heating and cooling, with a progressive increase in flexibility at the binding interface. In the first stage, 1000 conformations were calculated using a rigid-body docking protocol. The 200 structures with the lowest intermolecular energies were refined using semiflexible simulated annealing in the second stage. Semiflexible and fully flexible regions were automatically defined by HADDOCK. The 200 structures with the lowest intermolecular energy values were refined using explicit water molecules in the last stage. The clusters of structures were identified by fitting them with an average root mean square deviation (RMSD) value of 0.5–1.5 Å for the backbone atoms of all the amino acids. The selected structure cluster for the structural analysis was based on the lowest Z-score without any restraint violations. Hydrogen bonds and salt bridges interactions were analyzed with PISA software [49].

5.9. VEGF-Induced HUVEC Proliferation Assay

Human umbilical vein endothelial cells (HUVECs) were purchased from Cambrex-Lonza (East Rutherford, NJ, USA). They were maintained and expanded in endothelial cell growth medium-2 (EGM-2) and gradually replaced by 20% CCS/M199. For the proliferation assay, cells were seeded in 0.2% gelatin-coated 96-well plates at 1.5 × 104 cells/well, and then the cells were cultured in M199 serum-free medium at 37 °C overnight. The following day, 20 ng/mL recombinant human VEGF (R&D Systems, Minneapolis, MN, USA) was added, and the cells were treated with different concentrations of Ech for 48 h. After removing the medium, the cells of the plate were fixed with 120 μL of 3.7% paraformaldehyde and stained with 50 μL of 0.05% crystal violet solution for 20 min. Each well was then washed four times with 200 μL of distilled water and dried. Colorization was performed by adding 150 μL of colorizing solution containing 50% alcohol and 0.1% acetic acid. The absorbance was read at 600 nm by an automatic microplate reader. The reported values are the average of at least three separate experiments.

5.10. Cell Migration Assay

The transwell migration assay was conducted according to previously described protocols [13]. In brief, human A375 melanoma cells, human pancreatic Panc-1 cells, and U373 MG human glioblastoma cells were used to determine the inhibitory activities of recombinant Ech and its mutant. Then, 1 × 105 cells were added to the transwell inserts and incubated in 24-well plates. Cells in the inserts were cultured in serum free DMEM medium, whereas the lower wells were filled with DMEM medium containing 10% FBS as an attractant. After 6 h of incubation, the non-migrated cells on the upper side of the inserts were removed using a cotton swab. The insert wells were fixed with 3.7% paraformaldehyde and stained with 0.05% crystal violet. Cells that had transmigrated to the lower surface were imaged at 100× magnification, for five random fields by using OLYMPUS-IX71 microscopy. Migrated cells were counted using NIH ImageJ software.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/12/11/709/s1, Figure S1: Crystal packing of Ech, Figure S2: 1H-15N heteronuclear NOE values measured for individual amide sites along the backbones, Figure S3: The interactions of Ech with water molecules and sulfate, Table S1: Summary of the interactions between Ech (chain A and chain B) and sulfate ion, Table S2: Summary of the interactions in the RGD loop, C-terminus, and between the RGD loop and C-terminus in Ech chain A, Table S3: Summary of the interactions in the RGD loop, C-terminus, and between the RGD loop and C-terminus in Ech chain B, Table S4: Summary of the interactions in the RGD loop, C-terminus, and between the RGD loop and C-terminus in trimestatin, Table S5: Summary of the interactions between Ech chain A and αvβ3 integrin, Table S6: Summary of the interactions between Ech chain B and αvβ3 integrin.

Author Contributions

Conceptualization, Y.-C.C. and W.-J.C.; Validation, Y.-C.C., Y.-T.C., J.-F.C. and W.-J.C.; Data analysis, Y.-C.C., J.-H.S., C.-H.H. and W.-J.C.; Resources, Y.-C.C., Y.-T.C., J.-F.C., C.-Y.C., C.-H.C. and W.-J.C.; Data curation, Y.-C.C., Y.-T.C. J.-F.C. and W.-J.C.; Supervision, W.-J.C.; Writing—original draft, Y.-C.C. Writing—review & editing, Y.-C.C. and W.-J.C.; Funding acquisition, W.-J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the research grants from the Ministry of Science and Technology (MOST-108-2622-B006-005-CC2), Taiwan, Republic of China and by the Program for Translational Innovation of Biopharmaceutical Development - Technology Supporting Platform Axis (AS-KPQ-106-TSPA).

Acknowledgments

The authors would like to say special thanks to Wen-Yih Jeng and Yu-Fang Su for their professional assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gan, Z.R.; Gould, R.J.; Jacobs, J.W.; Friedman, P.A.; Polokoff, M.A. A potent platelet aggregation inhibitor from the venom of the viper, Echis carinatus. J. Biol. Chem. 1988, 263, 19827–19832. [Google Scholar] [PubMed]
  2. Savage, B.; Marzec, U.M.; Chao, B.H.; Harker, L.A.; Maraganore, J.M.; Ruggeri, Z.M. Binding of the snake venom-derived proteins applaggin and echistatin to the arginine-glycine-aspartic acid recognition site(s) on platelet glycoprotein IIb.IIIa complex inhibits receptor function. J. Biol. Chem. 1990, 265, 11766–11772. [Google Scholar]
  3. Atkinson, R.; Saudek, V.; Pelton, J. Echistatin: The refined structure of a disintegrin in solution by 1 H NMR and restrained molecular dynamics. Int. J. Pept. Protein Res. 2009, 43, 563–572. [Google Scholar] [CrossRef]
  4. Garsky, V.M.; Lumma, P.K.; Freidinger, R.M.; Pitzenberger, S.M.; Randall, W.C.; Veber, D.F.; Gould, R.J.; Friedman, P.A. Chemical synthesis of echistatin, a potent inhibitor of platelet aggregation from Echis carinatus: Synthesis and biological activity of selected analogs. Proc. Natl. Acad. Sci. USA 1989, 86, 4022–4026. [Google Scholar] [CrossRef] [Green Version]
  5. Cooke, R.M.; Carter, B.G.; Martin, D.M.A.; Murray-Rust, P.; Weir, M.P. Nuclear magnetic resonance studies of the snake toxin echistatin. 1H resonance assignments and secondary structure. JBIC J. Biol. Inorg. Chem. 1991, 202, 323–328. [Google Scholar] [CrossRef]
  6. Monleon, D.; Esteve, V.; Kovacs, H.; Calvete, J.J.; Celda, B. Conformation and concerted dynamics of the integrin-binding site and the C-terminal region of echistatin revealed by homonuclear NMR. Biochem. J. 2005, 387, 57–66. [Google Scholar] [CrossRef]
  7. Marcinkiewicz, C.; Vijay-Kumar, S.; McLane, M.A.; Niewiarowski, S. Significance of RGD loop and C-terminal domain of echistatin for recognition of alphaIIb beta3 and alpha(v) beta3 integrins and expression of ligand-induced binding site. Blood 1997, 90, 1565–1575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Yamada, T.; Kidera, A. Tailoring echistatin to possess higher affinity for integrin alpha (IIb)beta (3). FEBS Lett. 1996, 387, 11–15. [Google Scholar] [CrossRef] [Green Version]
  9. McLane, M.A.; Kowalska, M.A.; Silver, L.; Shattil, S.J.; Niewiarowski, S. Interaction of disintegrins with the αIIbβ3 receptor on resting and activated human platelets. Biochem. J. 1994, 301, 429–436. [Google Scholar] [CrossRef]
  10. Wierzbicka-Patynowski, I.; Niewiarowski, S.; Marcinkiewicz, C.; Calvete, J.J.; Marcinkiewicz, M.M.; McLane, M.A. Structural requirements of echistatin for the recognition of alpha(v)beta(3) and alpha(5)beta(1) integrins. J. Biol. Chem. 1999, 274, 37809–37814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Zannetti, A.; Del Vecchio, S.; Iommelli, F.; Del Gatto, A.; De Luca, S.; Zaccaro, L.; Papaccioli, A.; Sommella, J.; Panico, M.; Speranza, A.; et al. Imaging of v 3 Expression by a Bifunctional Chimeric RGD Peptide not Cross-Reacting with v 5. Clin. Cancer Res. 2009, 15, 5224–5233. [Google Scholar] [CrossRef] [Green Version]
  12. Minea, R.O.; Helchowski, C.M.; Zidovetzki, S.J.; Costa, F.K.; Swenson, S.D.; Markland, F. Vicrostatin – An Anti-Invasive Multi-Integrin Targeting Chimeric Disintegrin with Tumor Anti-Angiogenic and Pro-Apoptotic Activities. PLoS ONE 2010, 5, e10929. [Google Scholar] [CrossRef]
  13. Chen, Y.-C.; Cheng, C.-H.; Shiu, J.-H.; Chang, Y.-T.; Chang, Y.-S.; Chuang, W.-J.; Lee, J.-C. Expression in Pichia pastoris and characterization of echistatin, an RGD-containing short disintegrin. Toxicon 2012, 60, 1342–1348. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, C.-Y.; Shiu, J.-H.; Hsieh, Y.-H.; Liu, Y.-C.; Chen, Y.-C.; Chen, Y.-C.; Jeng, W.-Y.; Tang, M.-J.; Lo, S.J.; Chuang, W.-J. Effect of D to E mutation of the RGD motif in rhodostomin on its activity, structure, and dynamics: Importance of the interactions between the D residue and integrin. Proteins Struct. Funct. Bioinform. 2009, 76, 808–821. [Google Scholar] [CrossRef]
  15. Fujii, Y.; Okuda, D.; Fujimoto, Z.; Horii, K.; Morita, T.; Mizuno, H. Crystal Structure of Trimestatin, a Disintegrin Containing a Cell Adhesion Recognition Motif RGD. J. Mol. Biol. 2003, 332, 1115–1122. [Google Scholar] [CrossRef]
  16. Senn, H.; Klaus, W. The Nuclear Magnetic Resonance Solution Structure of Flavoridin, an Antagonist of the Platelet GP IIb-IIIa Receptor. J. Mol. Biol. 1993, 232, 907–925. [Google Scholar] [CrossRef]
  17. Shin, J.; Hong, S.Y.; Chung, K.; Kang, I.; Jang, Y.; Kim, D.S.; Lee, W. Solution Structure of a Novel Disintegrin, Salmosin, fromAgkistrondon halysVenom. Biochemistry 2003, 42, 14408–14415. [Google Scholar] [CrossRef]
  18. Chang, Y.-T.; Shiu, J.-H.; Chuang, W.-J.; Chen, Y.-C.; Chen, C.-Y.; Chang, Y.-S. Effects of the RGD loop and C-terminus of rhodostomin on regulating integrin αIIbβ3 recognition. PLoS ONE 2017, 12, e0175321. [Google Scholar] [CrossRef] [Green Version]
  19. Saudek, V.; Atkinson, R.A.; Pelton, J.T. Three-dimensional structure of echistatin, the smallest active RGD protein. Biochemistry 1991, 30, 7369–7372. [Google Scholar] [CrossRef]
  20. Rádis-Baptista, G. Integrins, cancer and snake toxins (mini-review). J. Venom. Anim. Toxins Incl. Trop. Dis. 2005, 11, 217–241. [Google Scholar] [CrossRef]
  21. Alday-Parejo, B.; Stupp, R.; Rüegg, C. Are Integrins Still Practicable Targets for Anti-Cancer Therapy? Cancers 2019, 11, 978. [Google Scholar] [CrossRef] [Green Version]
  22. Desgrosellier, J.S.; Cheresh, D.A. Integrins in cancer: Biological implications and therapeutic opportunities. Nat. Rev. Cancer. 2010, 10, 9–22. [Google Scholar] [CrossRef] [Green Version]
  23. Nieberler, M.; Reuning, U.; Reichart, F.; Notni, J.; Wester, H.-J.; Schwaiger, M.; Weinmüller, M.; Räder, A.; Steiger, K.; Kessler, H. Exploring the Role of RGD-Recognizing Integrins in Cancer. Cancers 2017, 9, 116. [Google Scholar] [CrossRef]
  24. Selistre-De-Araújo, H.S.; Pontes, C.L.S.; Montenegro, C.F.; Martin, A.C.B.M. Snake Venom Disintegrins and Cell Migration. Toxins 2010, 2, 2606–2621. [Google Scholar] [CrossRef]
  25. McLane, M.A. Disintegrins in health and disease. Front. Biosci. 2008, 13, 6617–6637. [Google Scholar] [CrossRef]
  26. McLane, M.A.; Sanchez, E.E.; Wong, A.; Paquette-Straub, C.; Perez, J.C. Disintegrins. Curr. Drug Targets Cardiovasc. Haematol. Disord. 2004, 4, 327–355. [Google Scholar] [CrossRef] [PubMed]
  27. Rivas-Mercado, E.; Garza-Ocañas, L. Disintegrins obtained from snake venom and their pharmacological potential. Med. Univ. 2017, 19, 32–37. [Google Scholar] [CrossRef]
  28. Macêdo, J.K.; Fox, J.W.; Castro, M.D.S. Disintegrins from Snake Venoms and their Applications in Cancer Research and Therapy. Curr. Protein Pept. Sci. 2015, 16, 532–548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. David, V.; Succar, B.B.; De Moraes, J.A.; Saldanha-Gama, R.F.G.; Barja-Fidalgo, C.; Zingali, R.B. Recombinant and Chimeric Disintegrins in Preclinical Research. Toxins 2018, 10, 321. [Google Scholar] [CrossRef] [Green Version]
  30. Hong, S.-Y.; Koh, Y.-S.; Chung, K.-H.; Kim, D.-S. Snake venom disintegrin, saxatilin, inhibits platelet aggregation, human umbilical vein endothelial cell proliferation, and smooth muscle cell migration. Thromb. Res. 2002, 105, 79–86. [Google Scholar] [CrossRef]
  31. Yeh, C.H.; Peng, H.C.; Yang, R.S.; Huang, T.F. Rhodostomin, a snake venom disintegrin, inhibits angiogenesis elicited by basic fibroblast growth factor and suppresses tumor growth by a selective alpha(v)beta (3) blockade of endothelial cells. Mol. Pharmacol. 2001, 59, 1333–1342. [Google Scholar] [CrossRef] [Green Version]
  32. Sheu, J.-R.; Yen, M.-H.; Kan, Y.-C.; Hung, W.-C.; Chang, P.-T.; Luk, H.-N. Inhibition of angiogenesis in vitro and in vivo: Comparison of the relative activities of triflavin, an Arg-Gly-Asp-containing peptide and anti-αvβ3 integrin monoclonal antibody. Biochim. Biophys. Acta (BBA)—Gen. Subj. 1997, 1336, 445–454. [Google Scholar] [CrossRef]
  33. Lazarovici, P.; Marcinkiewicz, C.; Lelkes, P.I. From Snake Venom’s Disintegrins and C-Type Lectins to Anti-Platelet Drugs. Toxins 2019, 11, 303. [Google Scholar] [CrossRef] [Green Version]
  34. Williams, C.J.; Headd, J.J.; Moriarty, N.W.; Prisant, M.G.; Videau, L.L.; Deis, L.N.; Verma, V.; Keedy, D.A.; Hintze, B.J.; Chen, V.B.; et al. MolProbity: More and better reference data for improved all-atom structure validation. Protein Sci. 2018, 27, 293–315. [Google Scholar] [CrossRef]
  35. Nishida, N.; Yano, H.; Nishida, T.; Kamura, T.; Kojiro, M. Angiogenesis in cancer. Vasc. Health Risk Manag. 2006, 2, 213–219. [Google Scholar] [CrossRef]
  36. Sumathipala, R.; Xu, C.; Seago, J.; Mould, A.P.; Humphries, M.J.; Craig, S.E.; Patel, Y.; Wijelath, E.S.; Sobel, M.; Rahman, S. The "Linker" Region (Amino Acids 38-47) of the Disintegrin Elegantin Is a Novel Inhibitory Domain of Integrin 5beta1-Dependent Cell Adhesion on Fibronectin: EVIDENCE FOR THE NEGATIVE REGULATION OF FIBRONECTIN SYNERGY SITE BIOLOGICAL ACTIVITY. J. Biol. Chem. 2006, 281, 37686–37696. [Google Scholar] [CrossRef] [Green Version]
  37. Swenson, S.D.; Minea, R.O.; Tuan, C.D.; Thein, T.-Z.; Chen, T.C.; Markland, F. A Novel Venom-Derived Peptide for Brachytherapy of Glioblastoma: Preclinical Studies in Mice. Molecules 2018, 23, 2918. [Google Scholar] [CrossRef] [Green Version]
  38. Tome, Y.; Kimura, H.; Kiyuna, T.; Sugimoto, N.; Tsuchiya, H.; Kanaya, F.; Bouvet, M.; Hoffman, R.M. Disintegrin targeting of an αvβ3 integrin-over-expressing high-metastatic human osteosarcoma with echistatin inhibits cell proliferation, migration, invasion and adhesion in vitro. Oncotarget 2016, 7, 46315–46320. [Google Scholar] [CrossRef] [Green Version]
  39. Fu, X.; Rivera, A.; Tao, L.; Zhang, X. Genetically modified T cells targeting neovasculature efficiently destroy tumor blood vessels, shrink established solid tumors and increase nanoparticle delivery. Int. J. Cancer 2013, 133, 2483–2492. [Google Scholar] [CrossRef]
  40. Guo, R.-T.; Chou, L.-J.; Chen, Y.-C.; Chen, C.-Y.; Pari, K.; Jen, C.J.; Lo, S.J.; Huang, S.-L.; Lee, C.-Y.; Chang, T.-W.; et al. Expression in Pichia pastoris and characterization by circular dichroism and NMR of rhodostomin. Proteins Struct. Funct. Bioinform. 2001, 43, 499–508. [Google Scholar] [CrossRef]
  41. Chen, C.-Y.; Cheng, C.-H.; Chen, Y.-C.; Lee, J.-C.; Choua, S.-H.; Huang, W.; Chuang, W.-J. Preparation of amino-acid-type selective isotope labeling of protein expressed in Pichia pastoris. Proteins Struct. Funct. Bioinform. 2005, 62, 279–287. [Google Scholar] [CrossRef]
  42. Shiu, J.-H.; Chen, C.-Y.; Chen, Y.-C.; Chang, Y.-T.; Chang, Y.-S.; Huang, C.-H.; Chuang, W.-J. Effect of P to A Mutation of the N-Terminal Residue Adjacent to the Rgd Motif on Rhodostomin: Importance of Dynamics in Integrin Recognition. PLoS ONE 2012, 7, e28833. [Google Scholar] [CrossRef] [Green Version]
  43. Kabsch, W. XDS. Acta Cryst. D 2010, 66, 125–132. [Google Scholar] [CrossRef] [Green Version]
  44. Emsley, P.; Cowtan, K. Coot: Model-building tools for molecular graphics. Acta Crystallogr. Sect. D Biol. Crystallogr. 2004, 60, 2126–2132. [Google Scholar] [CrossRef] [Green Version]
  45. Murshudov, G.N.; Skubák, P.; Lebedev, A.A.; Pannu, N.S.; Steiner, R.A.; Nicholls, R.A.; Winn, M.D.; Long, F.; Vagin, A.A. REFMAC5 for the refinement of macromolecular crystal structures. Acta Crystallogr. Sect. D Biol. Crystallogr. 2011, 67, 355–367. [Google Scholar] [CrossRef] [Green Version]
  46. De Vries, S.J.; Van Dijk, M.; Bonvin, A.M.J.J. The HADDOCK web server for data-driven biomolecular docking. Nat. Protoc. 2010, 5, 883–897. [Google Scholar] [CrossRef] [Green Version]
  47. Van Agthoven, J.F.; Xiong, J.-P.; Alonso, J.L.; Rui, X.; Adair, B.D.; Goodman, S.L.; Arnaout, M.A. Structural basis for pure antagonism of integrin αVβ3 by a high-affinity form of fibronectin. Nat. Struct. Mol. Biol. 2014, 21, 383–388. [Google Scholar] [CrossRef] [Green Version]
  48. CCP4i Software. Available online: http://structure.usc.edu/ccp4/ (accessed on 9 May 2020).
  49. PISA Server. Available online: https://www.ebi.ac.uk/pdbe/pisa/pi_link.html (accessed on 22 March 2020).
Toxins 12 00709 g001 550
Figure 1. Electron density map of the RGD loop and C-termini of Ech. The 2|Fo|- |Fc| electron density map at 1.80 Å resolution is contoured at 1.0 σ around the residues from the RGD loop (K21-D30) and C-terminus (P40-T49). (A) Electron density map of the RGD loop with residues K21 to D30 (chain A); (B) Electron density map of the RGD loop with residues K21 to D30 (chain B); (C) Electron density map of the C-terminus with residues P40 to K45 (chain A); (D) Electron density map of the C-terminus with residues P40 to T49 (chain B).
Figure 1. Electron density map of the RGD loop and C-termini of Ech. The 2|Fo|- |Fc| electron density map at 1.80 Å resolution is contoured at 1.0 σ around the residues from the RGD loop (K21-D30) and C-terminus (P40-T49). (A) Electron density map of the RGD loop with residues K21 to D30 (chain A); (B) Electron density map of the RGD loop with residues K21 to D30 (chain B); (C) Electron density map of the C-terminus with residues P40 to K45 (chain A); (D) Electron density map of the C-terminus with residues P40 to T49 (chain B).
Toxins 12 00709 g001
Toxins 12 00709 g002 550
Figure 2. Comparison of X-ray and NMR structures of Ech. The β-sheet secondary structure is shown in cyan. The side chains of the R and D residues in the RGD loop are shown in blue and red, and disulfide bond linkages are shown in red. (A) Crystal structure of Ech chain A; (B) Crystal structure of Ech chain B; (C) Superimposition of Ech chain A and B, shown in cyan and magenta, respectively; (D) Superimposition of Ech chain A and B RGD loop with residues K21 to D30; (E) Superimposition of Ech chain A and B C-terminus with residues P40 to T49; (F) NMR structure (2Ech) of Ech; (G) Superimposition of Ech chain A and NMR structure (2Ech). Chain A and the NMR structure are shown in cyan and gray, respectively.
Figure 2. Comparison of X-ray and NMR structures of Ech. The β-sheet secondary structure is shown in cyan. The side chains of the R and D residues in the RGD loop are shown in blue and red, and disulfide bond linkages are shown in red. (A) Crystal structure of Ech chain A; (B) Crystal structure of Ech chain B; (C) Superimposition of Ech chain A and B, shown in cyan and magenta, respectively; (D) Superimposition of Ech chain A and B RGD loop with residues K21 to D30; (E) Superimposition of Ech chain A and B C-terminus with residues P40 to T49; (F) NMR structure (2Ech) of Ech; (G) Superimposition of Ech chain A and NMR structure (2Ech). Chain A and the NMR structure are shown in cyan and gray, respectively.
Toxins 12 00709 g002
Toxins 12 00709 g003 550
Figure 3. Structural comparison of the RGD loops (K21-D30) and sequence alignment of Ech and trimestatin (1J2L). (A) The conformation of Ech’s RGD loop (chain A). The Cα to Cα distance between residues R24 and D26 is shown; (B) The conformation of Ech’s RGD loop (chain B). The Cα to Cα distance between residues R24 and D26 is shown; (C) The conformation of trimestatin’s RGD loop. The Cα to Cα distance between residues R49 and D51 is shown; (D) Hydrogen bonds (residues K21-D30, A23-D29, and R24-D27) of Ech’s RGD loop (chain A); (E) Hydrogen bonds (residues K21-D30, A23-D29, and R24-D27) of Ech’s RGD loop (chain B); (F) Hydrogen bonds (residues R46-D55 and A48-D54) of trimestatin’s RGD loop; (G) Sequence alignment of Ech and trimestatin. The conserved residues are boxed in gray.
Figure 3. Structural comparison of the RGD loops (K21-D30) and sequence alignment of Ech and trimestatin (1J2L). (A) The conformation of Ech’s RGD loop (chain A). The Cα to Cα distance between residues R24 and D26 is shown; (B) The conformation of Ech’s RGD loop (chain B). The Cα to Cα distance between residues R24 and D26 is shown; (C) The conformation of trimestatin’s RGD loop. The Cα to Cα distance between residues R49 and D51 is shown; (D) Hydrogen bonds (residues K21-D30, A23-D29, and R24-D27) of Ech’s RGD loop (chain A); (E) Hydrogen bonds (residues K21-D30, A23-D29, and R24-D27) of Ech’s RGD loop (chain B); (F) Hydrogen bonds (residues R46-D55 and A48-D54) of trimestatin’s RGD loop; (G) Sequence alignment of Ech and trimestatin. The conserved residues are boxed in gray.
Toxins 12 00709 g003
Toxins 12 00709 g004 550
Figure 4. Comparison of hydrogen bonding interactions between the RGD loop and C-terminus in Ech and trimestatin, as well as those within the C-terminus. (A) Hydrogen bond between RGD loop and C-terminus (chain A); The hydrogen bond between the side chain of the D30 residue in the RGD loop and the amide proton of the N42 residue in the C-terminus is shown; (B) Hydrogen bonds between the RGD loop and C-terminus (chain B). Three hydrogen bonds between the carboxyl group of residue M28 and the sidechain Nζ proton of residue K45, between the side chain of residue D30 and the amide proton of residue N42, and between the side chain of residue D30 and the sidechain Nζ proton of residue K45; (C) Hydrogen bond between the RGD loop and C-terminus in trimestatin. The hydrogen bond between the side chain of the D55 residue in the RGD loop and the amide proton of the W67 residue in the C-terminus is shown; (D) Interaction in Ech C-terminus (chain A). The interaction between the side chain of the N42 residue and the mainchain amide proton of residues P43 and H44; (E) Interaction in Ech C-terminus (chain B). The interaction between the side chain of the N42 residue and the mainchain amide proton of residues P43 and H44; (F) Hydrogen bond in the C-terminus of trimestatin. The hydrogen bond between the side chain of residue W67 and the backbone carboxyl group of residue N68.
Figure 4. Comparison of hydrogen bonding interactions between the RGD loop and C-terminus in Ech and trimestatin, as well as those within the C-terminus. (A) Hydrogen bond between RGD loop and C-terminus (chain A); The hydrogen bond between the side chain of the D30 residue in the RGD loop and the amide proton of the N42 residue in the C-terminus is shown; (B) Hydrogen bonds between the RGD loop and C-terminus (chain B). Three hydrogen bonds between the carboxyl group of residue M28 and the sidechain Nζ proton of residue K45, between the side chain of residue D30 and the amide proton of residue N42, and between the side chain of residue D30 and the sidechain Nζ proton of residue K45; (C) Hydrogen bond between the RGD loop and C-terminus in trimestatin. The hydrogen bond between the side chain of the D55 residue in the RGD loop and the amide proton of the W67 residue in the C-terminus is shown; (D) Interaction in Ech C-terminus (chain A). The interaction between the side chain of the N42 residue and the mainchain amide proton of residues P43 and H44; (E) Interaction in Ech C-terminus (chain B). The interaction between the side chain of the N42 residue and the mainchain amide proton of residues P43 and H44; (F) Hydrogen bond in the C-terminus of trimestatin. The hydrogen bond between the side chain of residue W67 and the backbone carboxyl group of residue N68.
Toxins 12 00709 g004
Toxins 12 00709 g005 550
Figure 5. The interactions between Ech and αvβ3 integrin. (A) Ech chain A (cyan) docked on αvβ3 (PDB code 1L5G, gray); the main chains and side chains of the RGD loop and C-terminal interacting residues are labeled; (B) Ech chain B (magenta) docked on αvβ3; the main chains and side chains of the RGD loop and C-terminal interacting residues are labeled.
Figure 5. The interactions between Ech and αvβ3 integrin. (A) Ech chain A (cyan) docked on αvβ3 (PDB code 1L5G, gray); the main chains and side chains of the RGD loop and C-terminal interacting residues are labeled; (B) Ech chain B (magenta) docked on αvβ3; the main chains and side chains of the RGD loop and C-terminal interacting residues are labeled.
Toxins 12 00709 g005
Toxins 12 00709 g006 550
Figure 6. Inhibition of VEGF-induced HUVEC proliferation by Ech. (A) Quantitative analysis of HUVECs is shown in the histogram. The results are shown as the percentage of inhibition, and three independent experiments were performed; (B) A representative fitting curve from one independent experiment is shown, and the resulting IC50 value was 103.2 nM. HUVECs were plated in 96-well plates in the presence or absence of Ech and allowed to grow for 48 h. The cells were incubated with 20 ng/mL VEGF and then treated with 50, 100, 500, 1000, and 5000 nM Ech. The numbers of cells were quantified using crystal violet stain, and the absorbance was determined by an ELISA reader at 600 nm.
Figure 6. Inhibition of VEGF-induced HUVEC proliferation by Ech. (A) Quantitative analysis of HUVECs is shown in the histogram. The results are shown as the percentage of inhibition, and three independent experiments were performed; (B) A representative fitting curve from one independent experiment is shown, and the resulting IC50 value was 103.2 nM. HUVECs were plated in 96-well plates in the presence or absence of Ech and allowed to grow for 48 h. The cells were incubated with 20 ng/mL VEGF and then treated with 50, 100, 500, 1000, and 5000 nM Ech. The numbers of cells were quantified using crystal violet stain, and the absorbance was determined by an ELISA reader at 600 nm.
Toxins 12 00709 g006
Toxins 12 00709 g007 550
Figure 7. Inhibition of human glioblastoma U373 MG cell migration by Ech and its K45E mutant. (A) Human glioblastoma U373 MG cells were plated on a transwell filter in the presence or absence of Ech proteins and allowed to migrate for 6 h. The cells were treated with 1.53, 5.1, 15.3, 51, 153, and 1530 nM recombinant Ech proteins. The numbers of migrated cells were quantified by counting stained cells microscopically; (B) The results are expressed as the percentage of inhibition, and the inhibition of U373 MG cell migration by Ech and its K45E mutant was dose dependent. The mean of three independent experiments is shown.
Figure 7. Inhibition of human glioblastoma U373 MG cell migration by Ech and its K45E mutant. (A) Human glioblastoma U373 MG cells were plated on a transwell filter in the presence or absence of Ech proteins and allowed to migrate for 6 h. The cells were treated with 1.53, 5.1, 15.3, 51, 153, and 1530 nM recombinant Ech proteins. The numbers of migrated cells were quantified by counting stained cells microscopically; (B) The results are expressed as the percentage of inhibition, and the inhibition of U373 MG cell migration by Ech and its K45E mutant was dose dependent. The mean of three independent experiments is shown.
Toxins 12 00709 g007
Toxins 12 00709 g008 550
Figure 8. Inhibition of human pancreatic Panc-1 cell migration by Ech and its K45E mutant. (A) Human pancreatic Panc-1 cells were plated on a transwell filter in the presence or absence of Ech proteins and allowed to migrate for 6 h. The cells were treated with 1.53, 5.1, 15.3, 51, 153, and 1530 nM recombinant Ech proteins. The numbers of migrated cells were quantified by counting stained cells microscopically; (B) The results are expressed as the percentage of inhibition, and the inhibition of Panc-1 cell migration by Ech and its K45E mutant was dose dependent. The mean of three independent experiments is shown.
Figure 8. Inhibition of human pancreatic Panc-1 cell migration by Ech and its K45E mutant. (A) Human pancreatic Panc-1 cells were plated on a transwell filter in the presence or absence of Ech proteins and allowed to migrate for 6 h. The cells were treated with 1.53, 5.1, 15.3, 51, 153, and 1530 nM recombinant Ech proteins. The numbers of migrated cells were quantified by counting stained cells microscopically; (B) The results are expressed as the percentage of inhibition, and the inhibition of Panc-1 cell migration by Ech and its K45E mutant was dose dependent. The mean of three independent experiments is shown.
Toxins 12 00709 g008
Table
Table 1. Data collection and structure refinement.
Table 1. Data collection and structure refinement.
Crystals
PDB Code		Echistatin
6LSQ
     Data collection
     Radiation sourceNSRRC TPS05A
     Wavelength (Å)0.99984
     Space groupP3221
     Unit cell parameters
           a (Å)33.79
           b (Å)33.79
           c (Å)120.90
           α90.00
           β90.00
           γ120.00
     Resolution (Å)30.00–1.80 (1.91–1.80)
     No. of reflections7695 (1245)
     Completeness (%)92.5 (90.0)
     I/σ (I)7.88 (2.96)
     CC1/299.1 (95.1)
     Refinement
     Resolution (Å)28.44–1.80
     Reflections (work)7541
     Reflections (free)419
     Rwork (%)21.2
     Rfree (%)25.5
     Geometry deviations
           Bond length (Å)0.009
           Bond angles (°)1.461
     Mean B values (Å2)28.851
     Ramachandran plot 1 (%)
         Favored95.7
         Allowed4.3
1 Categories were defined by MolProbity. Values in parentheses are for the highest-resolution shell.
Table
Table 2. Comparison of the Cα(Ri)-Cα(Di + 2), Cβ(Ri)-Cβ(Di + 2), Cζ(Ri)-Cγ(Di + 2), and Cα(Ri)-Cα(Xi + 3) distances (Å) of echistatin and trimestatin.
Table 2. Comparison of the Cα(Ri)-Cα(Di + 2), Cβ(Ri)-Cβ(Di + 2), Cζ(Ri)-Cγ(Di + 2), and Cα(Ri)-Cα(Xi + 3) distances (Å) of echistatin and trimestatin.
LigandsRGD MotifCα(Ri)-Cα(Di + 2)Cβ(Ri)-Cβ(Di + 2)Cζ(Ri)-Cγ(Di + 2)Cα(Ri)-Cα(Xi + 3)
Echistatin (NMR) 1ARGDDM6.2 ± 0.758.2 ± 1.112.3 ± 0.76.3 ± 0.8
Echistatin (A chain)ARGDDM6.58.612.46.2
Echistatin (B chain)ARGDDM5.86.04.85.4
Trimestatin 2ARGDFP6.79.211.16.8
1 The NMR structure of Echistatin; PDB code: 2Ech. 2 The X-ray structure of Trimestatin; PDB code: 1J2L.
Table
Table 3. Summary of the interactions between Ech chain A and αvβ3 integrin.
Table 3. Summary of the interactions between Ech chain A and αvβ3 integrin.
Echαvβ3
Chain Aαvβ3
R22 N313
R24D150, D218
G25D218
D26 S121, Y122, S123, N215, R216, D217
Y31 D126
H44 Y122
Table
Table 4. Summary of the interactions between Ech chain B and αvβ3 integrin.
Table 4. Summary of the interactions between Ech chain B and αvβ3 integrin.
Echαvβ3
Chain Bαvβ3
R24D150, D218
D26 S121, Y122, S123, N215, R216, D217
M28 D126
Y31 D126
K45 M180
Table
Table 5. Summary of inhibition of cell adhesion by Ech and its C-terminal truncated mutant.
Table 5. Summary of inhibition of cell adhesion by Ech and its C-terminal truncated mutant.
Protein Sequence IC50 (nM)
(Mutant) RGD LoopC-Terminusαvβ3αIIbβ3α5β1αvβ5
Ech23ARGDDM42NPHKGPAT20.7 ± 8.0 51.5 ± 3.9 132.6 ± 15.7 286.4 ± 53.6
Ech (P43∆)23ARGDDM42NP132.0 ± 12.0 362.3 ± 56.32469.5 ± 97.7 3358.3 ± 1040.4
Table
Table 6. Summary of inhibition of cell adhesion by Ech and its C-terminal alanine mutants.
Table 6. Summary of inhibition of cell adhesion by Ech and its C-terminal alanine mutants.
ProteinRGD LoopC-Terminus IC50 (nM)
(Mutant)SequenceSequenceαvβ3αIIbβ3α5β1
Ech23ARGDDM42NPHKGPAT 20.7 ± 8.0 51.5 ± 3.9 132.6 ± 15.7
Ech (P43A)23ARGDDM42NAHKGPAT13.8 ± 1.1 67.7 ± 35.0 58.5 ± 9.5
Ech (H44A)23ARGDDM42NPAKGPAT 16.4 ± 6.9 20.7 ± 4.5 30.0 ± 5.3
Ech (K45A)23ARGDDM42NPHAGPAT14.7 ± 4.0132.5 ± 53.384.9 ± 24.9
Table
Table 7. Summary of inhibition of angiogenesis activity by Ech and other disintegrins.
Table 7. Summary of inhibition of angiogenesis activity by Ech and other disintegrins.
Protein RGD Motif C-Terminal Sequence IC50 (nM)
Echistatin ARGDDM NPHKGPAT103.2
Saxatilin ARGDDM NPFHA100.0
Salmosin ARGDDL NPFHA130.0–270.0
Triflavin ARGDFP WNGL100.0–400.0
Rhodostomin PRGDMP YH108.0
Table
Table 8. Summary of inhibition of tumor cell migration by Ech and its C-terminal mutants.
Table 8. Summary of inhibition of tumor cell migration by Ech and its C-terminal mutants.
ProteinSequence IC50 (nM)
RGD LoopC-TerminusA375U373MGPanc-1
Ech23ARGDDM42NPHKGPAT1.5 5.7154.5
Ech (K45E)23ARGDDM42NPHEGPAT13.3 27.21603.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Y.-C.; Chang, Y.-T.; Chen, C.-Y.; Shiu, J.-H.; Cheng, C.-H.; Huang, C.-H.; Chen, J.-F.; Chuang, W.-J. Structural Insight into Integrin Recognition and Anticancer Activity of Echistatin. Toxins 2020, 12, 709. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12110709

AMA Style

Chen Y-C, Chang Y-T, Chen C-Y, Shiu J-H, Cheng C-H, Huang C-H, Chen J-F, Chuang W-J. Structural Insight into Integrin Recognition and Anticancer Activity of Echistatin. Toxins. 2020; 12(11):709. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12110709

Chicago/Turabian Style

Chen, Yi-Chun, Yao-Tsung Chang, Chiu-Yueh Chen, Jia-Hau Shiu, Chun-Ho Cheng, Chun-Hao Huang, Ju-Fei Chen, and Woei-Jer Chuang. 2020. "Structural Insight into Integrin Recognition and Anticancer Activity of Echistatin" Toxins 12, no. 11: 709. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12110709

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