An alternative conformation of human TrpRS suggests a role of zinc in activating non-enzymatic function.

Tryptophanyl-tRNA synthetase (TrpRS) in vertebrates contains a N-terminal extension in front of the catalytic core. Proteolytic removal of the N-terminal 93 amino acids gives rise to T2-TrpRS, which has potent anti-angiogenic activity mediated through its extracellular interaction with VE-cadherin. Zinc has been shown to have anti-angiogenic effects and can bind to human TrpRS. However, the connection between zinc and the anti-angiogenic function of TrpRS has not been explored. Here we report that zinc binding can induce structural relaxation in human TrpRS to facilitate the proteolytic generation of a T2-TrpRS-like fragment. The zinc-binding site is likely to be contained within T2-TrpRS, and the zinc-bound conformation of T2-TrpRS is mimicked by mutation H130R. We determined the crystal structure of H130R T2-TrpRS at 2.8 Å resolution, which reveals drastically different conformation from that of wild-type (WT) T2-TrpRS. The conformational change creates larger binding surfaces for VE-cadherin as suggested by molecular dynamic simulations. Surface plasmon resonance analysis indicates more than 50-fold increase in binding affinity of H130R T2-TrpRS for VE-cadherin, compared to WT T2-TrpRS. The enhanced interaction is also confirmed by a cell-based binding analysis. These results suggest that zinc plays an important role in activating TrpRS for angiogenesis regulation.

Introduction

As a member of class I aminoacyl-tRNA synthetase (aaRS), tryptophanyl-tRNA synthetase (TrpRS) uses a Rossmann-fold catalytic domain to catalyze the attachment of tryptophan to its cognate tRNA, with the help of an anticodon binding domain.[1] In vertebrates, TrpRS have an N-terminal extension of about 150 amino acids beyond the core catalytic architecture.[2,3] The extension is composed of a vertebrate-specific extension (VSE), also known as WHEP domain (named after the tRNA synthetases (i.e., tryptophanyl (W), histidyl (H), and glutamyl-prolyl (EP) tRNA synthetases) first identified to contain this domain), followed by an eukaryotic specific extension (ESE)[4] (Fig. 1A). The WHEP domain has a conserved helix-turn-helix conformation and is dispensable for aminoacylation but critical for regulating the non-enzymatic function of TrpRS in angiogenesis.[2,3] The natural splice variant mini-TrpRS, which lacks most of the WHEP domain but contains the ESE, was found to inhibit VEGF induced angiogenesis.[5,6] Further truncation by leukocyte elastase to remove the first 93 amino acid creates T2-TrpRS that is completely inactive for aminoacylation, but has enhanced anti-angiogenic activity compared to mini-TrpRS.[5]

Figure 1.

Zinc induces conformational change of human FL-TrpRS and the generation of T2-TrpRS. (A) Domain organization of human FL-TrpRS and its crystal structure in complex with Trp-AMP (PDB ID 1R6T). B, C) Time course of plasmin digestion of FL-TrpRS in the absence (B) and presence (C) of zinc (0.1 mM), followed by mass spectrometry analysis. After plasmin digestion, fragments were separated by SDS-PAGE, the bands were cut from the gel and digested with trypsin and AspN for LC-MS/MS analysis. The whole protein mass measurement was performed by MALDI_TOF MS analysis. The cleavage site of the fragments were confirmed by N-terminal sequencing, and the expected mass were calculated from the protein sequence of these fragments. A red box highlights two T2-TrpRS-like fragments generated in the presence of zinc.

T2-TrpRS was shown to inhibit angiogenesis through an interaction with VE-cadherin, which is an endothelial cell junction molecule that adheres apposed endothelial cells through a “strand-swapped” dimer formed by reciprocal swapping of two tryptophan residues of one VE-cadherin molecule into a hydrophobic pocket of the first extracellular cadherin (EC) domain of another VE-cadherin.[7,8,9] Our previous study has demonstrated that T2-TpRS binds the two tryptophan residues (Trp2 and Trp4) through its tryptophan and AMP binding pockets in the active site, respectively, and thereby inhibits the VE-cadherin mediated endothelial cell-cell junctions and disrupts the formation of new blood vessels.[8] Notably, human TrpRS manifests its anti-angiogenic activity only when the WHEP domain is removed. Full length (FL) TrpRS cannot inhibit angiogenesis,[5,6,10] due to the steric hindrance of the WHEP domain that blocks the access of Trp2 and Trp4 of VE-cadherin to the active site pocket.[8]

