Abstract
Ligand-independent activation of VEGFRs is a hallmark of diabetes and several cancers. Like EGFR, VEGFR2 is activated spontaneously at high receptor concentrations. VEGFR1, on the other hand, remains constitutively inactive in the unligated state, making it an exception among VEGFRs. Ligand stimulation transiently phosphorylates VEGFR1 and induces weak kinase activation in endothelial cells. Recent studies, however, suggest that VEGFR1 signaling is indispensable in regulating various physiological or pathological events. The reason why VEGFR1 is regulated differently from other VEGFRs remains unknown. Here, we elucidate a mechanism of juxtamembrane inhibition that shifts the equilibrium of VEGFR1 towards the inactive state, rendering it an inefficient kinase. The juxtamembrane inhibition of VEGFR1 suppresses its basal phosphorylation even at high receptor concentrations and transiently stabilizes tyrosine phosphorylation after ligand stimulation. We conclude that a subtle imbalance in phosphatase activation or removing juxtamembrane inhibition is sufficient to induce ligand-independent activation of VEGFR1 and sustain tyrosine phosphorylation.
Similar content being viewed by others
Introduction
The vascular endothelial growth factor receptors (VEGFR) are the key regulator of normal physiological and pathological angiogenesis and vasculogenesis1,2. The VEGFR family comprises three receptor tyrosine kinases (RTK): VEGFR1, VEGFR2, and VEGFR3. Among them, VEGFR1 is an elusive family member. Even after three decades of its discovery, the function and regulation of VEGFR1 remain poorly understood3,4,5. VEGFR2 is the primary receptor for VEGFs that regulates diverse cellular functions, including blood vessel development during embryogenesis, hematopoiesis, and tumor angiogenesis1,6. VEGFR3, on the other hand, is the primary receptor for the lymphangiogenic factor VEGF-C and VEGF-D2. During embryonic development, VEGFR1 acts as a decoy receptor. The VEGFR1 and VEGFR2 are activated by a common bivalent ligand (VEGF-A). VEGFR1 negatively regulates the VEGFR2 signaling by sequestering excess VEGF-A, preventing over-activation of VEGFR25,7,8 Compared to VEGFR2, VEGFR1 binds to its ligand VEGF-A with a ten-fold stronger affinity9,10. Yet, the ligand binding induces only a weak kinase activation in VEGFR1 and does not generate subsequent downstream signaling in endothelial cells, vascular smooth muscle, or fibroblast cells4,8,11,12. Although VEGFR1 and VEGFR2 share a high degree of sequence and structural homology, it is unclear why the two RTK are differently regulated.
VEGFR1 and VEGFR2 share similar structural architecture, comprising an extracellular ligand-binding domain (ECD) made up of seven immunoglobulin-like subdomains (D1 to D7), a single-passed transmembrane (TM) segment, a cytosolic juxtamembrane (JM) segment tethered to a kinase domain (KD) followed by a C-terminal tail (Fig. 1a)1,2,13. VEGFRs are activated by a canonical model of monomer to dimer (or multimer) transition upon ligand binding14,15,16. In the unligated state, the receptor exists predominantly as a monomer (Fig. 1b)17,18, and the kinase-domain adopts a platelet-derived growth factor receptor (PDGFR)-like JM-in inactive conformation19,20,21,22. In the inactive conformation, the N-terminal portion of the JM segment is buried inside the catalytic site, like a wedge, locking the kinase domain in an autoinhibited conformation. Binding to bivalent ligands leads to a ligand-dependent dimerization of the ECD (Fig. 1b)23,24. In the ligated state, the VEGFR forms a symmetric dimer in the outside (ECD) and an asymmetric dimer inside (kinase domain)23,24,25. The ECD dimerization rearranges the TM segment26, removes the JM-inhibition (to JM-out conformation), and brings two adjacent kinase-domain in close proximity, allowing autophosphorylation of multiple tyrosine residues in the C-terminal tail (Fig. 1b)19,21,23,27. The phosphotyrosine residues then function as a docking site for assembling downstream signaling modules. Structural analysis of the kinase domain suggests that VEGFR1 is not a pseudokinase. All the regulatory motifs (R-spine and C-spine) and the catalytic residues are conserved (Supplementary Fig. 1a)28,29. The lack of kinase activity of VEGFR1 was attributed to an inhibitory sequence in the JM segment30 and Asn1050 in the A-loop31, the molecular mechanism of which is unknown.
In contrast, VEGFR1-signaling is indispensable in regulating hematopoietic cell function and developing pathophysiological conditions5,7. Ligand-dependent activation and VEGFR1-mediated cell-signaling regulate diverse physiological functions32,33,34,35. Overexpression or downregulation of VEGFR1 is linked to several cancers and cancer-associated pain, retinopathy, tumor survival, and autoimmune disorders36,37,38,39,40,41,42. The mechanism of how VEGFR1 autoinhibition is released under pathological conditions is an open question.
To gain further insight, we investigated the ligand-independent and ligand-dependent activation of VEGFR1 and VEGFR2 on the plasma membrane by a single-cell assay using fluorescence microscopy. Our data revealed that, unlike VEGFR2, VEGFR1 does not show concentration-dependent autophosphorylation in the absence of a ligand and is transiently phosphorylated upon ligand stimulation. We decipher that an electrostatic latch in the JM-S and an H-bond between a tyrosine residue in the JM-B and C-helix in VEGFR1 together constitute a JM inhibition that likely stabilizes the inactive JM-in conformation. Slow release of the JM inhibition makes the VEGFR1 autophosphorylation inefficient. Finally, we proposed a mechanism explaining how a delicate balance in releasing JM inhibition maintains the VEGFR1 signaling constitutively off in the unligated state.
Results
VEGFR1 does not show concentration-dependent activation without ligands at the plasma membrane
Ligand-independent activation of RTKs is a key signature of several forms of cancer and manifestation of drug resistance43,44,45,46,47,48,49. Receptor density at the plasma membrane is an important determinant of ligand-independent activation of RTKs50,51,52,53,54,55,56,57. The density-dependent activation of RTK was explained by an equilibrium shift model between multiple receptor species58,59. Recent studies showed that VEGFR2 forms a ligand-independent dimer at a physiological concentration on the membrane and is able to autophosphorylate26. We ask, in the unligated state, if VEGFR1 autophosphorylates spontaneously on the plasma membrane.
We begin with a single-cell assay to comparatively study the concentration-dependent activation of VEGFR1 and VEGFR2 at the plasma membrane with and without ligand stimulation, respectively (Fig. 2, Supplementary Figs. 1, 2). We transiently transfected CHO cell lines with VEGFR1-mCherry or VEGFR2-mCherry constructs and stimulated them with VEGF165. The transient transfection generates a heterogeneous population of cells expressing a diverse concentration of receptors on the plasma membrane. Since the localization of VEGFR family kinases is not solely restricted to the plasma membrane60,61, in this study, we focused on the regions surrounding the plasma membrane (cell perimeter) (Supplementary Fig. 1d). The activation of the receptor at the membrane was probed by determining the phosphorylation level of Y1213 or Y1175 for VEGFR1 or VEGFR2, respectively, with specific antibodies62,63,64. We observed that the unligated VEGFR2 did not autophosphorylate Y1175 at low receptor concentrations but phosphorylates spontaneously at higher receptor concentrations59 (Fig. 2a, b, e). Our data suggests that a critical concentration of the receptors at the plasma membrane is required to activate the kinase domain ligand independently. Whereas, upon ligand stimulation, VEGFR2 linearly phosphorylates Y1175, suggesting the phosphorylation is independent of receptor concentration (i.e., phosphorylation level is proportional to the receptors at the plasma membrane).
In the single-cell assay, the VEGFR1 also linearly phosphorylates Y1213 when stimulated with VEGF165 (Fig. 2d, f). However, the fraction of tyrosine phosphorylated by VEGFR1 is significantly lower than VEGFR2 (Supplementary Fig. 2f). Unexpectedly, we observed that VEGFR1 did not show any ligand-independent autophosphorylation of Y1213, even at the highest receptor concentration measured in our studies (Fig. 2c, f). To rule out if the lack of ligand-independent activation of VEGFR1 is not cell-dependent, we repeat the assay by transiently transfecting VEGFR1 to COS-7 and a macrophage cell line (RAW264.7) (Supplementary Fig. 2a–e). We observed a similar phosphorylation profile of Y1213, as seen in the CHO cell line (Fig. 1f and Supplementary Fig. 2c–d). Suggesting that the lack of ligand-independent activation of VEGFR1 is an intrinsic property of the receptor and not an artifact. We then ask: why does VEGFR1 phosphorylate a lower fraction of tyrosine residues than VEGFR2 upon ligand stimulation? What is the molecular basis that constitutively inactivates VEGFR1 in the ligand-free state?
Phosphorylation of VEGFR1 is transiently stable compared to VEGFR2
To understand why VEGFR1 and VEGFR2 are phosphorylated differentially upon ligand stimulation, we next studied the phosphorylation kinetics and half-life of phosphotyrosine residue (Fig. 2g, h and Supplementary Fig. 2g, h). We determined the phosphorylation level of Y1213 or Y1175 in VEGFR1 or VEGFR2, respectively, by immunoblotting over a period of time after ligand stimulation. The phosphorylation kinetics (Fig. 2h, Supplementary Fig. 2g, and Supplementary Table 4) shows that VEGFR1 is phosphorylated slower (rate = 0.07 ± 0.01 arb. units/min) than VEGFR2 (rate = 0.17 ± 0.02 arb. units /min). We also note that the phosphotyrosine (Y1213) in VEGFR1 is transiently stable (t1/2 = 14 ± 4 min) compared to sustained phosphorylation of Y1175 in VEGFR2 (t1/2 > 60 ± 8 min) (Fig. 2h, Supplementary Fig. 2i and Supplementary Table 4). In VEGFR1, Y1213 phosphorylation peaks approximately at ten minutes, and the phosphorylation decays by forty minutes. In contrast, in VEGFR2, the phosphorylation of Y1175 peaks approximately at five minutes, and the phosphorylation does not decay to fifty percent by one hour. We next determine the phosphorylation kinetics of the total phosphotyrosines to find out if the phosphorylation pattern of VEGFR1 and VEGFR2 is not specific to the Y1213 and Y1175 but a general property of the respective VEGFR isotype (Supplementary Fig. 1b). We observed that the decay kinetics of the total phosphotyrosines follow the same pattern as of Y1213 and Y1175 in VEGFR1 and VEGFR2, respectively. We speculate that the slow phosphorylation rate and transient stability of phosphotyrosine residue in VEGFR1 may contribute to a lower fraction of phosphorylated tyrosine residue (Supplementary Fig. 2f). We ask why the VEGFR1 phosphorylation is transient.
Deletion of ECD does not constitutively activate the VEGFR1
The ligand-independent activation of VEGFR is obstructed by electrostatic repulsion between the Ig-like domain (D4-7) in the ECD dimer interface (Supplementary Fig. 3a)65,66,67. Despite that, substituting cysteine at 482 with arginine (C482R mutation) in the D5 of VEGFR2, linked to infantile hemangioma68, constitutively activates the kinase by stabilizing a ligand-independent dimer26. A similar pathogenic cysteine to arginine substitution was reported for fibroblast growth factor receptors (FGFR)54,69. This suggests that a conserved ligand-independent activation mechanism prevails in RTKs carrying similar Ig-like ECD fold. However, no such mutation has been reported for VEGFR1. We next investigated if mutating the homologous C482 to arginine constitutively activates VEGFR1. We replaced C471R in the D5 of VEGFR1 and determined its activation (Fig. 3a, b). As expected, the VEGFR2 C482R mutant is constitutively activated and linearly phosphorylates Y1175 even without a ligand (Fig. 3c and Supplementary Fig. 3c, e). Surprisingly, the VEGFR1 C471R mutant, in the unligated state, is constitutively autoinhibited (Fig. 3d and Supplementary Fig. 3b, d). We wonder if the inability to dimerize renders the VEGFR1 C471R mutant inactive.
We next study the oligomeric states of VEGFR1 and VEGFR2 by fluorescence recovery after photobleaching (FRAP) experiment (Fig. 3e and Supplementary Fig. 4)70. The oligomeric status of VEGFR constructs was determined from the diffusion coefficient (Dconfocal) derived from the rate of fluorescence recovery at the bleached spot on the plasma membrane. We considered that the dimeric state would have slower Dconfocal relative to the monomeric state55. In our experiment, we used two chimeric constructs of VEGFR1 as a monomer (named VEGFR1-GPA-G83I) and dimer (named VEGFR1-GPA) control, where the TM helix is replaced by glycophorin-A (GPA) G83I mutant and wild-type GPA, respectively (Fig. 3a)71,72. The relative increase in the Dconfocal for the VEGFR1-GPA-G83I mutant (0.033 ± 0.013 µm2s−1) confirms a dimer-to-monomer transition on mutating TM segment in the VEGFR1-GPA chimera (0.021 ± 0.005 µm2 s−1) (Fig. 3e, Supplementary Fig. 4g, and Supplementary Table 3). We first turned to VEGFR2 to determine the Dconfocal for the wild-type and C482R mutant in the presence and absence of VEGF165, respectively (Fig. 3e, Supplementary Fig. 4e, f, and Supplementary Table 3). Overall, Dconfocal for VEGFR2 agrees with the recently published data26. The wild-type VEGFR2 (Dconfocal = 0.021 ± 0.008 µm2 s−1) tends to form a ligand-independent dimer, which explains the concentration-dependent activation of VEGFR2 in the absence of ligand (Fig. 2e). The ligand binding reorients the ECD and induces oligomerization (Dconfocal = 0.011 ± 0.004 µm2 s−1) mediated by a homotypic interaction between D4, D5, and D7 (Fig. 3e, Supplementary Figs. 3a, and 4c–f)17,23,24. Our data shows that the VEGFR2 C482R mutant forms a stable ligand-independent dimer (Dconfocal = 0.017 ± 0.006 µm2 s−1) that spontaneously activates the kinase domain (Fig. 3c, e and Supplementary Table 3)26.
We then evaluated the dimerization propensity for the VEGFR1 constructs and made the following observations (Fig. 3e, Supplementary Figs. 4, 5a–b, and Supplementary Table 3). 1) The VEGFR1 does not dimerize in the absence of a ligand (Dconfocal = 0.038 ± 0.018 µm2 s−1), and the ligand binding induces receptor dimerization (Dconfocal = 0.018 ± 0.007 µm2 s−1) (Fig. 3e, and Supplementary Fig. 4a, b, e, f). 2) Formation of the ligand-dependent dimer is independent of kinase activity as suggested by the Dconfocal (0.016 ± 0.006 µm2 s−1) of the kinase-dead mutant (D1022N) (Supplementary Fig. 5a-b). 3) The C471R mutation in the ECD of the VEGFR1 does not induce spontaneous dimerization (Dconfocal = 0.035 ± 0.014 µm2 s−1). The mutant only dimerizes upon VEGF165 binding (Dconfocal = 0.017 ± 0.007 µm2 s−1). In summary, our data indicate that VEGFR1 remains predominantly an inactive monomer in the unligated state. Perhaps the polarized electrostatic surface of the ECD (D4-D7) prevents receptor dimerization in the absence of ligand24,65,67, and mutation of C471 does not induce autophosphorylation of the receptor. Therefore, we speculate that removing ECD inhibition might spontaneously activate the VEGFR1, as observed in many RTKs73,74,75,76 and is often linked to pathological manifestations77.
To test this, we measured the autophosphorylation of Y1213 and Y1175 in the ECD-deleted (ΔECD) construct of VEGFR1 and VEGFR2, respectively (Fig. 3f, g and Supplementary Fig. 3f, g). As shown previously76, the VEGFR2 ΔECD construct was constitutively activated (Supplementary Fig. 3f) and linearly phosphorylates Y1175 in the single-cell assay (Fig. 3f). Counter-intuitively, the deletion of ECD did not activate the VEGFR1 even at the higher receptor concentration (Fig. 3g and Supplementary Fig. 3g). We speculate that the TM-JM segment connecting the ECD and kinase domain (Fig. 1) may be constitutively inhibiting the ligand-independent activation of VEGFR1.
The transmembrane domain does not drive ligand-independent VEGFR1 activation
The TM segment is a major driving force for RTK dimerization. The dynamic equilibrium between receptor dimer and monomer is rotationally coupled to the orientation of the TM segment78,79. VEGFR2 TM segment adopts two dimer structures, ligand-independent and ligand-dependent (Fig. 4d)26,80. The sequence comparison between the VEGFR1 and VEGFR2 shows that the residues at the ligand-independent dimer interface are conserved (Fig. 4d). In contrast, the residues at the ligand-dependent dimer interface are not conserved. We speculate that the presence of T763 and C764 may bias the TM structure to the ligand-dependent dimer over the ligand-independent dimer81,82. Thus, in the absence of a ligand, structural incompatibility between the active TM dimer and the unligated ECD dimer may render the VEGFR1 to remain in the monomeric state83. Ligand binding favors the TM structure towards a ligand-dependent dimer. To test that, we mutated the T761, T763, or C764 individually in VEGFR1 to the corresponding residue in VEGFR2 and measured the Y1213 phosphorylation (Supplementary Fig. 5c). We also replaced the TM segment of VEGFR1 with VEGFR2 in the full-length and ΔECD construct of VEGFR1 (Fig. 4a–c, e, and Supplementary Fig. 5c). We observed that none of the single-mutant and TM chimeric constructs could activate the VEGFR1 ligand independently. The FRAP analysis suggests that the full-length VEGFR1-TMVEGFR2 chimera has a higher propensity to form a ligand-independent dimer (Dconfocal = 0.022 ± 0.007 µm2 s−1) than the wild-type VEGFR1 (Fig. 4f, Supplementary Fig. 5e, f, and Supplementary Table 3). Suggesting the TM segment of VEGFR1 is a weak dimerization motif compared to VEGFR2 in a ligand-free state. To determine if a stronger TM dimerization motif could spontaneously phosphorylate Y1213 independent of VEGF165 stimulation, we turned to the VEGFR1-TMGPA chimera (Supplementary Fig. 4i). We observed that even the VEGFR1-TMGPA chimera could not phosphorylate the Y1213 constitutively. The inability of the VEGFR1-TMVEGFR2 or VEGFR1-TMGPA dimer to activate the kinase-domain ligand independently is counterintuitive. These data indicate that the regulatory elements downstream of the TM segment may constitutively autoinhibit VEGFR1 in the unligated state. Therefore, we replaced the JM or TM-JM segment of VEGFR1 with the VEGFR2 in the ΔECD background. Replacing the JM segment spontaneously activates the kinase and linearly phosphorylates Y1213 (Fig. 4b, c and Supplementary Fig. 5d). Together, we conclude that the JM segment is a key regulator of VEGFR1 activation.
The electrostatic latch stabilizes the inactive conformation of the JM segment
A repressor sequence present in the JM segment of VEGFR1 is known to inhibit the downstream signaling and cell migration constitutively30. The JM segment of VEGFR1 and VEGFR2 are homologous and have minor differences in the amino acid sequence (Fig. 5a). In the autoinhibited state, the kinase domain of VEGFR1 (PDB ID: 3HNG) and VEGFR2 (PDB ID: 4AGC)21 adopts a JM-in-like inactive conformation, found in the PDGFR family of kinases (Supplementary Fig. 6a)19,20,21. The JM-B segment is buried deep into the catalytic site, stabilizing the folded conformation of the activation loop and preventing the rearrangement of the N- and C-lobe to an active state. The conformation of the JM-Z region sets the direction of the rest of the JM segment inward to the kinase domain. In spite of sharing a high degree of sequence and structural homology, it is not clear how the JM segment of VEGFR1 is differentially regulated from VEGFR2.
To find answers, we revisited the structure of the VEGFR and PDGFR family of kinases. In the crystal structure of VEGFR, the JM-S segment was unresolved. We model the JM-S segment of VEGFR1 and VEGFR2 based on the inactive structure of PDGFR (PDB ID: 5K5X)84 (Supplementary Fig. 6a). The structural evaluation revealed two key aspects: First, the JM-S segment carries an overall negative charge. Second, the C-lobe of VEGFR1 has a positive charge patch, which is absent in the other VEGFR and PDGFR family (Fig. 5a–c). In our model, the positive charge residues (K1142, K1079, or R1146) in the C-lobe of the VEGFR1 kinase domain form salt bridges with the negatively charged residues (D802 and E803) in the JM-S segment (Fig. 5b). In VEGFR2, the corresponding residues in the C-lobe do not form salt bridges with the JM-S segment (Fig. 5c and Supplementary Fig. 6g). We hypothesized that the unique salt bridge between the JM-S and the C-lobe acts as an electrostatic latch that stabilizes the JM-in conformation of VEGFR1, rendering it constitutively inactive in the ligand-free state.
To test our hypothesis, we use molecular dynamics simulation to find the relative stability of the electrostatic latch in VEGFR1 and VEGFR2, respectively (Supplementary Fig. 6f, g). Our analysis of the distance between the ion pairs suggests that the electrostatic latch in VEGFR1 may be more stable than in VEGFR2 (Supplementary Fig. 6f–i). We observed that the E803 and D802 in the JM-S segment of VEGFR1 maintained electrostatic contact with the respective C-lobe residues during the simulation. We speculate that the electrostatic latch may be an integral component of the autoinhibited VEGFR1 structure and may regulate the transition between an inactive to an active conformation.
Removing JM inhibition increases the ligand-independent activation of VEGFR1
To evaluate the structure and function of the electrostatic latch, we interrogate two VEGFR1 constructs, VEGFR1-JMVEGFR2 chimera, and triple mutant (3 M) (where positively charged residues in the C-lobe K1142S, K1079Q, and R1146T are mutated to the corresponding residues in VEGFR2) (Figs. 5d, e,6a). Our structural model shows that the electrostatic latch is broken in the VEGFR1-JMVEGFR2 and 3 M construct (Fig. 5d, e and Supplementary Fig. 6h–i). We speculate that perturbing the electrostatic latch may destabilize the autoinhibitory interaction of the JM-B. Thus, replacing the VEGFR1 JM segment with VEGFR2 or the triple mutant (3 M) may restore ligand-independent activation of VEGFR1. In the single-cell assay, replacing VEGFR1 JM or TM-JM segments with VEGFR2 (Fig. 6b, c) restores the concentration-dependent autophosphorylation in the ligand-free state. The 3 M mutant partially restored the Y1213 phosphorylation, suggesting a critical role for the electrostatic latch in stabilizing the inactive JM-in structure (Fig. 6d). However, the complete restoration of the ligand-independent VEGFR1 autophosphorylation might require additional JM restraint to be removed.
To investigate why perturbing the electrostatic latch does not fully restore the ligand-independent activation of VEGFR1, we revisited the JM-in structure of VEGFR2. We observed that the conserved Y801 in the JM-B region of VEGFR2, which forms an H-bond with the critical glutamic acid residue in C-helix in the PDGFR19,85, is moved out of the catalytic site and does not interact with the C-helix (Supplementary Fig. 7a, b). The JM-B segment of VEGFR1, which is shorter by one residue than VEGFR2, is unresolved in the crystal structure (PDB ID 3HNG) (Fig. 5a). Based on structural modeling, we predict that moving the corresponding Y794 in VEGFR1 to (−)1 position may place the Y794 in the catalytic site, allowing it to interact with the glutamic acid residue in C-helix (Supplementary Fig. 7a, b). We hypothesize that if the Y794 is moved out of the catalytic site and simultaneously removing the electrostatic latch may activate the VEGFR1 ligand independently. Using the single-cell assay, we determined the ligand-independent activation of VEGFR1 in a ΔS mutant (where three consecutive serine residues, 786-788 at the JM-B, are removed) and a ΔS mutant in the 3 M background (Fig. 6a). We observed that the ΔS mutant could partially restore the Y1213 phosphorylation without ligand (Fig. 6d). To our delight, the ΔS mutant introduced in the 3 M background restored the ligand-independent VEGFR1 activation to a level comparable to the VEGFR1-JMVEGFR2 chimera (Fig. 6d and Supplementary Fig. 7c). We ask if removing the JM inhibition is enough to remodel the phosphorylation kinetics of VEGFR1.
Removing JM inhibition remodels transient phosphorylation of VEGFR1 to sustained phosphorylation
The FRAP experiment shows an increased dimeric propensity for the VEGFR1-JMVEGFR2 (Dconfocal = 0.0205 ± 0.0059 µm2 s−1) and VEGFR1-TMJMVEGFR2 (Dconfocal = 0.0202 ± 0.007 µm2 s−1) chimera (Fig. 6e, Supplementary Fig. 7d-e, and Supplementary Table 3). Suggesting the JM-in conformation of VEGFR1 is incompatible with the ligand-independent TM dimer. Stimulating VEGFR1-JMVEGFR2 (Dconfocal = 0.0127 ± 0.0047 µm2 s−1) and VEGFR1-TMJMVEGFR2 (Dconfocal = 0.0131 ± 0.005 µm2 s−1) with VEGF165 induces the receptor oligomerization and also increases the relative fraction of Y1213 phosphorylation compared to wild-type VEGFR1 (Fig. 6e, f and Supplementary Table 3). Removing the JM inhibition in the VEGFR1-TMJMVEGFR2 chimera increases the rate of Y1213 phosphorylation (0.16 ± 0.03 arb. units/min). It also remodels the phosphorylation half-life from transient to sustained (t1/2 = 48 ± 13.8 min) (Fig. 6g, h, Supplementary Fig. 7f–h, and Supplementary Table 4). We conclude that the multiple interactions between the JM segment and the kinase core are the dominant force in autoinhibiting the VEGFR1 constitutively in the unligated state.
Discussion
To summarize, we presented a molecular mechanism explaining how VEGFR1 and VEGFR2 are differentially autophosphorylated in the unligated state and upon ligand stimulation (Fig. 7). An equilibrium shift model of multiple species including a monomeric receptor and various dimers, previously proposed for EGFR, explains the ligand-independent activation of VEGFR258,59,86. The ECD of EGFR in the unligated state forms a head-to-head dimer that keeps the TM segment separated, preventing the two kinase domains from adopting an asymmetric active dimer87,88,89. Consistent with the previous finding, our results show that the TM segment in VEGFR2 dimerizes in the absence of a ligand, allowing the two kinase domains to transphosphorylate at a critical receptor concentration26,54. In contrast, we found that VEGFR1 remains predominantly as a monomer in the absence of a ligand, and replacing the TM segment of VEGFR1 with VEGFR2 allows the receptor to dimerize independently of the ligand (Fig. 4f). However, forced dimerization of the VEGFR1 failed to activate the kinase domain (Fig. 4e and Supplementary Fig. 4i).
Unlike other RTKs59,76, removal of the ECD does not spontaneously activate the VEGFR1 (Fig. 4c). Replacing the JM segment of the VEGFR1 with that of VEGFR2 activates the kinase domain ligand independently (Fig. 6c). Our data supports a central role for the JM segment in balancing the equilibrium between an inactive and active state of VEGFR1. In the unligated state, the electrostatic latch and H-bond interaction steaming from Y794 in JM-B may shift the equilibrium towards the inactive state, thus making VEGFR1 an inefficient kinase. The JM-in conformation of VEGFR1 suppresses the spontaneous kinase activation (Figs. 5 and 6)26. We therefore proposed that slow kinase activation along with the cellular protein tyrosine (PTP) activity90,91 maintains the basal activity of VEGFR1 constitutively inhibited (Fig. 7). The role of the JM segment of VEGFR is in contrast to the EGFR, where the JM segment of the receiver kinase stabilizes the active asymmetric dimer by interacting with the C-lobe of the activator and couples to the TM dimer19,59,89,92.
The asymmetric dimer of the EGFR kinase domain generates asymmetry in the phosphorylation of tyrosine residues in the C-terminal tail, with a clear preference for activator93. A possible asymmetry in the phosphorylation of tyrosine residues in the C-terminal tail may occur in an asymmetric heterodimer of VEGFR1 and VEGFR225. Our data suggests that the heterodimer of two kinetically different kinases may generate distinctly different tyrosine phosphorylation patterns compared to the homodimer, hence generating an altered signaling output25,94,95.
The ligand bias by the ECD dimer generates differential signaling output in EGFR96,97 and is a crucial determinant for deciding the cell fate98. A transient versus sustained Erk activation generated by stimulating EGFR with two different ligands, EGF or NGF, switches the signaling outcomes to differentiation from the proliferation99,100. Recent studies suggest that the subtle structural difference in the ligand-induced ECD dimer may determine the differential signaling output in EGFR96,97. Here, we studied how a common ligand induces different outputs in two VEGF receptor isomers. We proposed that the slow removal of JM inhibition in VEGFR1 may generate transient phosphorylation of tyrosine residues (Fig. 7). We speculate that the bias signaling by the two VEGF receptor isomers may account for the differential response of VEGFR1 in promoting macrophage proliferation101, whereas VEGFR2 supports differentiation of endothelial cell (EC) progenitors102.
Methods
DNA constructs
The cDNA encoding wild-type human VEGFR1 (Residue number 1–1338; pDONR223-FLT1 was a gift from William Hahn & David Root; Addgene plasmid # 23912) was subcloned into pcDNA-mCherry vector between the HindIII and Kpn1 restriction sites. The plasmid encoding human wild-type VEGFR2 tagged with mCherry (cloned in pBE-vector; 108854) was a gift from Kalina Hristova (Johns Hopkins University, Baltimore, MD). All point mutations and deletion constructs were prepared by PCR-based mutagenesis strategy103. For the VEGFR1 ECD deletion constructs, Haemagglutinin (HA) tag was incorporated at the 5´ ends. The list of primers used is tabulated in Table S3. In brief, PCR-amplified product at the expected molecular weight was extracted and purified from 0.6% agarose gel and incubated with Dpn1 at 37 °C for 3 h. The Dpn1 digested PCR product was further purified from the reaction mixture and incubated with Polynucleotide Kinase 4 (PNK4) and Ligase in Tris buffer containing 10 mM MgCl2 and 2 mM ATP. Finally, the ligated mixture was transformed into the E.coli Top10 cell. Chimeric constructs were prepared using the Gibson assembly approach104. All the chimeric constructs were cloned into BamH1 and Xho1 restriction sites in a pcDNA-mCherry vector (mCherry2-N1 was a gift from Michael Davidson, Addgene Plasmid # 54517)
Cell culture, immunofluorescence
Chinese Hamster Ovary (CHO) (obtained from National Center for Cell Science- India), African green monkey kidney fibroblast-like cell line (COS-7) (obtained from National Center for Cell Science- India), and murine macrophage (RAW264.7) cell lines (obtained from ATCC, Catalog No. TIB-71) were cultured in DMEM supplemented with 10% FBS, 50 μg/ml penicillin, and streptomycin at 37 °C with 5% CO2. In the single-cell experiments and immunoblotting, VEGFR1 or VEGFR2 was transiently transfected in respective cell lines. In brief, the COS-7 or the CHO cells were grown on coverslips (10 mm, #1 thickness) up to a confluency of 80% and transfected with Lipofectamine 2000 in serum-free media (Opti-MEM). After 6 hours of incubation, cells were supplemented with complete media containing 10% FBS and allowed to grow for another 12 hours. Alternatively, RAW264.7 cells were nucleofected using the D032 program in the Lonza nucleofector. 4 × 107 cells were incubated with 4 μg of DNA in Amaxa buffer for 5 min at room temperature. Nucleofected cells were resuspended in DMEM supplemented with 15% FBS. 18 h post-transfection, cells were serum-starved by growing in serum-free opti-MEM for another 8 h. To activate VEGFR1 or VEGFR2, transfected cells were treated with 100 nM VEGF165 for 5 mins, then immediately washed with ice-cold PBS two times and fixed by treating with 4% paraformaldehyde (PFA) in PBS for 30 min at room temperature. After fixation, cells were washed five times with 1× PBS and permeabilized with 0.2% PBST (1X PBS and 0.2% Triton-X 100) for 5 min at room temperature. Cells were blocked with 1% BSA for 1 h, followed by staining with anti-phosphotyrosine antibodies for VEGFR2 or VEGFR1 overnight at 4 °C (dilution 1:200 for anti phospho-VEGFR1 and anti phospho-VEGFR2 antibody). After overnight incubation, coverslips were washed three times with PBST, followed by incubation with a secondary antibody conjugated with FITC for 2 h at room temperature. Finally, coverslips were mounted with prolonged gold. in between each step, the coverslips were washed three times with PBST and finally with PBS before mounting. All antibodies used, along with the catalog numbers and dilutions, are mentioned in Supplementary Table 1.
Microscopy and image analysis
Transiently transfected COS-7 cells were imaged by an Olympus IX81 epifluorescence microscope. For CHO and RAW264.7 cell line, all images were acquired with the Leica SP8 confocal platform using an oil immersion HC PL APO pinhole was set at 1 AU. Image scanning was done in bidirectional mode at 500 Hz. Expression of mCherry fused VEGFR constructs and the phosphorylation levels was imaged using 40 mW/552 nm and 20 mW/488 nm solid-state lasers, respectively. For mCherry, λex and λem were set to 552 nm and 576–651 nm, respectively. For FITC, λex and λem were set to 488 nm and 505-531 nm, respectively. Laser power was set at 0.08 mW and 0.04 mW for mCherry and FITC channels, respectively. The background fluorescence was corrected by subtracting the respective fluorescence measured for untransfected cells. Image analysis was done by selecting ROIs drawn manually using ImageJ freehand tool105, as explained previously59. Two sets of ROI were selected; ROI1 was drawn to calculate total cell intensity, and ROI2 was drawn at the peripheral region of the cell to calculate cytoplasmic intensity (Supplementary Fig. 1d). The fluorescence intensity of mCherry or FITC at the membrane (Im) was calculated by:
Im = IROI1-IROI2, where IROI1 and IROI2 are fluorescence of mCherry or FITC at the respective ROI.
To analyze the tyrosine phosphorylation as a function of receptor expression level, the intensity of the mCherry channel was binned in the range of 2500 arb. units, and the average value of mCherry intensity was plotted against the corresponding average value of FITC intensity. For measuring the normalized phosphorylation level, the intensity from the FITC channel was normalized by the corresponding mCherry intensity for the respective bin, the cells were binned at an intensity range of 10000 arb. units.
Immunoblotting
Activation of VEGFR1 or VEGFR2 was performed as mentioned in the above section. Briefly, CHO cells were cultured in DMEM and transfected with VEGFR1 and VEGFR2 constructs. 18 h post transfection, cells were serum-starved by growing in serum-free opti-MEM for another 8 h and activated with 100 nM VEGF165 for 5 mins. After activation, cells were washed with ice-cold PBS before being incubated in RIPA buffer (10 mM Tris-Cl, pH 8.0, 140 mM NaCl, 1 mM EDTA, 0.1% SDS, and 1% Triton X-100) containing protease and phosphatase inhibitors (2 mM Benzamidine, 1 mM PMSF), and then sonicated on ice. Protein samples were prepared by heating them with 5X loading buffer, resolving them on a 6% SDS-PAGE, and blotted onto a PVDF membrane. The blot was blocked with 5% skimmed milk in 1× TBS and 0.1% TWEEN-20 for 1 h at RT and then incubated with primary antibody at 4 °C overnight (dilution 1:1000 for anti VEGFR2, anti Phospho-VEGF Receptor 2, anti VEGFR1 and anti Phospho-VEGFR1 Antibody; 1:2000 for anti HA antibody). The unbound primary antibody was removed by washing the blot three times with 1× TBST (0.1% tween-20), followed by incubation with a secondary antibody diluted in 3% skimmed milk (dilution 1: 2000 for rabbit and goat HRP secondary antibody; 1:3000 for mouse HRP secondary antibody). Blot was washed three times with 1× TBST and two times with 1XTBS wash before developing with the ClarityTM Western ECL substrate kit (Bio-Rad). The images were acquired using the Bio-Rad Chemidoc system. The densitometric analysis was performed using ImageJ. All antibodies used, along with the catalog numbers and dilutions, are mentioned in Supplementary Table 1.
Cell-based tyrosine phosphorylation kinetic measurements
Chinese Hamster Ovary (CHO) cell line was cultured in DMEM and transfected with VEGFR constructs. 18 h post transfection, cells were serum-starved in Opti-MEM for another 8 h. The VEGFR tyrosine phosphorylation was stimulated by treating with VEGF165. Samples were collected at the indicated time points and analyzed by immunoblotting (dilution 1:1000 for Anti phosphotyrosine Antibody, anti VEGFR2, anti Phospho-VEGF Receptor 2, anti VEGFR1 and anti Phospho-VEGFR1 Antibody; 1:2000 for anti HA antibody). Densitometry analysis of immunoblots was performed using ImageJ software. All antibodies used, along with the catalog numbers and dilutions, are mentioned in Supplementary Table 1.
FRAP experiment
CHO cells were seeded in a 35 mm glass-bottomed Petri dish at a density of 0.5 × 106 cells per plate. 24 h post seeding, cells were transfected with VEGFR1 or VEGFR2 constructs. The cells were grown for 18 h post-transfection and were serum-starved for 6 h before ligand stimulation. Before activating with VEGF165, cells were washed with pre-warmed PBS(1×) and activated with 100 nM VEGF165 in Opti-MEM.
FRAP experiments were conducted using a Leica SP8 confocal microscope at room temperature. Imaging was performed using an HC PL APO CS2 63X/ 1.4 NA objective at 5× digital zoom. The confocal pinhole was set at 1 Airy unit. Image scanning was done in bidirectional mode at 500 Hz. The excitation and emission filters for the mCherry channel were set to 552 nm and 576–651 nm, respectively. A circular region of radius 1.2 μm (bleached spot) at the cell edge was bleached using a 40 mW solid-state laser (100% intensity) for 500 ms, with a pixel dwell time of 1.15 μs. The laser powers for pre-bleach and post-bleach imaging were set at 1% (0.4 mW). 10 pre-bleached and 120 post-bleached frames were recorded in 512 × 512 pixel format at 2frame/s.
FRAP analysis
The measured raw fluorescence intensity was corrected for background fluorescence and photofading effects. At first, the fluorescence intensity at each frame was measured for the ROI defined by the bleached spot [FROI(t)]. To account for background fluorescence [Fbg(t)], the fluorescence was measured in a cell-free region of the image in the respective frame. To adjust for observational photofading [Ftotal(t)], the fluorescence from the total cell, except the bleached spot, was measured. The raw FRAP data [FROI(t)], were corrected for background and photofading at each frame by:
Finally, the corrected FRAP (\({F}_{{corrected}}\,\left(t\right)\)) was normalized by corrected prebleach intensity (\({F^{\prime}}_{{corrected}}\)):
(Where \({F}_{{ROI}}\left({t}_{0}\right)\), \({F}_{{Total}}\left({t}_{0}\right)\), and \({F}_{{bg}}\left({t}_{0}\right)\) indicate the initial fluorescence intensity from the bleached spot, total cell, and background, respectively)
The effective bleached spot radius (re) was determined (Supplementary Fig. 4h)70. The half-life of fluorescence recovery (t1/2) was determined by fitting corrected fluorescence intensities at indicated time points to the first-order kinetic equation:
Finally, the diffusion coefficient (Dconfocal) was derived from the modified Soumpasis equation106:
where re is the effective bleached spot radius and the coefficient 0.25 was numerically determined70
We validated our FRAP experimental setup by comparing the diffusion coefficients (Dconfocal) of ligand-free EGFR-mCherry and VEGFR2-mCherry with previously reported diffusion coefficient measurements using single-molecule tracking. We observed that the diffusion coefficient of EGFR (0.033 μm2 s−1) and VEGFR2 (0.020 μm2 s−1) (Table S3) determined from Confocal FRAP is in agreement with the average diffusion coefficient value reported for the heterogeneous population of EGFR (0.036 μm2 s−1)96 and VEGFR2 (0.032 μm2 s−1)107, respectively.
Homology modelling
The models for the autoinhibited conformation of VEGFR1 (781-1158), VEGFR2 (787-1162), VEGFR1-JMVEGFR2 chimera, and VEGFR1 triple mutated (K1079Q, K1142S, and R1146D) were built using iTasser server (https://zhanggroup.org/I-TASSER/)108. In all our models, the kinase domain adopts an autoinhibited conformation of C-helix-in and DFG–out. The kinase domain of VEGFR-1 (825-1158) was modelled based on PDB ID 3HNG, and the juxtamembrane segment (781-824) was modelled based on PDGFR structure (PDB ID 5K5X) (Supplementary Fig. 6a). The linker region between the N-lobe and C-lobe (925-991) could not be modelled due to the lack of any homology structure and remained unstructured in the model. The final model, with the highest confidence score as well as a low root-mean-square-deviation (see Table S5, RMSD of 0.589 Å w.r.t 3HNG) with respect to VEGFR1 structure (PDB ID: 3HNG), was selected and the inhibitor N-(4-chlorophenyl)−2-[(pyridin-4-ylmethyl)amino]benzamide was docked at the ATP binding pocket based on 3HNG structure, using PyMOL (DeLano, W. L., 2009). The structure was further energy minimized and equilibrated before analyzing by MD simulation. A preliminary MD simulation of 100 ns was performed to determine the structural stability of the selected model. The lowest energy structure selected remained stable in the preliminary simulation and was thus used for further studies. The models for the VEGFR1-JMVEGFR2 chimera and the triple mutant construct of VEGFR1 were similarly constructed, energy minimized, and evaluated for further structural studies. Similarly, VEGFR2 (PDB ID 4AGC) was used to model the kinase domain (832-1162), and the juxtamembrane segment (787-831) was modelled based on PDGFR structure (PDB ID 5K5X). The model was selected as explained previously, and Axitinib was docked at the ATP binding pocket based on 4AGC structure, using PyMOL.
Molecular dynamics simulations
The MD simulations were run using GROMACS 2019.6109,110. System preparation was done using the CHARMM-GUI server (www.charmm-gui.org/) with TIP3p water molecules111 and 0.15 M NaCl. The model required 1000 steps approximately for energy minimization. The solvated system was equilibrated for 125 ps at a 1 fs time step, using H-bonds as constraints by implementing Linear Constraint Solver for Molecular Simulations (LINCS) algorithm112. At least three independent simulations of 2 µs each for VEGFR1 (wildtype), 500 ns each for VEGFR2 (wildtype), VEGFR1-JMVEGFR2 chimera, and VEGFR1 triple mutated were performed using the CHARMM36m force field113. To comparatively study the orientation of the juxtamembrane segment with respect to the kinase domain in the crystal structure of cFMS (PDB ID 2OGV)85 and FLT3 (PDB ID 1RJB)20, we performed three independent simulation of 200 ns each. The simulations were performed under constant pressure (1 bar)114 and temperature (300 K) (NPT)115 and a time step of 2 fs. Potential-shift-Verlet was used for electrostatic and van der Waals interactions with a 12-Å cutoff. The trajectories were visually analyzed using VMD116, and the structures were visualized in PyMol. RMSDs, RMSFs, and distances were measured using the tools provided in GROMACS.
Statistics and reproducibility
Statistical analyses were performed with GraphPad Prism 9 using an unpaired two-tailed student’s t-test. No statistical method was used to predetermine the sample size. No data were excluded from the analyses. All data were presented as mean ± SD from at least three biologically independent experiments. P < 0.05 was considered as statistically significant.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The kinase domain structures used for the analysis of the conserved signature motifs in VEGFR1 and VEGFR2 can be found in the PDB: 3VHK (Crystal structure of the VEGFR2 kinase domain in complex with a back pocket binder), 1VR2 (Human vascular endothelial growth factor receptor 2 (KDR) Kinase domain) and 3HNG (Crystal structure of VEGFR1 in complex with N-(4-Chlorophenyl)−2-(pyridin-4-ylmethyl)amino)benzamide). The structure used for comparative structural analysis of VEGFR1 and VEGFR2 extracellular domain can be found in the PDB: 5T89 (Structure of VEGF-A in complex with VEGFR-1 domains D1-6). The JM segment of the VGEFR2 kinase domain (PDB: 4ASD (Crystal Structure of VEGFR2 Juxtamembrane and Kinase Domains) in Complex with SORAFENIB (BAY 43-9006) was modeled based on PDGFR crystal structure (PDB: 5K5X (Crystal structure of human PDGFRA)). The JM boundary of VEGFR2 was defined using FMS-like tyrosine kinase 3 structure (PDB: 1RJB (Crystal Structure of FLT3)). The kinase domain structures used to analyze the inhibitory interaction between the JM segment and activation loop in cFMS and VEGFR2 can be found in the PDB: 2OGV (Crystal Structure of the Autoinhibited Human c-Fms Kinase Domain) and 4AGC (Crystal structure of VEGFR2 juxtamembrane and kinase domains in complex with Axitinib (AG-013736) (N-Methyl-2-(3-((E)−2-pyridin-2-yl- vinyl)−1H-indazol-6-ylsulfanyl)-benzamide). All the relevant data are contained within this article and in the supporting information. Source data are available in the Source Data file and as a Figshare deposition (https://doi.org/10.6084/m9.figshare.24241207). Uncropped blots are available in the Source Data file. Source data are provided with this paper.
References
Ferrara, N., Gerber, H. P. & LeCouter, J. The biology of VEGF and its receptors. Nat. Med. 9, 669–676 (2003).
Olsson, A. K., Dimberg, A., Kreuger, J. & Claesson-Welsh, L. VEGF receptor signalling - in control of vascular function. Nat. Rev. Mol. Cell Biol. 7, 359–371 (2006).
de Vries, C. et al. The fms-like tyrosine kinase, a receptor for vascular endothelial growth factor. Science 255, 989–991 (1992).
Maru, Y., Yamaguchi, S. & Shibuya, M. Flt-1, a receptor for vascular endothelial growth factor, has transforming and morphogenic potentials. Oncogene 16, 2585–2595 (1998).
Fischer, C., Mazzone, M., Jonckx, B. & Carmeliet, P. FLT1 and its ligands VEGFB and PlGF: drug targets for anti-angiogenic therapy? Nat. Rev. Cancer 8, 942–956 (2008).
Simons, M., Gordon, E. & Claesson-Welsh, L. Mechanisms and regulation of endothelial VEGF receptor signalling. Nat. Rev. Mol. Cell Biol. 17, 611–625 (2016).
Cao, Y. Positive and negative modulation of angiogenesis by VEGFR1 ligands. Sci. Signal 2, re1 (2009).
Hiratsuka, S., Minowa, O., Kuno, J., Noda, T. & Shibuya, M. Flt-1 lacking the tyrosine kinase domain is sufficient for normal development and angiogenesis in mice. Proc. Natl Acad. Sci. USA 95, 9349–9354 (1998).
Sawano, A., Takahashi, T., Yamaguchi, S., Aonuma, M. & Shibuya, M. Flt-1 but not KDR/Flk-1 tyrosine kinase is a receptor for placenta growth factor, which is related to vascular endothelial growth factor. Cell Growth Differ. 7, 213–221 (1996).
Shinkai, A. et al. Mapping of the sites involved in ligand association and dissociation at the extracellular domain of the kinase insert domain-containing receptor for vascular endothelial growth factor. J. Biol. Chem. 273, 31283–31288 (1998).
Waltenberger, J., Claesson-Welsh, L., Siegbahn, A., Shibuya, M. & Heldin, C. H. Different signal transduction properties of KDR and Flt1, two receptors for vascular endothelial growth factor. J. Biol. Chem. 269, 26988–26995 (1994).
Seetharam, L. et al. A unique signal transduction from FLT tyrosine kinase, a receptor for vascular endothelial growth factor VEGF. Oncogene 10, 135–147 (1995).
Lemmon, M. A. & Schlessinger, J. Cell signaling by receptor tyrosine kinases. Cell 141, 1117–1134 (2010).
Lemmon, M. A. et al. Two EGF molecules contribute additively to stabilization of the EGFR dimer. EMBO J. 16, 281–294 (1997).
Yarden, Y. & Schlessinger, J. Epidermal growth factor induces rapid, reversible aggregation of the purified epidermal growth factor receptor. Biochemistry 26, 1443–1451 (1987).
Huang, Y. et al. Molecular basis for multimerization in the activation of the epidermal growth factor receptor. Elife 5, e14107 (2016).
Yang, Y., Xie, P., Opatowsky, Y. & Schlessinger, J. Direct contacts between extracellular membrane-proximal domains are required for VEGF receptor activation and cell signaling. Proc. Natl Acad. Sci. USA 107, 1906–1911 (2010).
Zhala, A. et al. Fluorescent Resonance Energy Transfer Imaging of VEGFR Dimerization. Anticancer Res. 34, 2123 (2014).
Hubbard, S. R. Juxtamembrane autoinhibition in receptor tyrosine kinases. Nat. Rev. Mol. Cell Biol. 5, 464–471 (2004).
Griffith, J. et al. The structural basis for autoinhibition of FLT3 by the juxtamembrane domain. Mol. Cell 13, 169–178 (2004).
McTigue, M. et al. Molecular conformations, interactions, and properties associated with drug efficiency and clinical performance among VEGFR TK inhibitors. Proc. Natl Acad. Sci. USA 109, 18281–18289 (2012).
Solowiej, J. et al. Characterizing the effects of the juxtamembrane domain on vascular endothelial growth factor receptor-2 enzymatic activity, autophosphorylation, and inhibition by axitinib. Biochemistry 48, 7019–7031 (2009).
Ruch, C., Skiniotis, G., Steinmetz, M. O., Walz, T. & Ballmer-Hofer, K. Structure of a VEGF-VEGF receptor complex determined by electron microscopy. Nat. Struct. Mol. Biol. 14, 249–250 (2007).
Markovic-Mueller, S. et al. Structure of the Full-length VEGFR-1 extracellular domain in complex with VEGF-A. Structure 25, 341–352 (2017).
Opatowsky, Y. et al. Structure, domain organization, and different conformational states of stem cell factor-induced intact KIT dimers. Proc. Natl Acad. Sci. USA 111, 1772–1777 (2014).
Sarabipour, S., Ballmer-Hofer, K. & Hristova, K. VEGFR-2 conformational switch in response to ligand binding. Elife 5, e13876 (2016).
Bae, J. H. & Schlessinger, J. Asymmetric tyrosine kinase arrangements in activation or autophosphorylation of receptor tyrosine kinases. Mol. Cells 29, 443–448 (2010).
Jura, N. et al. Catalytic control in the EGF receptor and its connection to general kinase regulatory mechanisms. Mol. Cell 42, 9–22 (2011).
Kornev, A. P., Taylor, S. S. & Ten Eyck, L. F. A helix scaffold for the assembly of active protein kinases. Proc. Natl Acad. Sci. USA 105, 14377–14382 (2008).
Gille, H. et al. A repressor sequence in the juxtamembrane domain of Flt-1 (VEGFR-1) constitutively inhibits vascular endothelial growth factor-dependent phosphatidylinositol 3’-kinase activation and endothelial cell migration. EMBO J. 19, 4064–4073 (2000).
Meyer, R. D., Mohammadi, M. & Rahimi, N. A single amino acid substitution in the activation loop defines the decoy characteristic of VEGFR-1/FLT-1. J. Biol. Chem. 281, 867–875 (2006).
Shin, J. Y. et al. CD4+VEGFR1(HIGH) T cell as a novel Treg subset regulates inflammatory bowel disease in lymphopenic mice. Cell Mol. Immunol. 12, 592–603 (2015).
Tchaikovski, V., Fellbrich, G. & Waltenberger, J. The molecular basis of VEGFR-1 signal transduction pathways in primary human monocytes. Arterioscler Thromb. Vasc. Biol. 28, 322–328 (2008).
Eilken, H. M. et al. Pericytes regulate VEGF-induced endothelial sprouting through VEGFR1. Nat. Commun. 8, 1574 (2017).
Niida, S. et al. VEGF receptor 1 signaling is essential for osteoclast development and bone marrow formation in colony-stimulating factor 1-deficient mice. Proc. Natl Acad. Sci. USA 102, 14016–14021 (2005).
Murakami, M. et al. Signaling of vascular endothelial growth factor receptor-1 tyrosine kinase promotes rheumatoid arthritis through activation of monocytes/macrophages. Blood 108, 1849–1856 (2006).
Benhadou, F. et al. Epidermal autonomous VEGFA/Flt1/Nrp1 functions mediate psoriasis-like disease. Sci. Adv. 6, eaax5849 (2020).
Selvaraj, D. et al. A functional role for VEGFR1 expressed in peripheral sensory neurons in cancer pain. Cancer Cell 27, 780–796 (2015).
Abdelrahim, M. et al. Regulation of vascular endothelial growth factor receptor-1 expression by specificity proteins 1, 3, and 4 in pancreatic cancer cells. Cancer Res. 67, 3286–3294 (2007).
Cao, R. et al. VEGFR1-mediated pericyte ablation links VEGF and PlGF to cancer-associated retinopathy. Proc. Natl Acad. Sci. USA 107, 856–861 (2010).
Qian, B. Z. et al. FLT1 signaling in metastasis-associated macrophages activates an inflammatory signature that promotes breast cancer metastasis. J. Exp. Med. 212, 1433–1448 (2015).
Muramatsu, M., Yamamoto, S., Osawa, T. & Shibuya, M. Vascular endothelial growth factor receptor-1 signaling promotes mobilization of macrophage lineage cells from bone marrow and stimulates solid tumor growth. Cancer Res. 70, 8211–8221 (2010).
Hirsch, F. R., Varella-Garcia, M. & Cappuzzo, F. Predictive value of EGFR and HER2 overexpression in advanced non-small-cell lung cancer. Oncogene 28, S32–37 (2009).
Allerstorfer, S. et al. FGF5 as an oncogenic factor in human glioblastoma multiforme: autocrine and paracrine activities. Oncogene 27, 4180–4190 (2008).
Palmieri, D. et al. Her-2 overexpression increases the metastatic outgrowth of breast cancer cells in the brain. Cancer Res. 67, 4190–4198 (2007).
Maris, C. et al. IGF-IR: a new prognostic biomarker for human glioblastoma. Br. J. Cancer 113, 729–737 (2015).
Xu, L. et al. Epidermal growth factor receptor regulates MET levels and invasiveness through hypoxia-inducible factor-1alpha in non-small cell lung cancer cells. Oncogene 29, 2616–2627 (2010).
Lopez-Gines, C. et al. New pattern of EGFR amplification in glioblastoma and the relationship of gene copy number with gene expression profile. Mod. Pathol. 23, 856–865 (2010).
Morrison, R. S. et al. Fibroblast growth factor receptor gene expression and immunoreactivity are elevated in human glioblastoma multiforme. Cancer Res. 54, 2794–2799 (1994).
Du, Z. & Lovly, C. M. Mechanisms of receptor tyrosine kinase activation in cancer. Mol. Cancer 17, 58 (2018).
Carraway, K. L. 3rd & Sweeney, C. EGF receptor activation by heterologous mechanisms. Cancer Cell 1, 405–406 (2002).
Hubbard, S. R., Mohammadi, M. & Schlessinger, J. Autoregulatory mechanisms in protein-tyrosine kinases. J. Biol. Chem. 273, 11987–11990 (1998).
Saraon, P. et al. Receptor tyrosine kinases and cancer: oncogenic mechanisms and therapeutic approaches. Oncogene 40, 4079–4093 (2021).
Sarabipour, S. & Hristova, K. Mechanism of FGF receptor dimerization and activation. Nat. Commun. 7, 10262 (2016).
Chung, I. et al. Spatial control of EGF receptor activation by reversible dimerization on living cells. Nature 464, 783–787 (2010).
Imoukhuede, P. I. & Popel, A. S. Quantification and cell-to-cell variation of vascular endothelial growth factor receptors. Exp. Cell Res. 317, 955–965 (2011).
Imoukhuede, P. I. & Popel, A. S. Expression of VEGF receptors on endothelial cells in mouse skeletal muscle. PLoS One 7, e44791 (2012).
Schlessinger, J. Cell signaling by receptor tyrosine kinases. Cell 103, 211–225 (2000).
Endres, N. F. et al. Conformational coupling across the plasma membrane in activation of the EGF receptor. Cell 152, 543–556 (2013).
Basagiannis, D. & Christoforidis, S. Constitutive endocytosis of VEGFR2 protects the receptor against shedding. J. Biol. Chem. 291, 16892–16903 (2016).
Ash, D. et al. The P-type ATPase transporter ATP7A promotes angiogenesis by limiting autophagic degradation of VEGFR2. Nat. Commun. 12, 3091 (2021).
Ito, N., Wernstedt, C., Engstrom, U. & Claesson-Welsh, L. Identification of vascular endothelial growth factor receptor-1 tyrosine phosphorylation sites and binding of SH2 domain-containing molecules. J. Biol. Chem. 273, 23410–23418 (1998).
Takahashi, T., Yamaguchi, S., Chida, K. & Shibuya, M. A single autophosphorylation site on KDR/Flk-1 is essential for VEGF-A-dependent activation of PLC-gamma and DNA synthesis in vascular endothelial cells. EMBO J. 20, 2768–2778 (2001).
Sakurai, Y., Ohgimoto, K., Kataoka, Y., Yoshida, N. & Shibuya, M. Essential role of Flk-1 (VEGF receptor 2) tyrosine residue 1173 in vasculogenesis in mice. Proc. Natl Acad. Sci. USA 102, 1076–1081 (2005).
Tao, Q., Backer, M. V., Backer, J. M. & Terman, B. I. Kinase insert domain receptor (KDR) extracellular immunoglobulin-like domains 4-7 contain structural features that block receptor dimerization and vascular endothelial growth factor-induced signaling. J. Biol. Chem. 276, 21916–21923 (2001).
Stuttfeld, E. & Ballmer-Hofer, K. Structure and function of VEGF receptors. IUBMB Life 61, 915–922 (2009).
Yuzawa, S. et al. Structural basis for activation of the receptor tyrosine kinase KIT by stem cell factor. Cell 130, 323–334 (2007).
Jinnin, M. et al. Suppressed NFAT-dependent VEGFR1 expression and constitutive VEGFR2 signaling in infantile hemangioma. Nat. Med 14, 1236–1246 (2008).
Sarabipour, S. & Hristova, K. Pathogenic cysteine removal mutations in FGFR extracellular domains stabilize receptor dimers and perturb the TM dimer structure. J. Mol. Biol. 428, 3903–3910 (2016).
Kang, M., Day, C. A., Kenworthy, A. K. & DiBenedetto, E. Simplified equation to extract diffusion coefficients from confocal FRAP data. Traffic 13, 1589–1600 (2012).
Lemmon, M. A., Flanagan, J. M., Treutlein, H. R., Zhang, J. & Engelman, D. M. Sequence specificity in the dimerization of transmembrane. alpha.-helixes. Biochemistry 31, 12719–12725 (1992).
Lemmon, M. A. et al. Glycophorin A dimerization is driven by specific interactions between transmembrane alpha-helices. J. Biol. Chem. 267, 7683–7689 (1992).
Arevalo, J. C. et al. A novel mutation within the extracellular domain of TrkA causes constitutive receptor activation. Oncogene 20, 1229–1234 (2001).
Qiu, F. H. et al. Primary structure of c-kit: relationship with the CSF-1/PDGF receptor kinase family–oncogenic activation of v-kit involves deletion of extracellular domain and C terminus. EMBO J. 7, 1003–1011 (1988).
Uren, A., Yu, J. C., Karcaaltincaba, M., Pierce, J. H. & Heidaran, M. A. Oncogenic activation of the alphaPDGFR defines a domain that negatively regulates receptor dimerization. Oncogene 14, 157–162 (1997).
Gonzalez-Magaldi, M., McCabe, J. M., Cartwright, H. N., Sun, N. & Leahy, D. J. Conserved roles for receptor tyrosine kinase extracellular regions in regulating receptor and pathway activity. Biochem J. 477, 4207–4220 (2020).
Gan, H. K., Cvrljevic, A. N. & Johns, T. G. The epidermal growth factor receptor variant III (EGFRvIII): where wild things are altered. FEBS J. 280, 5350–5370 (2013).
Bell, C. A. et al. Rotational coupling of the transmembrane and kinase domains of the Neu receptor tyrosine kinase. Mol. Biol. Cell 11, 3589–3599 (2000).
Dosch, D. D. & Ballmer-Hofer, K. Transmembrane domain-mediated orientation of receptor monomers in active VEGFR-2 dimers. FASEB J. 24, 32–38 (2010).
Manni, S. et al. Structural and functional characterization of alternative transmembrane domain conformations in VEGF receptor 2 activation. Structure 22, 1077–1089 (2014).
Dawson, J. P., Weinger, J. S. & Engelman, D. M. Motifs of serine and threonine can drive association of transmembrane helices. J. Mol. Biol. 316, 799–805 (2002).
Moriki, T., Maruyama, H. & Maruyama, I. N. Activation of preformed EGF receptor dimers by ligand-induced rotation of the transmembrane domain. J. Mol. Biol. 311, 1011–1026 (2001).
Krimmer, S. G. et al. Cryo-EM analyses of KIT and oncogenic mutants reveal structural oncogenic plasticity and a target for therapeutic intervention. Proc. Natl Acad. Sci. USA 120, e2300054120 (2023).
Liang, L., Yan, X. E., Yin, Y. & Yun, C. H. Structural and biochemical studies of the PDGFRA kinase domain. Biochem. Biophys. Res. Commun. 477, 667–672 (2016).
Walter, M. et al. The 2.7 A crystal structure of the autoinhibited human c-Fms kinase domain. J. Mol. Biol. 367, 839–847 (2007).
Schlessinger, J. Receptor tyrosine kinases: legacy of the first two decades. Cold Spring Harb. Perspect. Biol. 6, a008912 (2014).
Zanetti-Domingues, L. C. et al. The architecture of EGFR’s basal complexes reveals autoinhibition mechanisms in dimers and oligomers. Nat. Commun. 9, 4325 (2018).
Zhang, X., Gureasko, J., Shen, K., Cole, P. A. & Kuriyan, J. An allosteric mechanism for activation of the kinase domain of epidermal growth factor receptor. Cell 125, 1137–1149 (2006).
Jura, N. et al. Mechanism for activation of the EGF receptor catalytic domain by the juxtamembrane segment. Cell 137, 1293–1307 (2009).
Bae, Y. S. et al. Epidermal growth factor (EGF)-induced generation of hydrogen peroxide: role in EGF receptor-mediated tyrosine phosphorylation. J. Biol. Chem. 272, 217–221 (1997).
Sundaresan, M., Yu, Z. X., Ferrans, V. J., Irani, K. & Finkel, T. Requirement for generation of H2O2 for platelet-derived growth factor signal transduction. Science 270, 296–299 (1995).
Red Brewer, M. et al. The juxtamembrane region of the EGF receptor functions as an activation domain. Mol. Cell 34, 641–651 (2009).
Kovacs, E. et al. Analysis of the role of the C-Terminal tail in the regulation of the epidermal growth factor receptor. Mol. Cell Biol. 35, 3083–3102 (2015).
Autiero, M. et al. Role of PlGF in the intra- and intermolecular cross talk between the VEGF receptors Flt1 and Flk1. Nat. Med. 9, 936–943 (2003).
Huang, K., Andersson, C., Roomans, G. M., Ito, N. & Claesson-Welsh, L. Signaling properties of VEGF receptor-1 and −2 homo- and heterodimers. Int J. Biochem. Cell Biol. 33, 315–324 (2001).
Freed, D. M. et al. EGFR ligands differentially stabilize receptor dimers to specify signaling kinetics. Cell 171, 683–695 e618 (2017).
Huang, Y. et al. A molecular mechanism for the generation of ligand-dependent differential outputs by the epidermal growth factor receptor. Elife 10, e73218 (2021).
Lemmon, M. A., Freed, D. M., Schlessinger, J. & Kiyatkin, A. The dark side of cell signaling: positive roles for negative regulators. Cell 164, 1172–1184 (2016).
Kiyatkin, A., van Alderwerelt van Rosenburgh, I. K., Klein, D. E. & Lemmon, M. A. Kinetics of receptor tyrosine kinase activation define ERK signaling dynamics. Sci. Signal 13, eaaz5267 (2020).
Nguyen, T. T. et al. Co-regulation of the mitogen-activated protein kinase, extracellular signal-regulated kinase 1, and the 90-kDa ribosomal S6 kinase in PC12 cells. Distinct effects of the neurotrophic factor, nerve growth factor, and the mitogenic factor, epidermal growth factor. J. Biol. Chem. 268, 9803–9810 (1993).
Weddell, J. C., Chen, S. & Imoukhuede, P. I. VEGFR1 promotes cell migration and proliferation through PLCgamma and PI3K pathways. NPJ Syst. Biol. Appl 4, 1 (2018).
Asahara, T. et al. Isolation of putative progenitor endothelial cells for angiogenesis. Science 275, 964–967 (1997).
Saiki, R. K. et al. Primer-directed enzymatic amplification of DNA with a thermostable DNA polymerase. Science 239, 487–491 (1988).
Gibson, D. G. et al. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat. Methods 6, 343–345 (2009).
Hartig, S. M. Basic image analysis and manipulation in ImageJ. Curr. Protoc. Mol. Biol. 14, Unit14 15 (2013).
Soumpasis, D. M. Theoretical analysis of fluorescence photobleaching recovery experiments. Biophys. J. 41, 95–97 (1983).
da Rocha-Azevedo, B. et al. Heterogeneity in VEGF Receptor-2 mobility and organization on the endothelial cell surface leads to diverse models of activation by VEGF. Cell Rep. 32, 108187 (2020).
Roy, A., Kucukural, A. & Zhang, Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat. Protoc. 5, 725–738 (2010).
Van Der Spoel, D. et al. GROMACS: fast, flexible, and free. J. Comput Chem. 26, 1701–1718 (2005).
Berendsen, H. J. C., Vanderspoel, D. & Vandrunen, R. Gromacs - a message-passing parallel molecular-dynamics implementation. Comput Phys. Commun. 91, 43–56 (1995).
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W. & Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79, 926–935 (1983).
Hess, B., Bekker, H., Berendsen, H. J. C. & Fraaije, J. G. E. M. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem. 18, 1463–1472 (1997).
Huang, J. et al. CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nat. Methods 14, 71–73 (2017).
Parrinello, M. & Rahman, A. Polymorphic transitions in single-crystals - a new molecular-dynamics method. J. Appl Phys. 52, 7182–7190 (1981).
Nose, S. & Klein, M. L. Constant pressure molecular-dynamics for molecular-systems. Mol. Phys. 50, 1055–1076 (1983).
Humphrey, W., Dalke, A. & Schulten, K. VMD: visual molecular dynamics. J. Mol. Graph 14, 33–38 (1996). 27-38.
Acknowledgements
The authors thank Dr. Arnab Gupta and his group members for access to the microscope facility. The authors thank Kaustav Gangopadhyay, Swarnendu Roy, and Subhankar Chowdhury for their help in the initial stage of the project and comments. The authors thank research funding from IISER Kolkata, infrastructural facilities supported by IISER Kolkata, and DST-FIST (SR/FST/LS-II/2017/93(c)). This work is supported by a grant from SERB (CRG/2020/000437). Fellowships from CSIR-UGC supported MPC and PM; Fellowship from KVPY supported DD.
Author information
Authors and Affiliations
Contributions
The manuscript was written through the contribution of all authors. All authors have approved the final version of the manuscript. R.D. and M.P.C. designed the experiments. M.P.C., P.M., P.K., S.B., and P.K.D. performed the Biochemical experiments, Imaging, and data analysis. D.D. did the molecular dynamics simulation and data analysis. R.D. and M.P.C. wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chakraborty, M.P., Das, D., Mondal, P. et al. Molecular basis of VEGFR1 autoinhibition at the plasma membrane. Nat Commun 15, 1346 (2024). https://doi.org/10.1038/s41467-024-45499-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-45499-2
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.