Identification of a carbohydrate recognition motif of purinergic receptors
Abstract
As a major class of biomolecules, carbohydrates play indispensable roles in various biological processes. However, it remains largely unknown how carbohydrates directly modulate important drug targets, such as G-protein coupled receptors (GPCRs). Here, we employed P2Y purinoceptor 14 (P2Y14), a drug target for inflammation and immune responses, to uncover the sugar nucleotide activation of GPCRs. Integrating molecular dynamics simulation with functional study, we identified the uridine diphosphate (UDP)-sugar-binding site on P2Y14, and revealed that a UDP-glucose might activate the receptor by bridging the transmembrane (TM) helices 2 and 7. Between TM2 and TM7 of P2Y14, a conserved salt bridging chain (K2.60-D2.64-K7.35-E7.36 [KDKE]) was identified to distinguish different UDP-sugars, including UDP-glucose, UDP-galactose, UDP-glucuronic acid, and UDP-N-acetylglucosamine. We identified the KDKE chain as a conserved functional motif of sugar binding for both P2Y14 and P2Y purinoceptor 12 (P2Y12), and then designed three sugar nucleotides as agonists of P2Y12. These results not only expand our understanding for activation of purinergic receptors but also provide insights for the carbohydrate drug development for GPCRs.
Editor's evaluation
This study describes a valuable model for the interaction of nucleotides, an important group of signalling molecules in health and disease, with their receptors. By combining experimental and theoretical methods, the authors provide compelling evidence for the mechanism of receptor activation. The results will be useful for drug design to target these receptors.
https://doi.org/10.7554/eLife.85449.sa0eLife digest
Sugars and other types of carbohydrates are biomolecules which play a range of key roles in the body. In particular, they are important messengers that help to coordinate immune responses. For example, a carbohydrate known as UDP-Glucose (a kind of UDP-sugar) can activate P2Y14, a receptor studded through the surface of many cells; this event then triggers a cascade of molecular events associated with asthma, kidney injury and lung inflammation. Yet it remains unclear how exactly UDP-Glucose recognizes P2Y14 – and, more broadly, how carbohydrates interact with purinergic receptors, the class of proteins that P2Y14 belongs to.
To examine this question, Zhao et al. combined functional experiments in the laboratory with molecular dynamics simulations, a computational approach. This work revealed that UDP-Glucose may activate P2Y14 by bridging its segments anchored within the cell membrane. A component of P2Y14, known as the KDKE chain, was found to have an important role in distinguishing between highly similar types of UDP-sugars. This allowed Zhao et al. to design three sugar molecules which could activate another purinergic receptor that also contained a KDKE chain.
Purinergic receptors are promising therapeutic targets. A finer understanding of how they recognise the molecules that activate them is therefore important to be able to identify and design new drug compounds.
Introduction
As significant components of the organism, carbohydrates play indispensable roles in energy supply, cell signaling, and immune responses (Gagneux and Varki, 1999). Dysregulation of carbohydrates has been proved to be associated with the development of various diseases (Reily et al., 2019). However, it is still elusive how carbohydrates directly act on major therapeutic targets, including G-protein coupled receptors (GPCRs) (Cheng and Jiang, 2019; Hauser et al., 2017). P2Y purinoceptor 14 (P2Y14) represents an outstanding model system for understanding carbohydrate modulation of GPCRs. It belongs to P2Y purinoceptor subfamily, consisting of receptors responding to nucleotides, including adenosine diphosphate (ADP) and UDP (Ralevic and Burnstock, 1998). Distinct from the other purinoceptors, P2Y14 is potently activated by UDP and a class of carbohydrates, that is, UDP-sugars (Abbracchio et al., 2006; Jacobson et al., 2020). UDP-sugars activate P2Y14 with a relative potency order of UDP-glucose (UDP-Glc), UDP-galactose (UDP-Gal), UDP-glucuronic acid (UDP-GlcA), and UDP-N-acetylglucosamine (UDP-GlcNAc) (Chambers et al., 2000; Hamel et al., 2011; Ko et al., 2009; Ko et al., 2007). These sugar nucleotides act as important signaling molecules via P2Y14 to mediate many physiological processes (Amison et al., 2017; Breton and Brown, 2018; Ferreira et al., 2017; Lazarowski, 2010; Müller et al., 2005; Sesma et al., 2016). Particularly, UDP-Glc regulates immune responses and associate with asthma, kidney injury, and lung inflammation (Amison et al., 2017; Breton and Brown, 2018; Ferreira et al., 2017; Müller et al., 2005; Sesma et al., 2016). As an isomer of UDP-Glc, UDP-Gal is present in various cell models, including physiologically relevant primary cultures of human bronchial epithelial cells (Lazarowski, 2010). It remains unknown how these sugar nucleotides are recognized by P2Y14.
As the closest homolog to P2Y14, P2Y purinoceptor 12 (P2Y12) has not been reported to be activated by any sugar nucleotide (Jacobson et al., 2020; Ralevic and Burnstock, 1998). P2Y12 is potently activated ADP. The reported agonist-bound structures of P2Y12 provide insights to understand the nucleotide activation of P2Y purinoceptors. The crystal structures of P2Y12 show that a full agonist 2-methylthio-adenosine-5’-diphosphate (2MeSADP, a close analogue of ADP) binds to an extracellular pocket consisting of transmembrane (TM) helices and extracellular loops (Zhang et al., 2014a). Since P2Y12 is highly similar to P2Y14 with 45.67% amino acid sequence identity, it would be interesting to investigate whether this receptor is also sensible to sugar nucleotides.
Here, we combined molecular docking, molecular dynamics (MD) simulations, and functional study to reveal the molecular mechanism how P2Y14 is activated by a sugar nucleotide. The ligand-binding models of different UDP-sugars (UDP-Glc, UDP-Gal, UDP-GlcA, and UDP-GlcNAc) were quantitatively characterized to identify the sugar recognition site of P2Y14. Both P2Y14 and P2Y12 were employed to unveil a conserved sugar-binding motif. Multiple carbohydrates were designed and validated as their agonists targeting the conserved functional motif.
Results
Identification of sugar-binding site in P2Y14
Both UDP and UDP-Glc potently activate P2Y14 with EC50 values of 50.9±6.1 nM and 40.3±1.5 nM, respectively (Figure 1A–C). Compared with UDP, UDP-Glc showed an increased potency on P2Y14 at high concentration (Figure 1B–C), suggesting that the sugar moiety of UDP-Glc contributes to activating P2Y14. To investigate how UDP-Glc regulates the P2Y14 via its sugar moiety, we used molecular docking to construct UDP-Glc-bound models of P2Y14 and compared them with UDP-bound P2Y14 models (Figure 1D–G). Because the protein structure of P2Y14 is unrevealed, we employed the X-ray structures of P2Y12 (Zhang et al., 2014a) as templates to constructed homology models of human P2Y14. The molecular docking showed that both UDP and UDP-Glc bind to an extracellular pocket consisting of TM helices 2–7 and extracellular loop (ECL) 2 (Figure 1D), which is corresponding to a known agonist binding pocket of P2Y12 (Zhang et al., 2014a). The docking score of UDP-Glc is –9.3 kcal/mol, and that of UDP is –8.5 kcal/mol, indicating that both UDP and UDP-Glc bind to P2Y14.
![](https://iiif.elifesciences.org/lax:85449%2Felife-85449-fig1-v1.tif/full/617,/0/default.jpg)
Identification of sugar-binding sites in P2Y purinoceptor 14 (P2Y14).
(A) Chemical structures of uridine diphosphate (UDP) and UDP-glucose (UDP-Glc). (B) Concentration-response curves of calcium mobilization for UDP or UDP-Glc in HEK293 cells transiently co-transfected with human P2Y14 and Gαqi5. Data are shown as mean ± SEM (n=3). See also Figure 1—source data 1 and Supplementary file 1. (C) Concentration-dependent comparison of calcium mobilization for UDP and UDP-Glc in HEK293 cells transiently co-transfected with human P2Y14 and Gαqi5 (n=3); **p<0.01, ***p<0.001, ****p<0.0001. (D) Ligand-bound model of P2Y14. Protein and compound are shown in cartoon and stick representations. (E, F) Docking models of UDP (E) and UDP-Glc (F) to P2Y14. Key residues are highlighted in yellow. Two sub-pockets for ligand binding are marked with circles. (G) Decomposition of ligand-binding energy for each receptor residue (n=10). (H) Calcium mobilization concentration-response curves for UDP or UDP-Glc in HEK293 expressing P2Y14 wild-type (WT) and mutants (n=3). See also Figure 1—source data 2 and Supplementary file 1. (I) Comparison of EC50s for UDP-Glc and UDP in HEK293 cells expressing P2Y14 mutants in calcium mobilization assay (n=3); **p<0.01, ***p<0.001. See also Figure 1—source data 3.
-
Figure 1—source data 1
Potency of uridine diphosphate (UDP) or UDP-glucose (UDP-Glc) in HEK293 cells expressing P2Y purinoceptor 14 (P2Y14).
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig1-data1-v1.docx
-
Figure 1—source data 2
Potency of uridine diphosphate (UDP) or UDP-glucose (UDP-Glc) in HEK293 cells expressing P2Y purinoceptor 14 (P2Y14) wild-type (WT) and mutants.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig1-data2-v1.docx
-
Figure 1—source data 3
Comparison of EC50s for uridine diphosphate-glucose (UDP-Glc) and UDP in HEK293 cells expressing P2Y purinoceptor 14 (P2Y14) mutants.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig1-data3-v1.docx
Compared with the UDP-bound receptor model (Figure 1E), the UDP-Glc-bound model showed extra interactions between the glucose moiety and the TM2, TM3, TM7, and ECL2 of P2Y14 (Figure 1F), enhancing the binding of UDP-Glc. Based on these molecular docking models, we further decomposed the ligand-binding energy to each receptor residue (Figure 1G). Five residues (K772.60, D812.64, N903.21, I170ECL2, and E2787.36; superscript indicates Ballesteros-Weinstein residue numbering; Ballesteros and Weinstein, 1995) were predicted to stabilize UDP-Glc binding (Figure 1F and G), while they made few energetic contributions (van der Waals energy >–0.25 Rosetta energy unit) to UDP binding (Figure 1E and G). As shown in Figure 1E and F, two sub-pockets of P2Y14 were unveiled for ligand binding. The sub-pocket 1 is formed by 16 residues of TMs 3–7 and ECL2 (Figure 1G, Figure 1—figure supplement 1A,B) and binds to the nucleotide moiety of the agonist, that is, UDP. The sub-pocket 2 is the specific sugar-binding site involving K772.60, D812.64, N903.21, I170ECL2, and E2787.36 (Figure 1F and G). These residues are primarily charged or polar amino acids, which could made hydrogen bonds with the glucose hydroxyl groups of UDP-Glc (Figure 1F). To validate the proposed sugar-binding sites, we designed single-point mutations of these five residues (K77A, D81A, N90A, I170A, and E278A). Among these mutants, D81A showed significantly reduced activities by UDP-Glc compared with the wild-type (WT) group (Figure 1H). However, substitution of D812.64 by alanine did not significantly affect the receptor activities by UDP (Figure 1H). Interestingly, K77A mutation diminished both UDP-Glc- and UDP-induced calcium mobilization (Figure 1H), but it showed greater impact on UDP-Glc-induced receptor responses than UDP-induced ones (Figure 1H, I), suggesting extra interactions between K772.60 and sugar moiety of UDP-Glc. Two mutations on TM3 and ECL2 (N90A and I170A) did not significantly affect the receptor responses by UDP or UDP-Glc (Figure 1—figure supplement 1C, D). These findings indicate that sub-pocket 2 residues of TM2 provide major contributions to stabilizing the sugar moiety of UDP-Glc.
UDP-Glc as a ‘glue’ for P2Y14 activation
The molecular docking employs rigid side chains of the receptor and does not include the influence of explicit environment on molecular interactions. To investigate how UDP-Glc interacts with P2Y14, we performed all-atom MD simulations of the P2Y14 receptor with and without UDP-Glc (Figure 2A). We used the molecular docking model of P2Y14 to construct the simulation systems. Each system was replicated to performed three independent simulations (see Materials and methods for detailed information). Apo P2Y14 and UDP-Glc-bound P2Y14 simulation models showed different conformations in TM6 and TM7 (Figure 2A and B). In UDP-Glc-bound P2Y14 simulations, the extracellular tip of TM6 shifted over 3 Å and TM7 over 4 Å toward the receptor core, compared with the apo P2Y14 simulations (Figure 2A and B). This inward shift of TM6 and TM7 allowed formation of polar and ionic interactions with the UDP-Glc (Figure 2B and C). During UDP-Glc-bound P2Y14 simulations, two charged residues K2777.35 and E2787.36 formed hydrogen bonds with the glucose 6’ hydroxyl group of UDP-Glc to keep TM7 close to the receptor core (Figure 2B and C), and an arginine residue (R2536.55) formed a salt bridge with the phosphate group of UDP-Glc to stabilize the inward shift of TM6 (Figure 2B and C). Consistently, compared with WT group (EC50 of 40.3±1.5 nM), single-point mutations (R253A and E278A) of TM6 and TM7 helices resulted to diminished UDP-Glc-induced calcium mobilization (EC50 of 808.6±43.6 nM for R253A and 60.2±3.6 nM for E278A) (Figure 1H, Figure 2D). In addition, at the extracellular side, the distance between TM5 and TM6 of UDP-Glc-bound P2Y14 was 5.9 Å shorter than that in the apo system (Figure 2—figure supplement 1A, B). Y1895.41 and T2576.59 made stable hydrophobic interactions to maintain the tight compact between TM5 and TM6 in UDP-Glc-bound-P2Y14 simulations, while TM6 did not interact with TM5 at the extracellular side in the apo simulations (Figure 2—figure supplement 1A, B). Compared with WT group (EC50 of 40.3±1.5 nM), a mutation of TM6 (T257A) showed significantly reduced UDP-Glc-induced responses (EC50 of 504.9±15.9 nM) (Figure 2—figure supplement 1C), fully agreeing with our simulation models. Collectively, these data suggest that UDP-Glc might serve as intramolecular ‘glue’ to make a tight helical bundle of P2Y14, involving TM6 and TM7.
![](https://iiif.elifesciences.org/lax:85449%2Felife-85449-fig2-v1.tif/full/617,/0/default.jpg)
Comparison of the apo P2Y purinoceptor 14 (P2Y14) and uridine diphosphate-glucose (UDP-Glc)-bound P2Y14 simulation systems.
(A–B) Side (A) and top (B) views of representative models of apo P2Y14 and UDP-Glc-bound P2Y14. Receptor is shown as cartoon. Ligand and key residues are shown as sticks. Movement of the extracellular tips of transmembrane (TM)6 and TM7 toward the receptor core is shown by arrows. See Supplementary file 2 for computational characterization of conformational changes. (C) Key interactions between P2Y14 and UDP-Glc. Hydrogen bonds and salt bridges are displayed as red dashed lines. See Supplementary file 2 for pairwise interaction details. (D) Concentration-response curves of calcium mobilization for UDP-Glc in HEK293 expressing P2Y14 wild-type (WT) and mutants. Data are shown as mean ± SEM (n=3). See also Figure 2—source data 1 and Supplementary file 1.
-
Figure 2—source data 1
Potency of uridine diphosphate-glucose (UDP-Glc) in HEK293 expressing P2Y purinoceptor 14 (P2Y14) wild-type (WT) and mutants.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig2-data1-v1.docx
Molecular recognition of P2Y14 via sugar-binding site
P2Y14 could be activated by different UDP-sugars with distinct potencies. With only one group substitution at the sugar moiety, UDP-Glc induced stronger activity on P2Y14 (EC50=40.3 ± 1.5 nM) than the other UDP-sugars, that is, UDP-Gal (EC50=78.3 ± 9.2 nM), UDP-GlcA (EC50=59.9 ± 4.8 nM), and UDP-GlcNAc (EC50=184.4 ± 11.8 nM) (Figure 3A and B). To investigate how P2Y14 recognizes different sugar moieties, we performed MD simulations of the human P2Y14 receptor complex with UDP-Gal, UDP-GlcA, and UDP-GlcNAc, respectively, and compared them with the UDP-Glc-bound P2Y14 simulations. We observed that UDP-Gal, UDP-GlcA, and UDP-GlcNAc bound to P2Y14 at the same pocket as UDP-Glc. Similar to UDP-Glc, their uridine groups occupied the sub-pocket 1 of P2Y14, while their sugar moieties bound to the sub-pocket 2 during simulations (Figure 3C–F, Figure 3—figure supplements 1 and 2). At the sub-pocket 2, a stable salt bridging chain formed by four charged residues (K772.60, D812.64, K2777.35, and E2787.36) were observed in all systems (Figure 3C–F). The negative charged glutamic acid residue E2787.36 linked TM2 and TM7 helices by forming salt bridges with K772.60 and K2777.35, while the other negative charged residue D812.64 forming a salt bridge with K772.60 to further stabilize these ionic interactions (Figure 3C–F).
![](https://iiif.elifesciences.org/lax:85449%2Felife-85449-fig3-v1.tif/full/617,/0/default.jpg)
Sugar moiety recognition of P2Y purinoceptor 14 (P2Y14).
(A) Chemical structures of uridine diphosphate-glucose (UDP-Glc), UDP-glucuronic acid (GlcA), UDP-N-acetylglucosamine (GlcNAc), and UDP-galactose (Gal). (B) Concentration-response curves of calcium mobilization for different UDP-sugars in HEK293 cells transiently co-transfected with human P2Y14 and Gαqi5. Data are shown as mean ± SEM (n=3). See also Figure 3—source data 1. (C–F) Molecular recognition of P2Y14 for UDP-Glc (C), UDP-GlcA (D), UDP-GlcNAc (E), and UDP-Gal (F). Receptor, ligands, and key residues are shown in cartoon and stick representations. Hydrogen bonds and salt bridges are displayed as red dashed lines. See Supplementary file 2 for pairwise interaction details. (G–I) Concentration-response curves of calcium mobilization for UDP-GlcA (G), UDP-GlcNAc (H), and UDP-Gal (I) in HEK293 expressing P2Y14 wild-type (WT) and mutants. Data are shown as mean ± SEM (n=3). See also Figure 3—source data 2.
-
Figure 3—source data 1
Potency of uridine diphosphate-glucuronic acid (UDP-GlcA), UDP-N-acetylglucosamine (GlcNAc), and UDP-galactose (Gal) in HEK293 expressing P2Y purinoceptor 14 (P2Y14).
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig3-data1-v1.docx
-
Figure 3—source data 2
Potency of uridine diphosphate-glucuronic acid (UDP-GlcA), UDP-N-acetylglucosamine (GlcNAc), and UDP-galactose (Gal) in HEK293 expressing P2Y purinoceptor 14 (P2Y14) wild-type (WT) and mutants.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig3-data2-v1.docx
In simulations, different sugar moieties bound to the K772.60-D812.64-K2777.35-E2787.36 salt bridging chain with distinct binding modes (Figure 3C–F, Figure 3—figure supplement 3). For UDP-Glc, both K2777.35 and E2787.36 could form hydrogen bonds with the glucose 6’ hydroxyl group to keep TM7 close to the receptor core (Figure 3C, Figure 3—figure supplement 4A). However, in UDP-GlcA-bound P2Y14 simulations, at the corresponding position, the 5’ carboxyl group of the sugar moiety failed to form hydrogen bond with the negatively charged E2787.36 (Figure 3D, Figure 3—figure supplement 4B). Compared with that of WT group (EC50=59.9 ± 4.8 nM), the single-point mutation of E278A did not decrease the UDP-GlcA-induced calcium mobilization (EC50=38.2 ± 2.2 nM) (Figure 3G), supporting with the proposed sugar-binding model (Figure 3D). Substitution of 2’ hydroxyl group by an acetamido group led to a rotation of the sugar moiety of UDP-GlcNAc in simulations (Figure 3E). Consequentially, the 6’ hydroxyl of N-acetylglucosamine group flipped to form hydrogen bonds with K772.60 and D812.64 instead of K2777.35 and E2787.36 (Figure 3E, Figure 3—figure supplement 4C). Consistently, single-point mutation of D81A made significant effect to reduce the UDP-GlcNAc-induced receptor activities (Figure 3H). Compared with the other three UDP-sugars, UDP-Gal has a different orientation of 4’ hydroxyl group. The 4’ hydroxyl group of galactose formed a stable hydrogen bond with K772.60 and disrupted the interaction between 6’ hydroxyl group with E2787.36 (Figure 3F, Figure 3—figure supplement 4D). Compared with UDP-Glc, UDP-Gal had more interactions with TM2 and less interactions with TM7 (E2787.36). Substitution of E2787.36 by alanine did not significantly affect the UDP-Gal-induced receptor response (Figure 3I), agreeing with the proposed UDP-Gal-binding model (Figure 3F). For all UDP-sugars, at least three residues of K772.60, D812.64, K2777.35, and E2787.36 participated in ligand binding (Figure 2, Figure 3). Both computational models and experimental data indicate the K772.60-D812.64-K2777.35-E2787.36 salt bridging chain as a sugar-binding site of P2Y14, which can recognize different sugar moieties. The interactions of ligands with the TM7 might determine the ligand potency on P2Y14.
Conserved sugar-binding motif for P2Y12 and P2Y14
In previous sections, we have identified K772.60-D812.64-K2777.35-E2787.36 salt bridging chain as an important functional site for sugar moiety recognition and UDP-sugar activation of P2Y14. These four residues (K2.60, D2.64, K7.35, and E7.36) are conserved between P2Y14 and its closest homolog, that is, P2Y12 (Figure 4A, Figure 4—figure supplement 1A). Therefore, we assumed that P2Y12 can also be activated by carbohydrate ligands. P2Y12 is activated by ADP (Herbert and Savi, 2003), but it has not been reported to be activated by any sugar nucleotide. To validate our assumption, we designed and synthesized three ADP-sugars, that is, ADP-glucose (ADP-Glc), ADP-glucuronic acid (ADP-GlcA), and ADP-mannose (ADP-Man), and then tested whether they can activate P2Y12 (Figure 4B and C). We docked ADP-Glc and ADP-Man to the X-ray structure of P2Y12 (Zhang et al., 2014a). The docking scores are –9.4 kcal/mol for ADP-Glc, for –10.0 kcal/mol for ADP-GlcA, and –9.3 kcal/mol for ADP-Man (Figure 4—figure supplement 1B), suggesting they stably bound to P2Y12. Consistently, in calcium mobilization assays, ADP-Glc, ADP-GlcA, and ADP-Man activated P2Y12 with EC50 values of 3.4±0.4 μM, 1.3±0.1 μM, and 12.3±0.9 μM, respectively (Figure 4C). Single-point mutations of K802.60, D842.64, K2807.35, and E2817.36 significantly diminished ADP-Glc-, ADP-GlcA-, and ADP-Man-induced responses, compared with WT P2Y12 (Figure 4C). These findings not only validate our assumption that P2Y12 can be activated by sugar nucleotides, but also indicate the conserved KDKE salt bridging chain as a functional motif for sugar binding. In addition, we found a less conserved phenylalanine residue F2777.32, which is adjacent to the KDKE motif of P2Y12 (Figure 4—figure supplement 2A). The alanine substitution of F2777.32 significantly reduced ADP-Glc-induced responses of P2Y12, compared with that of the WT group (Figure 4—figure supplement 2B, C). However, single-point mutations at the corresponding position of P2Y14 did not significantly affect the UDP-Glc-induced activation (Figure 4—figure supplement 2C–F). These findings indicate different nature of P2Y12 and P2Y14.
![](https://iiif.elifesciences.org/lax:85449%2Felife-85449-fig4-v1.tif/full/617,/0/default.jpg)
Adenosine diphosphate (ADP)-sugars binding to P2Y purinoceptor 12 (P2Y12).
(A) Sequence log of the alignment between P2Y12 and P2Y14. Sequences of P2Y12 and P2Y14 involving 359 species were collected for making sequence alignments. See Supplementary file 3 and Figure 4—source data 1 for species repertoire information. The height of a letter is proportional to the relative frequency of that residue at a particular site. Four residues of K2.60-D2.64-K7.35-E7.36 (KDKE) sugar-binding motif are marked by arrows, with the corresponding residues in P2Y12 and P2Y14. (B) Chemical structure of ADP-glucose (Glc), ADP-glucuronic acid (GlcA), and ADP-mannose (Man). (C) Calcium mobilization concentration-response curves for ADP-Glc, ADP-GlcA, and ADP-Man in HEK293 expressing P2Y12 wild-type (WT) and mutants. Data are shown as mean ± SEM (n=3). See also Figure 4—source data 1 and Supplementary file 1.
-
Figure 4—source data 1
Potency of adenosine diphosphate-glucose (ADP-Glc), ADP-glucuronic acid (GlcA), and ADP-mannose (Man) in HEK293 expressing P2Y purinoceptor 12 (P2Y12) wild-type (WT) and mutants.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig4-data1-v1.docx
-
Figure 4—source data 2
The neighbor joining tree of P2Y purinoceptor 12 (P2Y12) showing evolutionary range of species.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig4-data2-v1.docx
-
Figure 4—source data 3
The neighbor joining tree of P2Y purinoceptor 13 (P2Y13) showing evolutionary range of species.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig4-data3-v1.docx
-
Figure 4—source data 4
The neighbor joining tree of P2Y purinoceptor 14 (P2Y14) showing evolutionary range of species.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig4-data4-v1.docx
-
Figure 4—source data 5
The neighbor joining tree of GPR87 showing evolutionary range of species.
- https://cdn.elifesciences.org/articles/85449/elife-85449-fig4-data5-v1.docx
Discussion
Mediated by UDP-sugars, P2Y14 plays an important role in immune responses and inflammation (Arase et al., 2009; Barrett et al., 2013; Breton and Brown, 2018; Ferreira et al., 2017; Müller et al., 2005; Sesma et al., 2012; Sesma et al., 2016), and possibly insulin resistance (Wang et al., 2008). Breton et al. found that kidney collecting duct intercalated cells present high levels of P2Y14, which is activated by UDP-Glc to promote neutrophil infiltration and renal inflammation (Breton and Brown, 2018). Exerting excessive P2Y14-mediated inflammatory reactions, high concentration of UDP-sugars was observed in extracellular tissue surrounding airway epithelial cells and lung secretions of cystic fibrosis patients (Ferreira et al., 2017; Müller et al., 2005; Sesma et al., 2016). UDP-Glc is also released from liver cells in obese states, possibly via hepatocellular apoptosis, leading to liver inflammation and insulin resistance (Wang et al., 2008). All these results indicate importance of UDP-sugar regulation of P2Y14 in pathological progresses.
In this work, we built the molecular model of UDP-Glc-bound P2Y14 to answer the long-standing question of sugar nucleotide regulation of the purinergic receptor. Binding to an extracellular pocket involving TMs 2 and 7 (Figure 1), the UDP-Glc might serve as intramolecular ‘glue’ attaching to TM6 and TM7 to activate P2Y14 (Figure 2). The agonist-induced remarkable conformational changes of TM6 and TM7 are also reported for P2Y12 (Zhang et al., 2014a). Compared with the AZD1283-bound (antagonist-bound) P2Y12 structure (Zhang et al., 2014b), the extracellular part of TM6 in the 2MeSADP-bound (agonist-bound) P2Y12 structure shifts over 10 Å and TM7 over 5 Å toward the center of TM helix bundle (Zhang et al., 2014a). The close parallel of P2Y12 and P2Y14 in the agonist-induced conformational changes indicates a common ligand-induced activation mechanism shared by purinergic receptors. In addition, in the studies involving the other UDP-sugars, we also found that the interactions between sugar-moieties of agonists with TM7 (E2787.36) is determinant for UDP-sugars’ potencies (Figure 3).
The carbohydrate-binding site has not been fully characterized for GPCRs. Except for P2Y14, it has not been reported that the other members of P2Y12-like subfamily can be directly activated by carbohydrates. Integrated computational modeling with mutagenesis study, we identified a conserved carbohydrate-binding motif (KDKE) for both P2Y14 and P2Y12 (Figure 4). The KDKE motif not only participates in receptor activation by bridging TM2 and TM7 (K772.60, D812.64, K2777.35, and E2787.36) (Figure 2), but also recognize different sugar moieties, including glucose, galactose, glucuronic acid, and N-acetylglucosamine groups (Figure 3, Figure 4). Remarkably, this KDKE motif can distinguish isomers as UDP-Glc and UDP-Gal. Our MD simulations showed that KDKE motif attracted the 6’ hydroxyl group of glucose and interacted with the 4’ hydroxyl group of galactose (Figure 3C and F). Consistent with our observations, a previous structure-activity relationship study revealed that selective mono-fluorination of the 6’ hydroxyl group of the glucose moiety results to fourfold less potency on P2Y14 (Ko et al., 2009). As another member of P2Y12-like subfamily, P2Y13 also has the conserved KDKE site (Figure 4—figure supplement 3), suggesting it might be regulated by carbohydrates. GPR87 is a close homolog of P2Y14 with the sequence identity of 44.94%. GPR87 has K/R2.60, D2.64, K/E7.35, and E7.36 at the corresponding positions of the KDKE sugar-binding motif, indicating varied carbohydrate sensitivities of this receptor in different species (Figure 4—figure supplement 4). The species repertoires of these receptors consist of many amniotes but a few anamniotes (Supplementary file 3a). Distant relatives of the amniote receptors might show different sensitivities to carbohydrates. Similar to the KDKE motif of the receptors, the UDP-sugar-binding sites consisting of charged residues have been discovered for sugar transferases (Gerlach et al., 2018; Hao et al., 2021). In typical glycosylation transfers as TarP and SseK3, two aspartic acids and one positively charged residue (arginine or lysine) participate in recognition of 3’ or 4’ hydroxyl groups of GlcNAc or GalNAc moiety (Gerlach et al., 2018; Hao et al., 2021). However, a salt bridging chain has not been observed in these sugar-binding sites. The different arrangements of UDP-sugar-binding sites between P2Y14 and these sugar transferases might be determinant for their sugar selectivity.
In conclusion, we revealed a conserved carbohydrate-binding motif in both P2Y12 and P2Y14, extending our understanding of how carbohydrates regulate GPCRs. Our molecular models of different sugar nucleotides provide great details for carbohydrate activation and recognition of these receptors, which would inspire further carbohydrate drug development for GPCRs. Whether the other carbohydrate-binding motifs exist in GPCRs is currently unknown. Further investigations focused on carbohydrate regulation of GPCRs will continue to add both new concepts and physiological understanding to the field.
Materials and methods
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Gene (Homo sapiens) | P2Y14 | GenBank | NM_001081455.2 | |
Gene (Homo sapiens) | P2Y12 | GenBank | NM_022788.5 | Optimized |
Strain, strain background (Escherichia coli) | Trans5α Chemically Competent Cell | TransGen Biotech | Cat.No: CD201-01 | |
Cell line (Homo sapiens) | HEK293 | ATCC | CRL-1573 | |
Antibody | Anti-HA primary antibody (Rabit monoclonal) | Cell Signaling Technology | Cat.No: 3724. | FCM (1:800) |
Antibody | Goat anti-rabbit IgG(H+L) FITC conjugate secondary antibody (Goat monoclonal) | TransGen Biotech | Cat.No: HS111 | FCM (1:200) |
Recombinant DNA reagent | pCDNA3-HA(plasmid) | This paper | HA version of pCDNA3 | |
Sequence-based reagent | P2Y14-K77A-F | This paper | PCR primers | GACTTTTCCT TTCGCGATCC TTGGTGAC |
Sequence-based reagent | P2Y14-K77A-R | This paper | PCR primers | GTCACCAAGG ATCGCGAAAG GAAAAGTC |
Sequence-based reagent | P2Y14-D81A-F | This paper | PCR primers | CAAGATCCT TGGTGCCTC AGGCCTTGG |
Sequence-based reagent | P2Y14-D81A-R | This paper | PCR primers | GACCAAGGCC TGAGGCACCAA GGATCTTG |
Sequence-based reagent | P2Y14-N90A-F | This paper | PCR primers | GTCCCTGGCA GCTGGCCGTG TTTGTGTGCAG |
Sequence-based reagent | P2Y14-N90A-R | This paper | PCR primers | CTGCACACAAA CACGGCCAGCT GCCAGGGAC |
Sequence-based reagent | P2Y14-I170A-F | This paper | PCR primers | GAGGTTACACA AGCAAAATGTA TAGAACTG |
Sequence-based reagent | P2Y14-I170A-R | This paper | PCR primers | GTTCTATACATT TTGCTTGTGTAA CCTC |
Sequence-based reagent | P2Y14-R274A-F | This paper | PCR primers | CAAAAGAAA TCTTGGCGT ATATGAAAG AATTC |
Sequence-based reagent | P2Y14-R274A-R | This paper | PCR primers | GAATTCTTTC ATATACGCCA AGATTTCTTTT G |
Sequence-based reagent | P2Y14-K277A-F | This paper | PCR primers | CTTGCGGTAT ATGGCAGAAT TCACTCTG |
Sequence-based reagent | P2Y14-K277A-R | This paper | PCR primers | CAGAGTGAAT TCTGCCATAT ACCGCAAG |
Sequence-based reagent | P2Y14-E278A-F | This paper | PCR primers | GCGGTATAT GAAAGCATT CACTCTGCT AC |
Sequence-based reagent | P2Y14-E278A-R | This paper | PCR primers | GTAGCAGAG TGAATGCTTT CATATACCG |
Sequence-based reagent | P2Y12-K80A-F | This paper | PCR primers | CACATTCCC ATTCGCGAT CCTGTCAGA TG |
Sequence-based reagent | P2Y12-K80A-R | This paper | PCR primers | CATCTGACAG GATCGCGAAT GGGAATGTG |
Sequence-based reagent | P2Y12-D84A-F | This paper | PCR primers | CAAGATCCTGT CAGCTGCCAAG CTCGGTAC |
Sequence-based reagent | P2Y12-D84A-R | This paper | PCR primers | GTACCGAGCTT GGCAGCTGACA GGATCTTG |
Sequence-based reagent | P2Y12-F277A-F | This paper | PCR primers | GAGAACACTC TGGCCTACGT CAAGGAATC |
Sequence-based reagent | P2Y12-F277A-R | This paper | PCR primers | GATTCCTTGAC GTAGGCCAGAG TGTTCTC |
Sequence-based reagent | P2Y12-K280A-F | This paper | PCR primers | CTGTTCTACGT CGCGGAATCCA CATTG |
Sequence-based reagent | P2Y12-K280A-R | This paper | PCR primers | CAATGTGGATT CCGCGACGTAG AACAG |
Sequence-based reagent | P2Y12-E281A-F | This paper | PCR primers | GTTCTACGTCA AGGCATCCACA TTGTGGC |
Sequence-based reagent | P2Y12-E281A-R | This paper | PCR primers | GCCACAATGTG GATGCCTTGAC GTAGAAC |
Commercial assay or kit | KOD-plus-Ver.2 | TOYOBO | Cat.No: KOD-211 | |
Commercial assay or kit | MycoBlue Mycoplasma Detector | Vazyme | Cat.No: D101-01 | |
Software, algorithm | GraphPad Prism 6 | GraphPad Prism 6 |
Chemicals
UDP-GlcNAc was prepared from D-GlcNAc as reported previously (Zheng et al., 2022). UDP-Glc and UDP-GlcA were prepared from Sucrose (Wang et al., 2022). UDP-Gal was prepared from D-Gal (Muthana et al., 2012). ADP-Man was synthesized by a two-step strategy. In detail, Man-1-p was first synthesized from D-Man using NahK from Bifidobacterium longum (Nishimoto and Kitaoka, 2007) and ATP as phosphorylation donor. Man-1-p was purified from the reaction mixture by the silver nitrate precipitation method (Wen et al., 2016; Wen et al., 2015). Then, ADP-Man was synthesized from Man-1-p and ATP by a GDP-mannose pyrophosphorylase from Pyrococcus furiosus, which could take ATP as substrate.
Cell lines
Request a detailed protocolHEK293 were purchased from ATCC and the identity of the cell line was confirmed by carrying out fingerprinting (Shanghai Genening Biotechnologies Inc, Shanghai, China). In addition, HEK293 cells were tested using the MycoBlue Mycoplasma Detector Kit (Vazyme,China), indicating that the cells were not contaminated by Mycoplasma.
Cell culture and transient transfections
Request a detailed protocolHEK293 cells were cultured in Dulbecco’s modified Eagle’s medium with 10% fetal bovine serum. All cells were maintained at 37°C in humidified incubators with 5% CO2 and 95% air. Human P2Y14 or P2Y12 receptors and G protein α-subunit (Gαqi5) were transiently co-transfected with HEK293 cells using PolyJet In Vitro DNA Transfection Reagent (SignaGen) according to the manufacturer’s instructions. Thus, a mixture of 1 μg of receptor DNA and 1 μg of Gαqi5 DNA was used to transfect with the six-well plate cells at 90% confluency. HEK293 cells transiently expressing P2Y14 or P2Y12 receptor were subsequently used for the intracellular Ca2+ assays after 48 hr post-transfection.
Cell surface expression
Request a detailed protocolHuman P2Y14 or P2Y12 was cloned into a pcDNA3 vector with HA tag for expression in HEK293 cells. Mutants of P2Y14 or P2Y12 were constructed according to Fast Mutagenesis System (TransGen). Cell surface expression of P2Y14 or P2Y12 was analyzed by flow cytometry. HEK293 cells were transfected with pCDNA3-HA-P2Y14 or P2Y12 in six-well plate overnight. After having been incubated with rabbit anti-HA primary antibody (1:800, CST) for 1 hr at 4°C, the cells were incubated with goat anti-rabbit IgG(H+L) FITC conjugate secondary antibody (1:200, TransGen) for 50 min at 4°C. Data were collected with a flow cytometer (FACS Calibur, BD) and analyzed with FlowJo software.
Intracellular Ca2+ mobilization
Request a detailed protocolIntracellular Ca2+ assays were carried out as follows. HEK293 cells were seeded (80,000 cells/well) into Matrigel-coated 96-well plate 24 hr prior to assay. The cells were incubated with 2 μM Fluo-4 AM (Invitrogen) diluted in HBSS solution (0.4 g L−1 KCl, 0.12 g L−1 Na2HPO4·12H2O, 0.06 g L−1 KH2PO4, 0.35 g L−1 NaHCO3, 0.14 g L−1 CaCl2, 0.10 g L−1 MgCl2·6H2O, 0.05 g L−1 MgSO4, and 8.0 g L−1 NaCl) at 37°C for 50 min. After dye loading, the cells were treated with the compounds of interest. Then, calcium response (relative fluorescence unit) was measured using Flexstation 3 (Molecular Device) with fluorescence excitation made at 485 nm and emission at 525 nm.
Molecular modeling, docking, and energy decomposition
Request a detailed protocolUsing the crystal structures of agonist-bound P2Y12 (PDB codes 4PXZ, 4PY0) (Zhang et al., 2014a) as templates, we employed Modeller (Sali and Blundell, 1993) to construct the human P2Y14 models. The human P2Y12 models are also built using these P2Y12 crystal structures (PDB codes 4PXZ, 4PY0) (Zhang et al., 2014a). The models with the lowest root mean square deviations from their template structures were selected for further analysis. A ligand was docked to the receptor using Schodinger Glide software in SP mode with default parameters (Friesner et al., 2004). A pocket binding to the ligand with Glide G-scores below –6.5 kcal/mol was considered as a possible ligand-binding pocket. To involve receptor flexibility, we used RosettaLigand (Davis and Baker, 2009) to generate representative ligand-bound receptor models. After Rosetta-based docking, the top 1000 models with lowest binding energy score were selected. Then, they were further scored with the ligand-binding energy between ligand and receptor. The top 10 models with the lowest ligand-binding energy were selected for energy decomposition. The van der Waals energy of ligand binding was mapped to each receptor residue by residue_energy_breakdown utility (Davis and Baker, 2009). The model with the lowest ligand binding energy was used as the representative model.
Modeling and simulations
Request a detailed protocolTo build a simulation system, we place the molecular model into a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine lipid bilayer. The lipid-embedded complex model was solvated in periodic boundary condition box (80 Å × 80 Å × 120 Å) filled with TIP3P water molecules and 0.15 M KCl using CHARMM-GUI (Wu et al., 2014). Each system was replicated to perform three independent simulations. On the basis of the CHARMM36m all-atom force field (Guvench et al., 2011; Huang et al., 2017; MacKerell et al., 1998), MD simulations were conducted using GROMAS 5.1.4 (Hess et al., 2008; Van Der Spoel et al., 2005). After 100 ns equilibration, a 500 ns production run was carried out for each simulation. All productions were carried out in the NPT ensemble at temperature of 303.15 K and a pressure of 1 atm. Temperature and pressure were controlled using the velocity-rescale thermostat (Bussi et al., 2007) and the Parrinello-Rahman barostat with isotropic coupling (Aoki and Yonezawa, 1992), respectively. Equations of motion were integrated with a 2 fs time step; the LINCS algorithm was used to constrain bond length (Hess, 2008). Non-bonded pair lists were generated every 10 steps using distance cutoff of 1.4 nm. A cutoff of 1.2 nm was used for Lennard-Jones (excluding scales 1–4) interactions, which were smoothly switched off between 1 nm and 1.2 nm. Electrostatic interactions were computed using particle-mesh-Ewald algorithm with a real-space cutoff of 1.2 nm. The last 200 ns trajectory of each simulation was used to calculate average values.
Sequence analysis
Request a detailed protocolTo analyze the conservation of residual sites, we collected sequences of receptors from UniProt database involving 359 species. Among these, we only found five anamniotes, including Astyanax mexicanus, Xenopus tropicalis, Microcaecilia unicolor, Geotrypetes seraphini, and Xenopus laevis. See Supplementary file 3 for species repertoire information. The multiple sequence alignments were performed using Clustal Omega. Logplots generated for these alignments by WebLog. In each logplot, the height of a letter is proportional to the information content of an amino acid in bits, which was calculated by Equation 1 as follows:
where N is the number of all sequences and pi is the probability of the amino acid in all sequences. A large value of the unit bits indicates a high conservation of a particular site.
Statistics
Statistical analyses were performed using GraphPad Prism 6 (GraphPad Software). EC50 values for compounds were obtained from concentration-response curves by nonlinear regression analysis. Comparison of two compounds or two constructs was analyzed by unpaired t test to determine statistical difference. All statistical data are given as mean ± SEM of at least three independent experiments performed in duplicate or triplicate.
Data availability
All data generated or analyzed during this study are included in the manuscript and supplement files.
References
-
Lipopolysaccharide (LPS) induced pulmonary neutrophil recruitment and platelet activation is mediated via the P2Y1 and P2Y14 receptors in micePulmonary Pharmacology & Therapeutics 45:62–68.https://doi.org/10.1016/j.pupt.2017.05.005
-
Constant-pressure molecular-dynamics simulations of the crystal-smectic transition in systems of soft parallel spherocylindersPhysical Review. A, Atomic, Molecular, and Optical Physics 46:6541–6549.https://doi.org/10.1103/physreva.46.6541
-
Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptorsMethods Neurosci 25:366–428.
-
Novel proinflammatory function of renal intercalated cellsAnnals of Nutrition & Metabolism 72 Suppl 2:11–16.https://doi.org/10.1159/000488303
-
Canonical sampling through velocity rescalingThe Journal of Chemical Physics 126:014101.https://doi.org/10.1063/1.2408420
-
A G protein-coupled receptor for UDP-glucoseThe Journal of Biological Chemistry 275:10767–10771.https://doi.org/10.1074/jbc.275.15.10767
-
Allostery in drug developmentAdvances in Experimental Medicine and Biology 1163:1–23.https://doi.org/10.1007/978-981-13-8719-7_1
-
RosettaLigand docking with full ligand and receptor flexibilityJournal of Molecular Biology 385:381–392.https://doi.org/10.1016/j.jmb.2008.11.010
-
Gene-based analysis of regulatory variants identifies 4 putative novel asthma risk genes related to nucleotide synthesis and signalingThe Journal of Allergy and Clinical Immunology 139:1148–1157.https://doi.org/10.1016/j.jaci.2016.07.017
-
Glide: a new approach for rapid, accurate docking and scoringJournal of Medicinal Chemistry 47:1739–1749.https://doi.org/10.1021/jm0306430
-
CHARMM additive all-atom force field for carbohydrate derivatives and its utility in polysaccharide and carbohydrate–protein modelingJournal of Chemical Theory and Computation 7:3162–3180.https://doi.org/10.1021/ct200328p
-
Discovery of novel P2Y14 agonist and antagonist using conventional and nonconventional methodsJournal of Biomolecular Screening 16:1098–1105.https://doi.org/10.1177/1087057111415525
-
Structures and mechanism of human glycosyltransferase β1,3-N-acetylglucosaminyltransferase 2 (B3GNT2), an important player in immune homeostasisThe Journal of Biological Chemistry 296:100042.https://doi.org/10.1074/jbc.RA120.015306
-
Trends in GPCR drug discovery: new agents, targets and indicationsNature Reviews. Drug Discovery 16:829–842.https://doi.org/10.1038/nrd.2017.178
-
P2Y12, a new platelet ADP receptor, target of clopidogrelSeminars in Vascular Medicine 3:113–122.https://doi.org/10.1055/s-2003-40669
-
P-LINCS: a parallel linear constraint solver for molecular simulationJournal of Chemical Theory and Computation 4:116–122.https://doi.org/10.1021/ct700200b
-
GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulationJournal of Chemical Theory and Computation 4:435–447.https://doi.org/10.1021/ct700301q
-
Update of P2Y receptor pharmacology: IUPHAR Review 27British Journal of Pharmacology 177:2413–2433.https://doi.org/10.1111/bph.15005
-
Structure-activity relationship of uridine 5’-diphosphoglucose analogues as agonists of the human P2Y14 receptorJournal of Medicinal Chemistry 50:2030–2039.https://doi.org/10.1021/jm061222w
-
Molecular recognition in the P2Y(14) receptor: Probing the structurally permissive terminal sugar moiety of uridine-5’-diphosphoglucoseBioorganic & Medicinal Chemistry 17:5298–5311.https://doi.org/10.1016/j.bmc.2009.05.024
-
Quantification of extracellular UDP-galactoseAnalytical Biochemistry 396:23–29.https://doi.org/10.1016/j.ab.2009.08.022
-
All-atom empirical potential for molecular modeling and dynamics studies of proteinsThe Journal of Physical Chemistry. B 102:3586–3616.https://doi.org/10.1021/jp973084f
-
The P2Y14 receptor of airway epithelial cells: coupling to intracellular Ca2+ and IL-8 secretionAmerican Journal of Respiratory Cell and Molecular Biology 33:601–609.https://doi.org/10.1165/rcmb.2005-0181OC
-
Identification of N-acetylhexosamine 1-kinase in the complete lacto-N-biose I/galacto-N-biose metabolic pathway in Bifidobacterium longumApplied and Environmental Microbiology 73:6444–6449.https://doi.org/10.1128/AEM.01425-07
-
Glycosylation in health and diseaseNature Reviews. Nephrology 15:346–366.https://doi.org/10.1038/s41581-019-0129-4
-
Comparative protein modelling by satisfaction of spatial restraintsJournal of Molecular Biology 234:779–815.https://doi.org/10.1006/jmbi.1993.1626
-
The UDP-sugar-sensing P2Y(14) receptor promotes Rho-mediated signaling and chemotaxis in human neutrophilsAmerican Journal of Physiology. Cell Physiology 303:C490–C498.https://doi.org/10.1152/ajpcell.00138.2012
-
UDP-glucose promotes neutrophil recruitment in the lungPurinergic Signalling 12:627–635.https://doi.org/10.1007/s11302-016-9524-5
-
GROMACS: fast, flexible, and freeJournal of Computational Chemistry 26:1701–1718.https://doi.org/10.1002/jcc.20291
-
Facile synthesis of sugar nucleotides from common sugars by the cascade conversion strategyJournal of the American Chemical Society 144:9980–9989.https://doi.org/10.1021/jacs.2c03138
-
Facile enzymatic synthesis of ketosesAngewandte Chemie 54:12654–12658.https://doi.org/10.1002/anie.201505714
-
Facile enzymatic synthesis of phosphorylated ketopentosesACS Catalysis 6:1649–1654.https://doi.org/10.1021/acscatal.5b02234
-
CHARMM-GUI Membrane Builder toward realistic biological membrane simulationsJournal of Computational Chemistry 35:1997–2004.https://doi.org/10.1002/jcc.23702
-
Cofactor‐driven cascade reactions enable the efficient preparation of sugar nucleotidesAngewandte Chemie International Edition 61:e202115696.https://doi.org/10.1002/anie.202115696
Article and author information
Author details
Funding
Shanghai Municipal Science and Technology Major Project
- Xi Cheng
Lingang Laboratory grant (LG202102-01-01)
- Xi Cheng
Lingang Laboratory grant (LG-QS-202206-08)
- Liuqing Wen
Fund of Youth Innovation Promotion Association (2022077)
- Xi Cheng
National Key Research and Development Program of China (2021YFA1301900)
- Xi Cheng
National Natural Science Foundation of China (22007092)
- Liuqing Wen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was partially supported by Shanghai Municipal Science and Technology Major Project (to XC and LW), Lingang Laboratory grant (LG202102-01-01 to XC, LG-QS-202206-08 to LW), Fund of Youth Innovation Promotion Association (2022077 to XC), National Key Research and Development Program of China (2021YFA1301900 to XC), and National Natural Science Foundation of China (22007092 to LW).
Copyright
© 2023, Zhao et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 335
- views
-
- 72
- downloads
-
- 2
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
B-cell repertoires are characterized by a diverse set of receptors of distinct specificities generated through two processes of somatic diversification: V(D)J recombination and somatic hypermutations. B cell clonal families stem from the same V(D)J recombination event, but differ in their hypermutations. Clonal families identification is key to understanding B-cell repertoire function, evolution and dynamics. We present HILARy (High-precision Inference of Lineages in Antibody Repertoires), an efficient, fast and precise method to identify clonal families from single- or paired-chain repertoire sequencing datasets. HILARy combines probabilistic models that capture the receptor generation and selection statistics with adapted clustering methods to achieve consistently high inference accuracy. It automatically leverages the phylogenetic signal of shared mutations in difficult repertoire subsets. Exploiting the high sensitivity of the method, we find the statistics of evolutionary properties such as the site frequency spectrum and 𝑑𝑁∕𝑑𝑆 ratio do not depend on the junction length. We also identify a broad range of selection pressures spanning two orders of magnitude.
-
- Cell Biology
- Computational and Systems Biology
Circular RNAs represent a class of endogenous RNAs that regulate gene expression and influence cell biological decisions with implications for the pathogenesis of several diseases. Here, we disclose a novel gene-regulatory role of circHIPK3 by combining analyses of large genomics datasets and mechanistic cell biological follow-up experiments. Using time-course depletion of circHIPK3 and specific candidate RNA-binding proteins, we identify several perturbed genes by RNA sequencing analyses. Expression-coupled motif analyses identify an 11-mer motif within circHIPK3, which also becomes enriched in genes that are downregulated upon circHIPK3 depletion. By mining eCLIP datasets and combined with RNA immunoprecipitation assays, we demonstrate that the 11-mer motif constitutes a strong binding site for IGF2BP2 in bladder cancer cell lines. Our results suggest that circHIPK3 can sequester IGF2BP2 as a competing endogenous RNA (ceRNA), leading to target mRNA stabilization. As an example of a circHIPK3-regulated gene, we focus on the STAT3 mRNA as a specific substrate of IGF2BP2 and validate that manipulation of circHIPK3 regulates IGF2BP2-STAT3 mRNA binding and, thereby, STAT3 mRNA levels. Surprisingly, absolute copy number quantifications demonstrate that IGF2BP2 outnumbers circHIPK3 by orders of magnitude, which is inconsistent with a simple 1:1 ceRNA hypothesis. Instead, we show that circHIPK3 can nucleate multiple copies of IGF2BP2, potentially via phase separation, to produce IGF2BP2 condensates. Our results support a model where a few cellular circHIPK3 molecules can induce IGF2BP2 condensation, thereby regulating key factors for cell proliferation.