Zinc is known to have anti-angiogenic effects, which are mediated through multiple mechanisms, including its binding to transcription factors to influence gene expression.[11,12] For example, zinc can dramatically reduce the expression of pro-angiogenic growth factors such as IL-6, IL-8, VEGF and MMP-9, while promoting the production of anti-angiogenic factors such as endostatin.[13,14,15] Zn2+ can also directly bind to endostatin and the binding is essential for the anti-angiogenic activity of endostatin.[16] Interestingly, both bovine and human TrpRS were shown to bind zinc to enhance their aminoacylation activities.[17,18,19] However, no zinc was found in the crystal structures of TrpRS,[4,10,20,21,22,23] suggesting that the zinc-bound conformation may be different. Interestingly, Wakasugi has identified a single amino acid substitution in TrpRS (H130R) that eliminates zinc-induced stimulation and renders constitutively high enzymatic active,[24] suggesting H130R TrpRS may mimic a zinc-bound conformation.

We sought out to investigate the effect of zinc on TrpRS. We found that zinc can induce conformational change of TrpRS and facilitate the proteolytic production of T2-TrpRS. We also found that the H130R mutation induces an alternative conformation of T2-TrpRS that may mimic the zinc-bound conformation. We determined the crystal structure of H130R T2-TrpRS and revealed conformational changes that correlate with an enhanced VE-cadherin binding. Altogether, these results suggest that zinc can promote the production and the anti-angiogenic activity of T2-TrpRS.

Results

Zinc induces conformational changes in human TrpRS and facilitates proteolytic generation of T2-TrpRS

To probe the effect of zinc on the conformation of human TrpRS, we carried out limited protease digestion with plasmin in the presence and absence of zinc, followed by mass spectrometry analysis (Fig. 1B, C). In the absence of zinc, although some proteolytic fragments were generated, the FL-TrpRS (expressed with a C-terminal His-tag to give a total of 484 residues) is relatively stable during 22 hrs' incubation with plasmin (Fig. 1B). However, when in the presence of zinc (0.1 mM), majority of the TrpRS protein were digested by plasmin within 22 hrs, producing a large number of proteolytic fragments (Fig. 1C). This increased susceptibility to protease digestion is indicative of a more relaxed structure of TrpRS induced by zinc binding. Interestingly, two new and major fragments (Gly97-His484, Gly97-Lys472) generated in the presence of zinc were similar to T2-TrpRS (Ser94-Gln471) (Fig. 1C). Therefore, zinc facilitates the generation of a T2-TrpRS-like fragment, most likely through its conformational opening effect on TrpRS.

H130R cannot induce conformational change in FL-TrpRS

Both zinc and the H130R substitution can enhance the aminoacylation activity of human TrpRS,[24] suggesting H130R may mimic a zinc-bound conformation in the absence of zinc. To investigate the potential conformational change induced by H130R, we subjected WT and H130R FL-TrpRS to trypsin digestion in the absence of zinc; however, no obvious difference in the digestion pattern between the two proteins was detected (Fig. 2A). Also, no obvious difference between the two proteins was detected in the presence of zinc (Fig. 2B). Consistent with the previous result (Fig. 1B, C), TrpRS is more efficiently digested in the presence of zinc, and a T2-TrpRS-like fragment can be generated (Fig. 2B). (To rule out the possibility that the enhanced efficiency of digestion is due to enzymatic activity enhancement of trypsin in the presence of zinc, we used another tRNA synthetase as a control. No change in the digestion pattern of ProRS was observed with and without zinc (Fig. S1)). These results further support the idea that zinc induces a more relaxed conformation in FL-TrpRS; however, this conformational change cannot be induced by the H130R substitution at least under our experimental conditions.

Figure 2.

Trypsin digestion analyses suggest H130R mutation induces a zinc-bound conformation in the context of T2-TrpRS. (A, B) Trypsin digestion of wild-type (WT) and H130R FL-TrpRS with increasing concentrations of the protease in absence (A) and presence (B) of zinc (0.1 mM). The concentration of trypsin is labeled as [protease] to [TrpRS] ratio. (C, D) Trypsin digestion of wild-type (WT) and H130R T2-TrpRS in absence (C) and presence (D) of zinc (0.1 mM). Red boxes highlight the resemblance of the digestion pattern of H130R T2-TrpRS in the absence of zinc with that of WT T2-TrpRS in the presence of zinc.

We were able to obtain crystals of H130R FL-TrpRS bound with Trp-AMP but not crystals of apo H130R FL-TrpRS. The structure of the complex was solved at 2.1 Å resolution (Table S1). Consistent with the result of the trypsin digestion experiment, the complex structure resembles that of wild-type FL-TrpRS in complex with Trp-AMP (PDB ID 1R6T), with a small root-mean-square deviation (RMSD) value of 0.64 Å for the main chain Cα atoms (Fig. S2). Therefore, we conclude that H130R substitution does not induce significant conformational change in FL-TrpRS, at least when it is bound with Trp-AMP.

H130R mutation mimics zinc and induces conformational change in T2-TrpRS

Next, we tested the effect of the H130R substitution in the context of T2-TrpRS. In the presence of zinc, the trypsin digestion pattern is the same between WT T2-TrpRS and the H130R mutant (Fig. 2D). However, in the absence of zinc, there is a difference between WT and H130R T2-TrpRS in their fragmentation pattern at low trypsin concentration (trypsin to TrpRS concentration ratio as 1:100,000 and 1:10,000) (Fig. 2C). Remarkably, the pattern of H130R T2-TrpRS in the absence of zinc resembles that of WT T2-TrpRS in the presence of zinc (Fig. 2C, D), suggesting that in the context of T2-TrpRS the H130R substitution does induce conformational change, which mimics the effect of zinc.

Overall conformational changes of H130R T2-TrpRS

To visualize the conformational change, we determined the crystal structure of H130R T2-TrpRS at 2.8 Å resolution by molecular replacement using WT human T2-TrpRS structure (PDB 2QUI) as the search model.[25] Like the wild-type protein (Fig. 3A), H130R T2-TrpRS is composed of a truncated ESE (Gly97-Lys154), a Rossmann fold catalytic domain (Pro155-Ser353) and an anticodon binding domain (Asp354-Gln471) (Fig. 3B). A homodimer of H130R T2-TrpRS was found in the asymmetric unit, with two virtually identical subunits A and B interacting at the catalytic domain to form a “W” form dimer (Fig. 3B). The structure is well folded overall with clear electron density, except for the disordered regions Gln344-Pro355 and Gln389-Gly392 in subunit A, Met350-Ser353 and Gly379-Gly391 in subunit B. Superposition of the structure with that of WT T2-TrpRS (PDB 2QUI) gives a large main chain RMSD value of 2.191 Å. Major conformational changes were observed at the entire ESE, KMSAS signature motif containing loop (Thr338-Thr362), and loop Ala376-Asp396 in the anticodon binding domain (Fig. 3A, B).

Figure 3.

Conformational change induced by H130R substitution in T2-TrpRS. (A) Overall crystal structure of wild-type (WT) T2-TrpRS homo-dimer. The eukaryotic specific extension (ESE, Glu82-Lys154) and the Ala376-Asp397 loop in the anticodon binding domain are colored in blue, whereas the KMSAS motif containing loop (Thr338-Thr362) is colored in green. Secondary structural elements of the ESE in subunit A are labeled. (B) The crystal structure of H130R T2-TrpRS homo-dimer. The ESE (Gly97-Lys153) and the Ala376-Asp397 loop are colored in magenta, whereas the KMSAS loop (Thr338-Thr362) is in orange. (C) Stereo view of the ESE structure in H130R (magenta) versus WT (blue) T2-TrpRS. (D) Stereo view of the conformational changes surround active site. WT T2-TrpRS in complex with Trp-AMP (PDB 1R6U) is aligned with H130R mutant structure. The α11, α12 helices and the Ala376-Asp397 loop are colored in blue in WT structure and in magenta in the H130R structure. KMSAS loop (Thr338-Thr362) is colored in green in WT and in orange in the H130R mutant.

The ESE undergoes dramatic structural changes (Fig. 3A-C). Compared to the wild-type structure, the α1′ helix in the H130R mutant no longer covers the substrate binding pocket, and is rotated 90° and shifted away about 13.7 Å from the catalytic domain to the anticodon binding domain. Following the α1′ helix, α2′ helix is also rotated for 37°, while β1′, β2′ strands, and α3′ helix all become disordered. (The H130R substitution site is located in α3′; see below.) The α4′ helix is the last part of the ESE and hinged for 34° from the end of the helix (Fig. 3C).

The catalytic domain of the H130R mutant superimposes well with that of the WT T2-TrpRS, except for the loop containing the KMSAS motif. HIGH and KMSKS (170HVGH173 and 349KMSAS353 in human TrpRS) are two conserved sequence motifs of class I tRNA synthetases involved in ATP recognition and catalysis; a third and less conserved motif in the catalytic domain (310AIDQ313 in human TrpRS) also participates in ATP binding.[22,23,25] The KMSAS region (Thr338-Thr362) is resolved in the WT T2-TrpRS crystal structure as flexible loops and turns;[4,22,25] however, it becomes completely disordered in subunit A of H130R T2-TrpRS, while in subunit B it folds into an α helix (Thr338-Ser351) (Fig. 3A, B). As a result of these conformational changes, the Trp-AMP binding pocket is more “open” in the H130R mutant (Fig. S3).

In the structure of the WT T2-TrpRS, α11 is located at the tip of anticodon binding domain and is involved in binding the tRNA anticodon loop.[4,22,26] However, in the H130R T2-TrpRS structure, it is folded as a loop (Ala376-Asp397) and flipped 180° onto the active site pocket of subunit A, while this loop is disordered from Gly379 to Gly391 in subunit B (Fig. 3B). Most interestingly, the Ala376-Asp397 loop in subunit A occupies a similar location as that of the KMSAS loop (Thr338-Thr362) in the WT T2-TrpRS, partially compensating for the disordered KMSAS loop in the mutant (Fig. 3D). (α11 is disordered in subunit B of the mutant.) Following α11, helix α12 swings ∼36° towards the active site pocket in both subunits (Fig. 3D).

Conformational changes around the H130R substitution site

In the structure of WT T2-TrpRS, His129, His130 and Phe131 are neighboring residues located on α3′ helix of the ESE (Fig. 4A). Particularly, NE2 of His130 forms a weak hydrogen bond with the main chain carbonyl oxygen of Phe405 (3.12 Å); side chain of Phe131 and His129 protrude into a hydrophobic patch formed by Ile175, Ile178, Phe179, Trp182 from α1 and Phe405, Phe406 from α12 of the anticodon binding domain. In addition, the α3′ helix is further stabilized by interactions with residues from α12 (Asp409-Ser422) and α13 (Asp397-Leu407). For example, Arg127 located on the loop between α2′ and α3′ interacts with Asp410 through a salt bridge (2.92 Å), whereas Arg134 on α3′ forms a hydrogen bond (2.56 Å) with Leu407 located on α12 (Fig. 4A).

Figure 4.

Local conformational changes around the H130R substitution site. (A, B) Conformational changes around residue 130 in WT (A) and H130R (B) T2-TrpRS. (C, D) Interactions between the ESE and catalytic domain in WT (C) and H130R (D) T2-TrpRS.

In contrast, H130R T2-TrpRS shows strikingly different local conformation around the substitution site (Fig. 4B). The α3′ helix becomes unwound. The Cα position of the 130 residue shifts 6.94 Å; the neighboring Phe131 is shifted 3.5 Å. More dramatically, the Cα position of His129 moves 5.74 Å, and its imidazole side chain flips 180° to become solvent exposed. Arg127 is no longer salt bridged with Asp410 and a large gap (10.58 Å) is opened up in between the two residues. Arg134 is also moved away from Leu407 (Fig. 4A, B).

These local conformational changes abolish the interactions of the α2′ and α3′ helices from the ESE with α1 of the catalytic domain and α12 and α13 from the anticodon binding domain (Fig. 4A, B). Many of the secondary structures such as β1′, β2′ and α3′ in the ESE become disordered. The side chains of Phe137 and Phe138 are flipped 90° and 180°, respectively, while the side chain of Arg141 is also flipped away from the catalytic domain and becomes solvent exposed (Fig. 4C, D). As a result, the entire ESE dissociates away from the catalytic domain. These and other conformational changes within the catalytic domain, make the Trp-AMP binding pocket more open in H130R T2-TrpRS (Fig. 4C, D; Fig. S3).

Molecular simulations of the interaction between T2-TrpRS and VE-cadherin

Considering the importance of the Trp and AMP binding pockets for the T2-TrpRS and VE-cadherin interaction,[8] we evaluated how the H130R mutation would affect the interaction using molecular docking and dynamic simulations (MDS). We used the crystal structure of the EC1-2 domain of chicken VE-cadherin[27] (PDB ID 3PPE) that share 58% sequence identity with human VE-cadherin to docked onto WT and H130R T2-TrpRS structures, respectively. Considering the expected low concentration of the extracellular T2-TrpRS, we used the monomer form of T2-TrpRS for the docking, followed by dynamic simulations to obtain final stable complexes (Fig. 5). As suggested by the small fluctuations of root-mean-square deviations (RMSDs) during the MDS process, the docking of WT and H130R T2-TrpRS are rather stable (Fig. S4A, B).

Figure 5.

Molecular dynamic simulations (MDS) of the interaction between extracellular domain EC1-2 of VE-cadherin and the monomer form of WT (A) and H130R (B) T2-TrpRS. The eukaryotic specific extension (ESE) and the Ala376-Asp397 loop in the anticodon binding domain of T2-TrpRS are colored in blue (WT) and magenta (H130R), whereas the KMSAS loop (Thr338-Thr362) is colored in green (WT) and orange (H130R). The EC1-2 domain of VE-cadherin is colored in yellow, the Trp2 and Trp4 residues are shown in sticks.

Interestingly, the docking model of WT and H130R T2-TrpRS are dramatically different. As shown in Fig. 5A and B, residues Trp2 and Trp4 of EC1-2 are docked in the Trp and AMP binding pocket of the WT T2-TrpRS, whereas the two Trp residues are no longer inserted in the Trp and AMP binding pockets of H130R T2-TrpRS. However, the contact areas between H130R T2-TrpRS and EC1-2 (1447.45 Å) are substantially larger than that in the WT complex (1276.4 Å). The EC1-2 domain adopts an extended conformation in the WT complex, while the two domains are more compact and has a different orientation in the H130R complex.

H130R mutation enhances the interaction of T2-TrpRS with VE-cadherin

The larger contact areas in the H130R complex suggest that the mutation may enhance the T2-TrpRS and VE-cadherin interaction. We test the interaction by Surface Plasmon Resonance (Fig. 6A, B). The EC1-2 domain of VE-cadherin was expressed from baculovirus transfected insect cells. Purified EC1-2 domain of VE-cadherin was immobilized on a CM5 chip, and the E. coli expressed human WT or H130R T2-TrpRS was injected as analyte. Interestingly, the binding affinity of the H130R mutant is more than 50-fold stronger than WT T2-TrpRS for EC1-2 (Fig. 6A, B), and the binding can be detected only in the presence of 0.1 mM ZnCl2 (see discussion). The stronger affinity of H130R T2-TrpRS is mostly due to a faster association rate (ka of 2.71 × 107 M−1s−1 for H130R versus 1.76 × 105 M−1s−1 for WT), consistent with the model that the H130R structure provides a larger binding surface for VE-cadherin.

Figure 6.

H130R substitution enhances the interaction between T2-TrpRS and VE-cadherin. (A) The fitted sensorgram of the binding between wild-type (WT) T2-TrpRS and EC1-2 domain of VE-cadherin by surface plasmon resonance. The response unit is plotted over 400 s that illustrates the association and dissociation of the complex. The concentration of WT T2-TrpRS is 2-fold increased from 12.5 nM to 200 nM. (B) The fitted sensorgram of the binding between H130R T2-TrpRS and EC1-2 domain of VE-cadherin. The response unit of H130R T2-TrpRS at concentration from 0.78 nM to 25 nM is plotted over 400 s. (C) In vitro pull down assay to test the interaction between TrpRS and endogenous VE-cadherin. The melanoma C8161 cells expressed VE-cadherin was pulled down by 6xHis tagged FL- and T2-TrpRS (WT and H130R), and immune-blotted by rabbit anti VE-cadherin. The TrpRS variants are blotted with HRP-conjugated 6xHis antibodies.

To further test the effect of H130R on VE-cadherin binding, we used a cell-based assay. Human melanoma C8161 cells express VE-cadherin and were used in a pull-down assay to detect the amount of the endogenous VE-cadherin that can be pulled down by various human TrpRS proteins (Fig. 6C). As expected from previous studies,[8] T2-TrpRS pulled down substantially larger amounts of VE-cadherin than did the FL-TrpRS. Interestingly, the H130R mutation enhanced the VE-cadherin interaction, especially in the context of T2-TrpRS. Notably, the H130R mutation also enhanced the VE-cadherin interaction with FL-TrpRS. Possibly, the presence of an interaction partner helps H130R to manifest a conformational opening effect on FL-TrpRS.

Discussion

Zinc has long been implicated in the function of mammalian TrpRS.[18] Bovine TrpRS was found to be zinc-bound, and elimination of zinc affected the tertiary structure and abolished the enzymatic activity.[18] Later studies using human TrpRS confirmed the zinc-binding capacity and the effect of zinc binding in stimulating the synthetase activity.[19] Interestingly, bovine TrpRS differs from human TrpRS at position 130 (according to the human sequence): bovine TrpRS has Arg instead of His as found in human TrpRS at this position. Remarkably, the H130R substitution in human TrpRS abolished zinc-induced stimulation and provided constitutively high enzymatic activity,[24] suggesting that human H130R TrpRS mimics a zinc-bound form, while the Arg130 residue in bovine TrpRS is not sufficient to induce the zinc-bound conformation. Therefore, it seems that human and bovine TrpRS have different affinity or sensitivity towards zinc and the position 130 is a determinant for this difference.

The fact that bovine TrpRS (with Arg130) can still bind zinc[18] and undergoes conformational change indicates that the H130R mutant may represent an intermediate state between zinc-free and zinc-bound conformations in human TrpRS. The conformation of this intermediate state may be shifted toward the zinc-free form by binding to Trp-AMP and towards the zinc-bound form by removal of the N-terminal residues. This speculation helps to explain the lack of conformational change in full-length human H130R TrpRS both in solution (Fig. 2) and in crystals (Fig. S2), and to interpret the observation that the enhanced VE-cadherin binding by H130R mutation in the context of T2-TrpRS is only observed in the presence of zinc (Fig. 6). It is possible that under certain conditions (such as the condition used for our SPR experiment) the H130R substitution alone is not sufficient to induce the zinc-bound conformation (like in bovine TrpRS); nevertheless, it can facilitate the conformational change induced by zinc binding. We should note that we have attempted to crystallize WT FL- and T2-TrpRS in the presence of zinc but did not yield any crystal. This suggests that the zinc-bound conformation is not as stable as the zinc-free conformation, which is consistent with our limited protease digestion analyses (Fig. 1B, C and Fig. 2).

Our previous work has demonstrated the critical role of Trp and AMP binding pockets in T2-TrpRS for VE-cadherin binding.[8] However, in the H130R mutant structure, large conformational changes occur around the active site pockets. In the docking model (Fig. 5), Trp2 and Trp4 of VE-cadherin are no longer bound to the active site of H130R T2-TrpRS. It remains an open question for whether the active site of T2-TrpRS is still involved in the VE-cadherin interaction upon zinc binding.

A recent paper reported that secreted TrpRS is involved in priming innate immunity upon pathogen infection through mediating an interaction between toll-like receptor 4 (TLR4) and myeloid differentiation factor 2 (MD2).[28] Interestingly, the activity of the FL-TrpRS can be recapitulated by the N-terminal 154 residues, which contain the WHEP domain and the ESE (Fig. 1A). In particular, it was proposed that the last helix of ESE is involved in the TLR4 interaction. As shown in Fig. 3, the ESE including the last helix of ESE (α4′) undergoes large conformational change in the H130R mutant that mimics a zinc-bound conformation. Therefore, it is possible that zinc binding may also influence the role of TrpRS in priming innate immunity.

Overall, our work has provided direct evidence for an alternative conformation of human TrpRS that mimics a zinc-bound conformation. This conformation may enhance the role of TrpRS in angiogenesis regulation by facilitating the proteolytic production of the anti-angiogenic fragment T2-TrpRS and by strengthening the interaction between T2-TrpRS and its receptor VE-Cadherin. For the first time, this work links zinc to regulatory functions of tRNA synthetases.

Accession numbers

The atomic coordinates and structure factors of FL-TrpRS in complex with Trp-AMP and the H130R T2-TrpRS have been deposited in the Protein Data Bank with the accession code 5UJI and 5UJJ, respectively.

Material and methods

Protein expression and purification

The genes encoding human H130R FL- and T2 (Gly97 to Gln417) TrpRS were amplified by PCR and constructed into pET20b vector (Novagen) with a C-terminal 6 × His tag. The proteins were expressed in transformed Escherichia coli BL21 (DE3) cells and induced by 0.1 mM Isopropyl β-D-1-thiogalactopyranoside (IPTG) at room temperature for 8 hrs. The cells were disrupted by M-110L Microfluidizer Processor (Microfluidics) in 25 mM Hepes-Na pH 7.5, 300 mM NaCl, 20 mM imidazole and 5% glycerol, the cell lysates were centrifuged at 25,000 rpm for 30 min. The soluble fraction was applied on a Ni-NTA agarose column (Qiagen) for affinity purification. T2-TrpRS was eluted by 250 mM imidazole and further purified by ion-exchange chromatography by Resouece Q column (GE Healthcare), protein was eluted by 25 mM Hepes-Na pH7.5 and sodium chloride gradient to 1 M. The fractions were concentrated and applied to HiLoad 16/60 Superdex 200 prep grade column (GE Healthcare) for size exclusion purification, purified human H130R T2-TrpRS was concentrated to 10 mg mL−1 and set up for crystallization. H130R FL-TrpRS was purified by same procedure as H130R T2-TrpRS.

Crystallization

Crystallization was performed by micro-sitting drop vapor diffusion method using Mosquito liquid transfer robot (TTP Labtech). H130R T2-TrpRS was pre-incubated with 2 mM MgCl2 and 4 mM β-Mercaptoethanol for crystallization. Crystals were obtained by mixing 0.1 uL of protein with 0.1 uL reservoir solution (20% PEG4000, 0.1 M Na Citrate pH5.0, 0.1 M MgCl2) equilibrated against 70 uL reservoir solution at 4°C.

H130R FL-TrpRS was diluted in a buffer containing 25 mM Hepes-Na pH7.5, 100 mM KCl, 2 mM MgCl2 and 4 mM β-Mercaptoethanol, and pre-incubated with 500 uM Trp-AMP for crystallization. The crystals were obtained in 1:1 ration of the protein and reservoir solution containing 20% PEG1500, 0.1 M Bis-Tris pH6.5.

Data collection

The crystals were cryoprotected in a solution containing the reservoir solution with 15% glycerol and flash-cooled with liquid nitrogen. A 2.8 Å resolution data set of human H30R T2-TrpRS and 2.1 Å resolution data set of H30R FL-TrpRS were collected at 100 K on beamline BL7-1 at the Stanford Synchrotron Radiation Laboratory using an ADSC Q315 detector. The crystal of H130R T2-TrpRS belongs to the space group C2221, with unit cell dimensions a = 84.89 Å, b = 139.06 Å, c = 146.72Å, α = β = γ = 90° (Table 1). Whereas the crystal of H130R FL-TrpRS in complex with Trp-AMP belongs to space group C121, and with unit cell parameters a = 136.35, b = 95.28, c = 99.63, α = γ = 90°, β = 130° (Table S1). Processing of all diffraction images and scaling of the integrated intensities were performed using the HKL2000.[29] There are two molecules of H130R T2-TrpRS in one asymmetric unit with Matthews coefficient VM of 2.5 Å3 Da−1, corresponding to a solvent content of 51%.[30]

Table 1.

X-ray crystallography data collection and refinement statistics.

Parameters	T2-TrpRS H130R
Data collection
Cell parameters (Å)	a = 84.89, b = 139.06,
	c = 146.72
	α = β = γ = 90°
Space group	C2221
Resolution (Å)	50 (2.85) – 2.80
No. of all reflections	108967
No. of unique reflections	21588
Completeness (%)	98.9 (99.4)
Redundacy	5.0 (5.1)
I/σI	6.5 (2.5)
Rmerge (%)	9.3 (47.1)
Refinement
Resolution (Å)	50 – 2.79
Total No. of reflections	20462
No. of reflections used	19358
Rwork / Rfree (%)	22.0 / 26.2
No. of atoms (Protein)	5762
R.m.s. deviations
Bond lengths (Å)	0.007
Bond angle (°)	1.15
Average B-factors (Å2)	73.0
Ramachandran Plot
Most favored region (%)	96.0
Allowed regions (%)	4.0
Disallowed regions (%)	0

Note: Values in parentheses are for the highest resolution shell.

Structure determination

The structure of H130R T2-TrpRS was determined by molecular replacement using wild-type human T2-TrpRS structure (PDB 2QUI)[25] as search model in program PHASER[31] from the CCP4 package.[32] One homodimer of H130R T2-TrpRS was found in one asymmetric unit. Iterative model building and refinement were performed using Coot,[33]and Refmac5,[34] the final model with R of 22.0% and R of 26.2% at 2.8 Å resolution. The H130R FL-TrpRS structure was determined using the structure of WT human FL-TrpRS (PDB 1R6T) as serch model and the method molecular replacement. Final model was refined to R of 21.5% and R of 23.2% (Table S1).

Structure and molecular docking

To derive the binding conformations of the wild-type (PDB ID 2QUI) and H130R T2-TrpRS with EC1-2 of VE-cadherin (PDB ID 3PPE), molecular docking was performed using ZDOCK (version 3.0.2) software.[35] Using a fast Fourier transform (FFT)-based initial-stage rigid-body molecular docking algorithm,[36] ZDOCK performs a global search in the translational and rotational space to produce all possible binding configurations between proteins. With a composite scoring function that combines pairwise shape complementarity with desolvation and electrostatics, the performance of ZDOCK in critical assessment of prediction of interaction (CAPRI) challenge proves that it is among the best protein-protein docking algorithms.[37,38] ZDOCK created 3600 and retained 2000 putative binding configurations of WT / H130R T2-TrpRS and EC1-2, the complex models with the top default ZDOCK scores were chosen to be the most possible binding configurations. To remove bad atomic contacts around the binding interface, each binding model was then subjected to 2000 steps steepest descent and 2000 steps conjugate gradient energy minimizations using SANDER module in the AMBER 12 suit of programs.[39]

Molecular dynamics simulation

Molecular Dynamics simulation (MDS) was performed using the AMBER 12 software and FF14SB force field.[39,40] The complex models of WT / H130R T2-TrpRS bound with EC1-2 domain of chicken VE-cadherin were independently immersed into the center of a truncated octahedron box of TIP3P water molecules with a margin distance of 12.0 Å, potassium counterions were added by using the AMBER XLEAP module to keep system in electric neutrality.[39] Each binding complex was firstly energy minimized by the steepest descent method for 2000 steps with WT / H130R T2-TrpRS and EC1-2 restricted by a harmonic constraint of 100 kcal·mol−1Å−2. A further conjugate gradient minimization of 2000 steps was performed with no constraint. Then the system was gradually heated from 0 K to 300 K under the NVT ensemble over a period of 200 ps, during which the Langevin thermostat with a coupling coefficient of 1.0 ps and a weak constraint of 10 kcal·mol−1Å−2 on proteins was applied. Each model was subsequently subjected to an equilibrium simulation for 200 ps by removing all constraints. Finally, a 20 ns production MDS for each model was conducted under NPT ensemble. Periodic boundary conditions were applied. System temperature was kept 300 K using the Berendsen thermostat with a time constant of 1 ps. Isotropic constant pressure was maintained by Berendsen pressure coupling algorithm with a time constant of 1 ps. Hydrogen involved covalent bonds were constrained by the SHAKE algorithm.[41] The long-range electrostatic interactions were treated by the Particle Mesh Ewald (PME) method.[42] The cutoffs for long-range electrostatic and Van der Waals interactions were both set to 12.0 Å. The time step was set to 2 fs in all simulations, and the coordinates were saved every 1 ps to record the MD trajectories.

Pull-down assays

Lysates of E. coli cells that transformed with genes expressing 6 × His tagged full-length or T2-TrpRS (WT and H130R) were mixed with lysates from melanoma C8161 cells expressing VE-cadherin. The mixtures were incubated at 4°C overnight, followed by incubation with Ni-NTA resin for 1 hour to capture histidine-tagged TrpRS proteins and VE-cadherin complex. Input samples were taken out prior to the addition of Ni-NTA resin. The resin was washed 3 times with wash buffer (20 mM Tris-HCl, pH 7.5, 150 mM NaCl, 1 mM DTT, 20 mM imidazole). VE-cadherin pull downed by TrpRS was eluted with SDS-PAGE reducing sample buffer and subjected to SDS-PAGE followed by western blot with rabbit anti VE-cadherin (Cell Signaling Technology) and HRP-conjugated 6xHis (Proteintech) antibodies to detect VE-cadherin and TrpRS proteins, respectively.

Surface plasmon resonance

The EC1-2 domain of VE-cadherin was expressed in a Bac-to-Bac baculovirus system using Sf9 insect cells. The EC1-2 domain was purified and immobilized on a CM5 chip as bait ligand though amine coupling. The WT or H130R T2-TrpRS expressed and purified from E.coli was injected as prey analyte in HBS-P buffer containing 20 mM Hepes pH7.5, 150 mM NaCl, 0.02% surfactant P-20. The response unit is plotted over injection time to probe the association and dissociation of T2-TrpRS with immobilized EC1-2 domain. Interestingly, the binding only can be detected when 0.1 mM ZnCl2 was present in the HBS-P buffer. After optimizing the association and dissociation parameters, the binding sensorgram, kinetic analysis was performed by injecting the WT or H130R T2-TrpRS at series concentrations. The binding constant was calculated from the sensorgram of WT T2-TrpRS ranges from 12.5 nM to 200 nM and the H130R T2-TrpRS from 0.78 nM to 25 nM in a 2 fold dilution.