Bacteria-Catalyzed Arginine Glycosylation in Pathogens and Host
- 1Institute of Infection and Immunity, Taihe Hospital, Hubei University of Medicine, Shiyan, China
- 2College of Biomedicine and Health, Huazhong Agricultural University, Wuhan, China
- 3College of Life Science and Technology, Huazhong Agricultural University, Wuhan, China
In recent years, protein glycosylation in pathogenic bacteria has attracted more and more attention, and accumulating evidence indicated that this type of posttranslational modification is involved in many physiological processes. The NleB from several enteropathogenic bacteria species as well as SseK from Salmonella enterica are type III secretion system effectors, which have an atypical N-acetylglucosamine (N-GlcNAc) transferase activity that specifically modified a conserved arginine in TRADD, FADD, and RIPK1. NleB/SseKs GlcNAcylation of death domain proteins abrogates homotypic and heterotypic death receptors/adaptors interactions, thereby blocking an important antimicrobial host response. Interestingly, NleB/SseKs could also GlcNAcylate themselves, and self-GlcNAcylation of NleB, SseK1, and SseK3 are crucial for their biological activity during infection. In addition, EarP (EF-P specific arginine rhamnosyl transferase for Posttranslational activation) catalyzes arginine rhamnosylation of translation elongation factor P (EF-P). Importantly, this kind of N-linked protein glycosylation is not only important for EF-P dependent rescue of polyproline stalled ribosomes but also for pathogenicity in Pseudomonas aeruginosa and other clinically relevant bacteria. Glycosylation of arginine is unique because the guanidine group of arginine has a high acid dissociation constant value and representing an extremely poor nucleophile. Recently, the crystal structures of NleB, SseKs, EarP, arginine GlcNAcylated death domain-containing proteins, NleB/FADD-DD, and EarP/EF-P/dTDP-β-L-rhamnose were solved by our group and other groups, revealing the unique catalytic mechanisms. In this review, we provide detailed information about the currently known arginine glycosyltransferases and their potential catalytic mechanisms.
Glycosylation is one of the most abundant and complex posttranslational modifications of proteins and involved in diverse processes such as cell differentiation and growth, signaling cascades, tumorigenesis, as well as host-pathogen interactions (Haltiwanger and Lowe, 2004; Marth and Grewal, 2008; Lu et al., 2015; Stowell et al., 2015). Additionally, glycosylation affects protein properties including stability, folding, and solubility (Moremen et al., 2012). Glycosyltransferases (GTs) are enzymes that establish natural glycosidic linkages, which catalyze the transfer of a sugar moiety from a glycosyl donor to a substrate, and are classified on the basis of their structure as GT-A, GT-B, GT-C, or GT-D fold (Lairson et al., 2008; Gloster, 2014; Zhang et al., 2014; Liang et al., 2015; Park et al., 2018). The GT-A type GT is composed of two β/α/β Rossmann domains, and the two domains are closely interlinked. Out of this, GT-A enzymes possess a conserved Asp-X-Asp (X is any amino acid) motif in which the carboxylates coordinate Mn2+/Mg2+. Just as the GT-A type, the structure of GT-B enzymes is made up of two β/α/β Rossmann-like domains as well. However, in contrast with the architecture of GT-A enzymes, the two β/α/β Rossmann domains in GT-B are linked flexibly. In addition, the metal ions are not required in GT-B enzymes either. The GT-C fold enzymes have numbers of transmembrane helices, and they also have a long-loop region within an active site. One study showed that DUF1792 has a Rossmann-like fold, but the sequence and structure are quite different with the currently annotated type GT-A, GT-B, or GT-C folds, which is why it was defined as a GT-D glycosyltransferase fold.
According to the glycosidic linkage, protein glycosylation can also be divided into O-linked or N-linked glycosylation (Jensen et al., 2012; Rojas-Macias et al., 2019). Typically, O-linked glycans are conjugated to the hydroxyl oxygen of serine and threonine (Van Den Steen et al., 1998; Yang and Qian, 2017). In contrast, N-linked glycans are usually linked to the amide side chain of an asparagine residue in an N-X-S/T (X can be any amino acid, but not proline) motif (Shakin-Eshleman et al., 1996; Mellquist et al., 1998). Unlike with the popular and fully understood N-linked glycosylation of asparagine, N-acetylglucosaminylation (N-GlcNAcylation) that occurs on arginine is quite rare, and the special molecular mechanism of arginine N-GlcNAcylation is completely unknown. It has so far only three reported cases of arginine glycosylation. Self-β-arginine glycosylation was discovered first in sweet corn protein amylogenin in 1995 (Singh et al., 1995). However, there have been no follow up reports. Secondly, the NleB from several enteropathogenic bacteria species as well as SseK from Salmonella enterica are type III secretion system (T3SS) effectors, which were shown to inactivate host death receptors/adaptors by an unprecedented N-GlcNAcylation of a conserved arginine (Li et al., 2013; Pearson et al., 2013). NleB/SseK manipulate host death receptor signaling pathways facilitate the pathogens infection and evade host immune defenses. NleB homologs are present in pathogenic Escherichia coli, Citrobacter rodentium (NleBc), and S. enterica (SseK1/2/3) (Deng et al., 2004; Araujo-Garrido et al., 2020). It should be noticed that enteropathogenic and enterohemorrhagic E. coli (EPEC and EHEC) have two copies of NleB, termed NleB1 and NleB2, and share about 61% amino acid sequence homology (Perna et al., 2001; Iguchi et al., 2009). In addition, EPEC NleB1, EPEC NleB2, EHEC NleB1, EHEC NleB2, SseK1, SseK2, and SseK3 is about 89, 60, 89, 60, 57, 53, and 52% identical to C. rodentium NleB, respectively (Araujo-Garrido et al., 2020). Interestingly, when compared with NleB1/NleBc/SseK1/3, NleB2, and SseK2 possessed a much lower GlcNAcylation activity (Li et al., 2013; Pearson et al., 2013; El Qaidi et al., 2017; Gunster et al., 2017; Newson et al., 2019). In the third case, a conserved arginine of the bacterial translation elongation factor P (EF-P) is rhamnosylated by EarP (EF-P specific arginine rhamnosyl transferase for Posttranslational activation) (Lassak et al., 2015; Rajkovic et al., 2015; Yanagisawa et al., 2016). Notably, this unique modification is important for EF-P dependent rescue of polyproline stalled ribosomes in clinically relevant bacteria such as Pseudomonas aeruginosa and Neisseria meningitides (Lassak et al., 2015; Yanagisawa et al., 2016). Moreover, several studies have shown that EF-P and EarP contribute to the pathogenicity of P. aeruginosa and N. meningitidis by controlling the translation of proline stretch-containing proteins critical for modulating motility, antibiotic resistance, and other traits that play key roles in establishing virulence (Lassak et al., 2015; Rajkovic et al., 2015; Yanagisawa et al., 2016).
Here we provide a summary of bacterial arginine glycosyltransferases and their targets in recent research progress, the unique catalytic mechanisms for arginine glycosylation are discussed as well.
Arginine N-Acetylglucosamine Transferase in Pathogenic E. COLI and C. Rodentium
EPEC is an attaching/effacing (A/E) pathogen that usually leads to severe watery diarrhea, which remains a serious health issue in developing countries (Kotloff et al., 2013). A related E. coli pathotype, EHEC, is the predominant pathogen of bloody diarrhea and hemolytic uremic syndrome (HUS) (Nguyen and Sperandio, 2012). These human bacterial pathogens, together with C. rodentium, a natural murine intestinal bacterium that behaved as the related human pathogens EPEC and EHEC, usually translocate a core set of effectors into host cells to antagonize host defense (Gaytan et al., 2016; Pearson et al., 2016; Pinaud et al., 2018; Shenoy et al., 2018).
Pathogenic E. coli and C. rodentium T3SS effectors, such as Tir (Ruchaud-Sparagano et al., 2011), EspL (Pearson et al., 2017), NleB (Nadler et al., 2010; Newton et al., 2010), NleC (Yen et al., 2010; Baruch et al., 2011; Muhlen et al., 2011; Pearson et al., 2011; Shames et al., 2011; Sham et al., 2011), NleD (Baruch et al., 2011; Creuzburg et al., 2017), and NleE (Nadler et al., 2010; Newton et al., 2010; Zhang et al., 2011), all of which could manipulate the host innate immune system, including the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling and death receptor signaling, via several different mechanisms. It should be noticed that NleB is required for virulence of C. rodentium in vivo (Kelly et al., 2006; Wickham et al., 2006). More importantly, several studies suggested that NleB, to some extent, is associated with the prevalence of human EHEC outbreaks and the outcome of infection (Wickham et al., 2006).
In 2010, Nadler et. al and Newton et. al reported that both NleE and NleB could inhibit NF-κB activation (Nadler et al., 2010; Newton et al., 2010). However, the inhibition activity of NleE and NleB is different, NleE could inhibit both TNFα and IL-1β stimulated NF-κB activation, whereas NleB effector could only inhibit the TNF signaling pathway (Newton et al., 2010; Ruchaud-Sparagano et al., 2011). Although it is well-known that NleB plays an important role in the suppresses NF-κB activation, but the underlying mechanisms are poorly understood. In 2013, one study proposed that glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was the target of C. rodentium NleB, and NleB acted as an O-GlcNAc transferase that modified GAPDH (Gao et al., 2013). Glycosylation of GAPDH inhibited the activity of the tumor necrosis factor receptor-associated factors 2 (TRAF2) and TRAF3, thereby leading to reduced NF-κB signaling and type I IFN signaling (Gao et al., 2013, 2016). Just a few months later, two research groups discovered independently that a critical function of NleB in hijacking of the host death receptor signaling and interfering with host defense (Li et al., 2013; Pearson et al., 2013). Surprisingly, NleB has an unusual N-GlcNAcylation activity toward a conserved arginine (Arg235 in TRADD, Arg117 in FADD, and Arg603 in RIPK1) in host death domain-containing proteins (Li et al., 2013; Pearson et al., 2013). NleB GlcNAcylation of these proteins abrogated homotypic and heterotypic death receptors/adaptors interactions, resulting in disrupting TNF signaling in EPEC or C. rodentium infected cells (Figure 1) (Li et al., 2013; Pearson et al., 2013; Ding et al., 2019). NleB could also block Fas ligand and TNF-associated apoptosis-inducing ligand (TRAIL)-induced cell death by preventing assembly of the canonical death inducing signaling complex (DISC) (Figure 1) (Li et al., 2013; Pearson et al., 2013; Ding et al., 2019). It should be noticed that other host death domain-containing adaptors, such as interleukin-1 receptor-associated kinase1 (IRAK1) and myeloid differentiation primary response 88 (MYD88), located downstream of the of IL-1 receptor (IL-1R), and lack the conserved arginine were not GlcNAcylated by NleB (Figure 1) (Li et al., 2013). This finding perfectly explained the previous observation that NleB could selectively inhibit TNF-α but not IL-1β activation of the NF-κB signaling (Newton et al., 2010; Ruchaud-Sparagano et al., 2011). Further, hijacking of host death receptor signaling pathways by NleB was required for C. rodentium colonization in the mouse model (Li et al., 2013; Pearson et al., 2013; Wong Fok Lung et al., 2016). In addition, several studies suggested that GAPDH is another specific target of EHEC NleB1, and EHEC NleB1-mediated GAPDH GlcNAcylation at Arg197 and Arg200 (Gao et al., 2013; El Qaidi et al., 2017, 2018; Park et al., 2018).
Figure 1. Inhibition of NF-κB signaling and death receptor signaling by NleB and SseK1/3. GlcNAcylation of GAPDH by NleB/SseK1 would suppress TRAF2 polyubiquitination and NF-κB activation. GlcNAcylation of death domain (DD) proteins by NleB and SseK1/3 abrogates homotypic and heterotypic death receptors/adaptors interactions and the assembly of TNFR1 complex, leading to disrupting TNF signaling in EPEC or Citrobacter rodentium infected cells, including NF-κB signaling, apoptosis, and necroptosis. NleB also blocked Fas ligand and TNF-associated apoptosis-inducing ligand (TRAIL)-induced cell death by preventing assembly of the canonical death inducing signaling complex (DISC). In contrast, interleukin-1 receptor-associated kinase1 (IRAK1), and myeloid differentiation primary response 88 (MYD88), lacking the conserved arginine were not GlcNAcylated by NleB and SseK1/3. Besides, the site-directed RA mutants of NleB (NleBArg13/53/159/293Ala), SseK1 (SseK1Arg30/158/339Ala), and SseK3 (SseK3Arg153/184/305/335Ala) abolished or attenuated the capability of enzyme activity toward their death domain-containing targets during infection, and loss of self-GlcNAcylation of NleB, SseK1, and SseK3 couldn't inhibit TNFα- or TRAIL-induced cell death.
Hypoxia-inducible factor 1-alpha (HIF-1α), which act as a key regulator of cellular O2 homeostasis, plays an important role in regulating oxidative glucose metabolism and glycolytic gene expression in the glucose metabolism pathway (Iyer et al., 1998; Cheng et al., 2014; Shukla et al., 2017). Remarkably, in more recent work, Wuhan Xiao and co-workers discovered that HIF-1α was GlcNAcylated at a conserved arginine (Arg18) during C. rodentium or EPEC infection, and the modification enhanced HIF-1α transcriptional activity, thus inducing downstream glucose metabolism-associated gene [such as glucose transporter 1 (GLUT1) gene] expression to alter host glucose metabolism (Xu et al., 2018). Additionally, a more recently study has shown that bacterial glutathione synthetase (GshB) was GlcNAcylated by C. rodentium NleB on Arg256 (El Qaidi et al., 2020). Further, NleB-mediated GlcNAcylation of GshB contributed to C. rodentium survival in oxidative stress conditions (El Qaidi et al., 2020).
Arginine N-Acetylglucosamine Transferase in S. Enterica
S. enterica is a motile, non-spore-forming, and intracellular Gram-negative pathogen that causes both localized and systemic diseases in a wide range of mammals (Gal-Mor et al., 2014). Unlike EPEC, which possesses one T3SS, pathogenic serovars of Salmonella possesses two T3SSs, T3SS1 (encoded by Salmonella pathogenicity islands 1, SPI-1), and T3SS2 (encoded by SPI-2), that inject numerous of effectors into host cells to benefit bacterial invasion and survival (Pearson et al., 2016; Pinaud et al., 2018). Interestingly, SPI-1 and SPI-2, which are activated at different infection stages. T3SS1 is activated upon contact with intestinal epithelial cells, and responsible for the invasion (Zhang et al., 2018; Lou et al., 2019). In contrast, the T3SS2 is expressed after Salmonella has entered host cell (Figueira et al., 2013; Jennings et al., 2017). Moreover, T3SS2 effectors manipulate vesicular trafficking, thereby enhancing Salmonella intracellular survival (Jennings et al., 2017). Three closely related Salmonella T3SS effectors, SseK1, SseK2, and SseK3, are translocated by the SPI-2 T3SS, behave as NleB-like arginine glycosyltransferase, although they displayed distinct differences in host substrate specificity (Brown et al., 2011; Li et al., 2013; Pearson et al., 2013; El Qaidi et al., 2017). Data by Günster et. al suggested SseK1 caused the GlcNAcylation of TRADD and FADD, whereas SseK3 resulted in weak GlcNAcylation of TRADD but not FADD (Gunster et al., 2017). However, another study reported that SseK1 only glycosylated GAPDH but not FADD in vitro (El Qaidi et al., 2017). Outside of this, a recent report coupled with our study provided evidence that endogenous SseK1 modified TRADD, but not FADD, and endogenous SseK3 modified TNFR1 (at Arg376) and TRAILR (at Arg293) during Salmonella infection (Newson et al., 2019; Pan et al., 2019). Significantly, in line with the biological activity of NleB effector, SseK1 and SseK3 inhibit NF-κB activity as well as TNFα-induced cell death (Figure 1) (Gunster et al., 2017; Newson et al., 2019). In contrast with SseK1 and SseK3, SseK2 result in inhibition of TNFα-induced NF-κB reporter activation in vitro only (Gunster et al., 2017; Newson et al., 2019). Additionally, SseK2 didn't induce detectable arginine GlcNAcylation during Salmonella (ΔsseK1/3 or ΔsseK1/3+pSseK2) infection (Gunster et al., 2017; Newson et al., 2019). Therefore, SseK2 may have a much weaker enzyme-substrate interaction network when compared with SseK1 and SseK3. Remarkably, TRIM32, an E3 ubiquitin ligase, was the first potential binding partner identified for SseK3 (Yang et al., 2015). However, SseK3 could not glycosylate TRIM32. Besides, TRIM32 was not required for SseK3 to inhibit NF-κB signaling either (Yang et al., 2015).
Interestingly, two studies have noticed that in addition to GlcNAcylation of TRADD, FADD, and RIPK1, NleB/SseKs could also GlcNAcylate themselves when over-expressed, though the functional significance of this modification is almost completely unknown (Park et al., 2018; Newson et al., 2019). Fortunately, in a more recent study, our study revealed that Arg13/53/159/293 in NleB, Arg30/158/339 in SseK1, and Arg153/184/305/335 in SseK3 were the self-GlcNAcylation sites, which is consistent with the results from one previous study (Newson et al., 2019; Pan et al., 2019). Moreover, the site-directed mutants, i.e., NleBArg13/53/159/293Ala (NleB RA mutant), SseK1Arg30/158/339Ala (SseK1 RA mutant), and SseK3Arg153/184/305/335Ala (SseK3 RA mutant), abolished or attenuated the enzyme activity toward their death domain-containing targets (Figure 1) (Pan et al., 2019). Importantly, the NleB RA mutant, the SseK1 RA mutant, and the SseK3 RA mutant could not inhibit TNFα- or TRAIL-induced cell death (Pan et al., 2019).
Arginine Rhamnosyltransferase Earp of P. Aeruginosa and Other Clinically Relevant Bacterial Species
Ribosomes are the workplaces of protein biosynthesis, the process of translating mRNA into protein. Surprisingly, the ribosomes form peptide bonds among different amino acids with various efficiency. Proline is the least efficient one, both as a donor in the peptidyl-tRNA binding site (P-site) and as an acceptor in the tRNA exiting site (E-site) (Muto and Ito, 2008; Pavlov et al., 2009; Johansson et al., 2011). In this case, ribosomes stall when met with an XPP/PPX motif (Elgamal et al., 2014). Ribosome stalling at polyproline motifs is rescued by the eukaryotic and archaeal elongation factor 5A (e/aIF5A) and its bacterial ortholog, the EF-P, but only when e/aIF5A and EF-P is post-translationally activated (Saini et al., 2009; Zou et al., 2012; Gutierrez et al., 2013; Hersch et al., 2013; Peil et al., 2013; Ude et al., 2013). In bacteria such as E. coli, EF-P is activated when the conserved Lys34 (K34 is conserved in the EF-Ps among about 80% of bacteria) is (R)-β-lysinylated and hydroxylated (Zou et al., 2012; Hersch et al., 2013; Peil et al., 2013; Ude et al., 2013). Comparably, in eukaryotes, conserved Lys50 in e/aIF5A is extended to hypusine (Saini et al., 2009; Gutierrez et al., 2013).
Unusually, several recent studies discovered that the EF-P (efp and earP are always directly adjacent) proteins from about 9% of bacteria, such as Shewanella oneidensis, P. aeruginosa, N. meningitidis, N. gonorrhoeae, and Bordetella pertussis, have a conserved Arg at the position corresponding to Lys34 (Lassak et al., 2015; Rajkovic et al., 2015; Yanagisawa et al., 2016). Significantly, Arg32 in EF-P is rhamnosylated by EarP, and this type of posttranslational modification strategy is crucial for EF-P dependent rescue of polyproline stalled ribosomes in P. aeruginosa and N. meningitidis (Figure 2) (Lassak et al., 2015; Rajkovic et al., 2015; Yanagisawa et al., 2016). P. aeruginosa is a Gram-negative, rod-shaped, asporogenous, and monoflagellated bacterium, and is classified as an opportunistic pathogen (De Lorenzo, 2015). It has been reported that lots of P. aeruginosa virulence factors, such as rhamnolipids and pyocyanin, are polyproline-containing proteins, showing a dependence on EF-P for their translation (Lassak et al., 2015). Moreover, proteins involved with motility, protein synthesis, and DNA replication also act as putative EF-P targets in P. aeruginosa (Rajkovic et al., 2015). Therefore, EF-P and EarP are crucial for P. aeruginosa pathogenicity. In total, this is the first example of N-linked glycosylation occurring on arginine in bacteria which is important for its own biological function.
Figure 2. EF-P arginine rhamnosylation and mode of action. A conserved arginine (Arg32) of EF-P is rhamnosylated by EarP using dTDP-β-L-rhamnose as a sugar donor, and this type of posttranslational modification is crucial for EF-P dependent rescue of polyproline stalled ribosomes. EF-P, translation elongation factor P; EarP, EF-P specific arginine rhamnosyl transferase for posttranslational activation.
The Catalytic Mechanism for NleB/SseKs-Mediated Arginine GlcNAcylation
Several previous mutagenesis analyses showed that NleB, SseK1, SseK2, and SseK3 contain a putative catalytic DXD motif, and an exchange of these residues to alanine completely inhibited their glycosyltransferase activity (Li et al., 2013; Pearson et al., 2013; Wong Fok Lung et al., 2016). Moreover, present structural data indicated that all of the NleB, SseK1, SseK2, and SseK3 proteins share the GT-A fold with the conserved DXD motif (Esposito et al., 2018; Park et al., 2018; Ding et al., 2019; Newson et al., 2019). Typical for the GT-A family of the glycosyltransferase, the DXD motif in NleB and SseKs is involved in a divalent cation (Mn2+) coordination. Therefore, the NleB and SseK family proteins are GT-A type glycosyltransferase (Esposito et al., 2018; Park et al., 2018; Ding et al., 2019; Newson et al., 2019).
Glycosyltransferases can be divided into “retaining” or “inverting” enzymes based on the stereochemistry of the glycosyl donor's anomeric bond during the glycosylation (Figure 3) (Lairson et al., 2008). Recently, the structures of GlcNAcylated TRADD-DD and RIPK1-DD were solved by one of our groups, revealing a β-configuration glycosidic linkage (Ding et al., 2019). Consistently, we previously reported the first synthesis of arginine GlcNAcylated peptides with a β-glycosidic linkage. Using these glycopeptides, we produced a monoclonal antibody that can specifically recognize GlcNAcylated TRADD-DD, FADD-DD, and RIPK1-DD (Pan et al., 2014). As the α-anomeric linkage in UDP-GlcNAc, these data strongly suggest that NleB is an inverting glycosyltransferase in catalyzing arginine GlcNAcylation (Figure 3). However, two retaining mechanisms have been proposed in the nuclear magnetic resonance (NMR) study of SseK1 and SseK3-catalyzed GlcNAc transfer (Esposito et al., 2018; Park et al., 2018). The major reasons for the discrepancy may due to the hydrolase activity of SseKs (Esposito et al., 2018). In the hydrolysis reaction, SseKs were function as hydrolases, a water molecule rather than the arginine acceptor executes the nucleophilic attack, thus generating an α-GlcNAc group (Esposito et al., 2018).
Figure 3. NleB/SseKs catalyze glycosyl group transfer with either inversion or retention of the anomeric stereochemistry. Glycosyltransferases (NleB/SseKs) can be segregated into “retaining” or “inverting” enzymes according to whether the stereochemistry of the glycosyl donor's anomeric bond is retained or inverted during the transfer.
Glycosylation of arginine is unique because the guanidine group of arginine has a high acid dissociation constant value and is intrinsically an extremely poor nucleophile at physiological pH. The crystal structure of the NleB/FADD-DD complex showed that His182, His281, Tyr283, Tyr284, Trp329, and the negatively charged Glu253 might promote deprotonation of the arginine and therefore facilitate catalysis (Ding et al., 2019). Furthermore, the geometry of the UDP-GlcNAc and the position of the Arg117 in FADD-DD suggest NleB adopts a direct-displacement SN2 (substitution, nucleophilic, bimolecular)-like reaction to transfer the GlcNAc from the UDP-GlcNAc to the Arg117, and the Glu253 is proposed to act as the catalytic base (Figure 4) (Ding et al., 2019). Different from this, some other studies showed that SseK1, SseK2, and SseK3 may be retaining glycosyltransferases and adopt an SNi or orthogonal mechanism for the arginine GlcNAcylation reaction (Esposito et al., 2018; Park et al., 2018).
Figure 4. The potential catalytic mechanism (SN2-like) for NleB-catalyzed arginine GlcNAcylation. Glu253 of NleB acts as the base to deprotonate the guanidinium in TRADD (at Arg235), FADD (at Arg117), and RIPK1 (at Arg603). The arginine then nucleophilically attacks the C1 atom of UDP-GlcNAc, forming an oxocarbenium ion-like transition state that progresses to SN2 (substitution, nucleophilic, bimolecular)-like displacement of the UDP.
The Catalytic Mechanism for EarP-Mediated Arginine Rhamnosylation
Compared with the NleB family, the crystal structure of EarP, revealed a GT-B fold and acted as an inverting arginine rhamnosyltransferase (Krafczyk et al., 2017; Sengoku et al., 2018). In the arginine rhamnosylation reaction, the stereochemistry of the dTDP-β-L-rhamnose's anomeric bond is reversed, resulting in α-rhamnosyl on EF-P (Krafczyk et al., 2017; Wang et al., 2017; Sengoku et al., 2018). Upon successful inverting glycosyl transfer from dTDP-β-L-rhamnose to Arg32 in EF-P, the X-ray crystal structure of EarP/EF-P/dTDP-β-L-rhamnose complex together with the NMR data indicated that EarP probably performed an SN2 displacement reaction, with Asp20 as the general base (Figure 5) (Sengoku et al., 2018). Additionally, EarP-mediated arginine rhamnosylation requires the rhamnose ring of the dTDP-β-L-rhamnose to undergo a suitable structural change to expose the β-anomeric face of the rhamnose C1 atom (Sengoku et al., 2018).
Figure 5. The potential catalytic mechanism (SN2-like) for EarP-catalyzed arginine GlcNAcylation. Asp20 of EarP acts as the base to deprotonate the guanidinium in EFP (at Arg32). The arginine then nucleophilically attacks the C1 atom of TDP-rhamnose, forming an oxocarbenium ion-like transition state that progresses to SN2-like displacement of the TDP.
Arginine glycosylation, the attachment of sugar moieties (GlcNAc or rhamnose) to protein's arginine residue, is a novel type of posttranslational modification. There are three types of bacterial arginine glycosyltransferases that make this kind of modification known to date, including NleB homologs in pathogenic E. coli and C. rodentium, SseK homologs in Salmonella, and EarP family in P. aeruginosa and other relevant bacterial species. Interestingly, NleB/SseKs could also GlcNAcylate themselves, and self-GlcNAcylation of NleB, SseK1, and SseK3 is crucial for their biological activity during infection. Excitingly, the success in generating anti-ArgGlcNAc and anti-ArgRha antibodies provides a powerful tool toward the discovery of novel arginine GlcNAcylated or rhamnosylated proteins, and this is the first step toward a comprehensive understanding of arginine glycosylation in nature (Pan et al., 2014; Li et al., 2016; Wang et al., 2017). Besides, a recent report provided evidence that 100066N and 102644N were acted as arginine-GlcNAc glycosyltransferase inhibitors (El Qaidi et al., 2018). Thus, these two compounds may have utility as reagents to further study arginine GlcNAcylation. In addition, the crystal structures of NleB, SseKs, EarP, arginine GlcNAcylated death domain-containing proteins, NleB/FADD-DD, and EarP/EF-P/dTDP-β-L-rhamnose provide compelling evidence of the catalytic mechanism for arginine glycosylation, which will be advantageous to us for designing NleB, SseKs, and EarP inhibitors against the certain pathogens.
XP, JL, and SL conceptualized the study. XP and JL prepared the original draft. All authors listed made a substantial intellectual contribution to the manuscript and edited and approved it for publication.
This work was supported by National Key Research and Development Programs of China 2018YFA0508000, the Foundation for Innovative Research Team of Hubei Provincial Department of Education T201713, Fundamental Research Funds for the Central Universities (2662017PY011, 2662018PY028, and 2662018JC001), and Huazhong Agricultural University Scientific & Technological Self-Innovation Foundation 2017RC003 to SL.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors apologize to colleagues whose work could not be cited owing to the focus of the review.
Araujo-Garrido, J. L., Bernal-Bayard, J., and Ramos-Morales, F. (2020). Type III secretion effectors with arginine N-glycosyltransferase activity. Microorganisms 8:357. doi: 10.3390/microorganisms8030357
Baruch, K., Gur-Arie, L., Nadler, C., Koby, S., Yerushalmi, G., Ben-Neriah, Y., et al. (2011). Metalloprotease type III effectors that specifically cleave JNK and NF-kappaB. Embo J. 30, 221–231. doi: 10.1038/emboj.2010.297
Brown, N. F., Coombes, B. K., Bishop, J. L., Wickham, M. E., Lowden, M. J., Gal-Mor, O., et al. (2011). Salmonella phage ST64B encodes a member of the SseK/NleB effector family. PLoS ONE 6:e17824. doi: 10.1371/journal.pone.0017824
Cheng, S. C., Quintin, J., Cramer, R. A., Shepardson, K. M., Saeed, S., Kumar, V., et al. (2014). mTOR- and HIF-1α-mediated aerobic glycolysis as metabolic basis for trained immunity. Science 345:1250684. doi: 10.1126/science.1250684
Creuzburg, K., Giogha, C., Wong Fok Lung, T., Scott, N. E., Muhlen, S., Hartland, E. L., et al. (2017). The type III effector NleD from enteropathogenic Escherichia coli differentiates between host substrates p38 and JNK. Infect. Immun. 85:e00620-16. doi: 10.1128/IAI.00620-16
Deng, W., Puente, J. L., Gruenheid, S., Li, Y., Vallance, B. A., Vazquez, A., et al. (2004). Dissecting virulence: systematic and functional analyses of a pathogenicity island. Proc. Natl. Acad. Sci. U.S.A. 101, 3597–3602. doi: 10.1073/pnas.0400326101
Ding, J., Pan, X., Du, L., Yao, Q., Xue, J., Yao, H., et al. (2019). Structural and functional insights into host death domains inactivation by the bacterial arginine GlcNAcyltransferase effector. Mol. Cell 74, 922–935.e926. doi: 10.1016/j.molcel.2019.03.028
El Qaidi, S., Chen, K., Halim, A., Siukstaite, L., Rueter, C., Hurtado-Guerrero, R., et al. (2017). NleB/SseK effectors from Citrobacter rodentium, Escherichia coli, and Salmonella enterica display distinct differences in host substrate specificity. J. Biol. Chem.292, 11423–11430. doi: 10.1074/jbc.M117.790675
El Qaidi, S., Scott, N. E., Hays, M. P., Geisbrecht, B. V., Watkins, S., and Hardwidge, P. R. (2020). An intra-bacterial activity for a T3SS effector. Sci. Rep. 10:1073. doi: 10.1038/s41598-020-58062-y
El Qaidi, S., Zhu, C., Mcdonald, P., Roy, A., Maity, P. K., Rane, D., et al. (2018). High-throughput screening for bacterial glycosyltransferase inhibitors. Front. Cell Infect. Microbiol. 8:435. doi: 10.3389/fcimb.2018.00435
Elgamal, S., Katz, A., Hersch, S. J., Newsom, D., White, P., Navarre, W. W., et al. (2014). EF-P dependent pauses integrate proximal and distal signals during translation. PLoS Genet. 10:e1004553. doi: 10.1371/journal.pgen.1004553
Esposito, D., Gunster, R. A., Martino, L., El Omari, K., Wagner, A., Thurston, T. L. M., et al. (2018). Structural basis for the glycosyltransferase activity of the salmonella effector SseK3. J. Biol. Chem. 293, 5064–5078. doi: 10.1074/jbc.RA118.001796
Figueira, R., Watson, K. G., Holden, D. W., and Helaine, S. (2013). Identification of salmonella pathogenicity island-2 type III secretion system effectors involved in intramacrophage replication of S. enterica serovar typhimurium: implications for rational vaccine design. mBio 4:e00065. doi: 10.1128/mBio.00065-13
Gal-Mor, O., Boyle, E. C., and Grassl, G. A. (2014). Same species, different diseases: how and why typhoidal and non-typhoidal Salmonella enterica serovars differ. Front. Microbiol. 5:391. doi: 10.3389/fmicb.2014.00391
Gao, X., Pham, T. H., Feuerbacher, L. A., Chen, K., Hays, M. P., Singh, G., et al. (2016). Citrobacter rodentium NleB protein inhibits Tumor Necrosis Factor (TNF) receptor-associated factor 3 (TRAF3) ubiquitination to reduce host type I interferon production. J. Biol. Chem. 291, 18232–18238. doi: 10.1074/jbc.M116.738278
Gao, X., Wang, X., Pham, T. H., Feuerbacher, L. A., Lubos, M. L., Huang, M., et al. (2013). NleB, a bacterial effector with glycosyltransferase activity, targets GAPDH function to inhibit NF-κB activation. Cell Host Microbe 13, 87–99. doi: 10.1016/j.chom.2012.11.010
Gaytan, M. O., Martinez-Santos, V. I., Soto, E., and Gonzalez-Pedrajo, B. (2016). Type three secretion system in attaching and effacing pathogens. Front Cell Infect. Microbiol. 6:129. doi: 10.3389/fcimb.2016.00129
Gunster, R. A., Matthews, S. A., Holden, D. W., and Thurston, T. L. M. (2017). SseK1 and SseK3 type III secretion system effectors inhibit NF-κB signaling and necroptotic cell death in salmonella-infected macrophages. Infect. Immun. 85:e00010-17. doi: 10.1128/IAI.00242-17
Gutierrez, E., Shin, B. S., Woolstenhulme, C. J., Kim, J. R., Saini, P., Buskirk, A. R., et al. (2013). eIF5A promotes translation of polyproline motifs. Mol. Cell 51, 35–45. doi: 10.1016/j.molcel.2013.04.021
Hersch, S. J., Wang, M., Zou, S. B., Moon, K. M., Foster, L. J., Ibba, M., et al. (2013). Divergent protein motifs direct elongation factor P-mediated translational regulation in Salmonella enterica and Escherichia coli. mBio 4, e00180–e00113. doi: 10.1128/mBio.00180-13
Iguchi, A., Thomson, N. R., Ogura, Y., Saunders, D., Ooka, T., Henderson, I. R., et al. (2009). Complete genome sequence and comparative genome analysis of enteropathogenic Escherichia coli O127:H6 strain E2348/69. J. Bacteriol. 191, 347–354. doi: 10.1128/JB.01238-08
Iyer, N. V., Kotch, L. E., Agani, F., Leung, S. W., Laughner, E., Wenger, R. H., et al. (1998). Cellular and developmental control of O2 homeostasis by hypoxia-inducible factor 1 alpha. Genes Dev. 12, 149–162. doi: 10.1101/gad.12.2.149
Jennings, E., Thurston, T. L. M., and Holden, D. W. (2017). Salmonella SPI-2 type III secretion system effectors: molecular mechanisms and physiological consequences. Cell Host Microbe. 22, 217–231. doi: 10.1016/j.chom.2017.07.009
Johansson, M., Ieong, K. W., Trobro, S., Strazewski, P., Aqvist, J., Pavlov, M. Y., et al. (2011). pH-sensitivity of the ribosomal peptidyl transfer reaction dependent on the identity of the A-site aminoacyl-tRNA. Proc. Natl. Acad. Sci. U.S.A. 108, 79–84. doi: 10.1073/pnas.1012612107
Kelly, M., Hart, E., Mundy, R., Marches, O., Wiles, S., Badea, L., et al. (2006). Essential role of the type III secretion system effector NleB in colonization of mice by Citrobacter rodentium. Infect. Immun. 74, 2328–2337. doi: 10.1128/IAI.74.4.2328-2337.2006
Kotloff, K. L., Nataro, J. P., Blackwelder, W. C., Nasrin, D., Farag, T. H., Panchalingam, S., et al. (2013). Burden and aetiology of diarrhoeal disease in infants and young children in developing countries (the Global Enteric Multicenter Study, GEMS): a prospective, case-control study. Lancet 382, 209–222. doi: 10.1016/S0140-6736(13)60844-2
Krafczyk, R., Macosek, J., Jagtap, P. K. A., Gast, D., Wunder, S., Mitra, P., et al. (2017). Structural basis for EarP-mediated arginine glycosylation of translation elongation factor EF-P. mBio 8:e01412-17. doi: 10.1128/mBio.01412-17
Lairson, L. L., Henrissat, B., Davies, G. J., and Withers, S. G. (2008). Glycosyltransferases: structures, functions, and mechanisms. Annu. Rev. Biochem. 77, 521–555. doi: 10.1146/annurev.biochem.76.061005.092322
Lassak, J., Keilhauer, E. C., Furst, M., Wuichet, K., Godeke, J., Starosta, A. L., et al. (2015). Arginine-rhamnosylation as new strategy to activate translation elongation factor P. Nat. Chem. Biol. 11, 266–270. doi: 10.1038/nchembio.1751
Li, S., Zhang, L., Yao, Q., Li, L., Dong, N., Rong, J., et al. (2013). Pathogen blocks host death receptor signalling by arginine GlcNAcylation of death domains. Nature 501, 242–246. doi: 10.1038/nature12436
Li, X., Krafczyk, R., Macosek, J., Li, Y. L., Zou, Y., Simon, B., et al. (2016). Resolving the alpha-glycosidic linkage of arginine-rhamnosylated translation elongation factor P triggers generation of the first Arg(Rha) specific antibody. Chem. Sci. 7, 6995–7001. doi: 10.1039/C6SC02889F
Liang, D. M., Liu, J. H., Wu, H., Wang, B. B., Zhu, H. J., and Qiao, J. J. (2015). Glycosyltransferases: mechanisms and applications in natural product development. Chem. Soc. Rev. 44, 8350–8374. doi: 10.1039/C5CS00600G
Mellquist, J. L., Kasturi, L., Spitalnik, S. L., and Shakin-Eshleman, S. H. (1998). The amino acid following an asn-X-Ser/Thr sequon is an important determinant of N-linked core glycosylation efficiency. Biochemistry 37, 6833–6837. doi: 10.1021/bi972217k
Muhlen, S., Ruchaud-Sparagano, M. H., and Kenny, B. (2011). Proteasome-independent degradation of canonical NFkappaB complex components by the NleC protein of pathogenic Escherichia coli. J. Biol. Chem. 286, 5100–5107. doi: 10.1074/jbc.M110.172254
Nadler, C., Baruch, K., Kobi, S., Mills, E., Haviv, G., Farago, M., et al. (2010). The type III secretion effector NleE inhibits NF-kappaB activation. PLoS Pathog. 6:e1000743. doi: 10.1371/journal.ppat.1000743
Newson, J. P. M., Scott, N. E., Yeuk Wah Chung, I., Wong Fok Lung, T., Giogha, C., Gan, J., et al. (2019). Salmonella effectors SseK1 and SseK3 target death domain proteins in the TNF and TRAIL signaling pathways. Mol. Cell Proteomics 18, 1138–1156. doi: 10.1074/mcp.RA118.001093
Newton, H. J., Pearson, J. S., Badea, L., Kelly, M., Lucas, M., Holloway, G., et al. (2010). The type III effectors NleE and NleB from enteropathogenic E. coli and OspZ from shigella block nuclear translocation of NF-kappaB p65. PLoS Pathog. 6:e1000898. doi: 10.1371/journal.ppat.1000898
Pan, M., Li, S., Li, X., Shao, F., Liu, L., and Hu, H. G. (2014). Synthesis of and specific antibody generation for glycopeptides with arginine N-GlcNAcylation. Angew. Chem. Int. Ed. Engl. 53, 14517–14521. doi: 10.1002/anie.201407824
Pan, X., Xue, J., Du, L., Zhuang, X., Cai, X., and Li, S. (2019). Arginine-GlcNAcylation of death domain and NleB/SseK proteins is crucial for bacteria pathogenesis by regulating host cell death. bioRxiv 746883. doi: 10.1101/746883
Park, J. B., Kim, Y. H., Yoo, Y., Kim, J., Jun, S. H., Cho, J. W., et al. (2018). Structural basis for arginine glycosylation of host substrates by bacterial effector proteins. Nat. Commun. 9:4283. doi: 10.1038/s41467-018-06680-6
Pavlov, M. Y., Watts, R. E., Tan, Z., Cornish, V. W., Ehrenberg, M., and Forster, A. C. (2009). Slow peptide bond formation by proline and other N-alkylamino acids in translation. Proc. Natl. Acad. Sci. U.S.A. 106, 50–54. doi: 10.1073/pnas.0809211106
Pearson, J. S., Giogha, C., Muhlen, S., Nachbur, U., Pham, C. L., Zhang, Y., et al. (2017). EspL is a bacterial cysteine protease effector that cleaves RHIM proteins to block necroptosis and inflammation. Nat. Microbiol. 2:16258. doi: 10.1038/nmicrobiol.2016.258
Pearson, J. S., Giogha, C., Ong, S. Y., Kennedy, C. L., Kelly, M., Robinson, K. S., et al. (2013). A type III effector antagonizes death receptor signalling during bacterial gut infection. Nature 501, 247–251. doi: 10.1038/nature12524
Pearson, J. S., Giogha, C., Wong Fok Lung, T., and Hartland, E. L. (2016). The Genetics of enteropathogenic Escherichia coli virulence. Annu. Rev. Genet. 50, 493–513. doi: 10.1146/annurev-genet-120215-035138
Pearson, J. S., Riedmaier, P., Marches, O., Frankel, G., and Hartland, E. L. (2011). A type III effector protease NleC from enteropathogenic Escherichia coli targets NF-κB for degradation. Mol. Microbiol. 80, 219–230. doi: 10.1111/j.1365-2958.2011.07568.x
Peil, L., Starosta, A. L., Lassak, J., Atkinson, G. C., Virumae, K., Spitzer, M., et al. (2013). Distinct XPPX sequence motifs induce ribosome stalling, which is rescued by the translation elongation factor EF-P. Proc. Natl. Acad. Sci. U.S.A. 110, 15265–15270. doi: 10.1073/pnas.1310642110
Perna, N. T., Plunkett, G. III., Burland, V., Mau, B., Glasner, J. D., Rose, D. J., Mayhew, G. F., et al. (2001). Genome sequence of enterohaemorrhagic Escherichia coli O157:H7. Nature 409, 529–533. doi: 10.1038/35054089
Rajkovic, A., Erickson, S., Witzky, A., Branson, O. E., Seo, J., Gafken, P. R., et al. (2015). Cyclic rhamnosylated elongation factor p establishes antibiotic resistance in pseudomonas aeruginosa. mBio 6:e00823. doi: 10.1128/mBio.00823-15
Rojas-Macias, M. A., Mariethoz, J., Andersson, P., Jin, C., Venkatakrishnan, V., Aoki, N. P., et al. (2019). Towards a standardized bioinformatics infrastructure for N- and O-glycomics. Nat. Commun. 10:3275. doi: 10.1038/s41467-019-11131-x
Ruchaud-Sparagano, M. H., Muhlen, S., Dean, P., and Kenny, B. (2011). The enteropathogenic E. coli (EPEC) tir effector inhibits NF-κB activity by targeting TNF α receptor-associated factors. PLoS Pathog. 7:e1002414. doi: 10.1371/journal.ppat.1002414
Sengoku, T., Suzuki, T., Dohmae, N., Watanabe, C., Honma, T., Hikida, Y., et al. (2018). Structural basis of protein arginine rhamnosylation by glycosyltransferase EarP. Nat. Chem. Biol. 14, 368–374. doi: 10.1038/s41589-018-0002-y
Shakin-Eshleman, S. H., Spitalnik, S. L., and Kasturi, L. (1996). The amino acid at the X position of an Asn-X-Ser sequon is an important determinant of N-linked core-glycosylation efficiency. J. Biol. Chem. 271, 6363–6366. doi: 10.1074/jbc.271.11.6363
Sham, H. P., Shames, S. R., Croxen, M. A., Ma, C., Chan, J. M., Khan, M. A., et al. (2011). Attaching and effacing bacterial effector NleC suppresses epithelial inflammatory responses by inhibiting NF-κB and p38 mitogen-activated protein kinase activation. Infect. Immun. 79, 3552–3562. doi: 10.1128/IAI.05033-11
Shames, S. R., Bhavsar, A. P., Croxen, M. A., Law, R. J., Mak, S. H., Deng, W., et al. (2011). The pathogenic Escherichia coli type III secreted protease NleC degrades the host acetyltransferase p300. Cell Microbiol. 13, 1542–1557. doi: 10.1111/j.1462-5822.2011.01640.x
Shukla, S. K., Purohit, V., Mehla, K., Gunda, V., Chaika, N. V., Vernucci, E., et al. (2017). MUC1 and HIF-1alpha signaling crosstalk induces anabolic glucose metabolism to impart gemcitabine resistance to pancreatic cancer. Cancer Cell 32, 71–87.e77. doi: 10.1016/j.ccell.2017.06.004
Singh, D. G., Lomako, J., Lomako, W. M., Whelan, W. J., Meyer, H. E., Serwe, M., et al. (1995). beta-Glucosylarginine: a new glucose-protein bond in a self-glucosylating protein from sweet corn. FEBS Lett. 376, 61–64. doi: 10.1016/0014-5793(95)01247-6
Ude, S., Lassak, J., Starosta, A. L., Kraxenberger, T., Wilson, D. N., and Jung, K. (2013). Translation elongation factor EF-P alleviates ribosome stalling at polyproline stretches. Science 339, 82–85. doi: 10.1126/science.1228985
Wang, S., Corcilius, L., Sharp, P. P., Rajkovic, A., Ibba, M., Parker, B. L., et al. (2017). Synthesis of rhamnosylated arginine glycopeptides and determination of the glycosidic linkage in bacterial elongation factor P. Chem. Sci. 8, 2296–2302. doi: 10.1039/C6SC03847F
Wickham, M. E., Lupp, C., Mascarenhas, M., Vazquez, A., Coombes, B. K., Brown, N. F., et al. (2006). Bacterial genetic determinants of non-O157 STEC outbreaks and hemolytic-uremic syndrome after infection. J. Infect. Dis. 194, 819–827. doi: 10.1086/506620
Wong Fok Lung, T., Giogha, C., Creuzburg, K., Ong, S. Y., Pollock, G. L., et al. (2016). Mutagenesis and functional analysis of the bacterial arginine glycosyltransferase effector NleB1 from enteropathogenic Escherichia coli. Infect. Immun. 84, 1346–1360. doi: 10.1128/IAI.01523-15
Xu, C., Liu, X., Zha, H., Fan, S., Zhang, D., Li, S., et al. (2018). A pathogen-derived effector modulates host glucose metabolism by arginine GlcNAcylation of HIF-1α protein. PLoS Pathog. 14:e1007259. doi: 10.1371/journal.ppat.1007259
Yanagisawa, T., Takahashi, H., Suzuki, T., Masuda, A., Dohmae, N., and Yokoyama, S. (2016). Neisseria meningitidis translation elongation factor P and its active-site arginine residue are essential for cell viability. PLoS ONE 11:e0147907. doi: 10.1371/journal.pone.0147907
Yang, Z., Soderholm, A., Lung, T. W., Giogha, C., Hill, M. M., Brown, N. F., et al. (2015). SseK3 is a salmonella effector that binds TRIM32 and modulates the host's NF-κB signalling activity. PLoS ONE 10:e0138529. doi: 10.1371/journal.pone.0138529
Yen, H., Ooka, T., Iguchi, A., Hayashi, T., Sugimoto, N., and Tobe, T. (2010). NleC, a type III secretion protease, compromises NF-κB activation by targeting p65/RelA. PLoS Pathog 6:e1001231. doi: 10.1371/journal.ppat.1001231
Zhang, H., Zhu, F., Yang, T., Ding, L., Zhou, M., Li, J., et al. (2014). The highly conserved domain of unknown function 1792 has a distinct glycosyltransferase fold. Nat. Commun. 5:4339. doi: 10.1038/ncomms5339
Zhang, K., Riba, A., Nietschke, M., Torow, N., Repnik, U., Putz, A., et al. (2018). Minimal SPI1-T3SS effector requirement for salmonella enterocyte invasion and intracellular proliferation in vivo. PLoS Pathog. 14:e1006925. doi: 10.1371/journal.ppat.1006925
Keywords: arginine glycosylation, T3SS effectors, death receptor signaling, NleB, SseK, EarP, posttranslational modification, glycosyltransferase
Citation: Pan X, Luo J and Li S (2020) Bacteria-Catalyzed Arginine Glycosylation in Pathogens and Host. Front. Cell. Infect. Microbiol. 10:185. doi: 10.3389/fcimb.2020.00185
Received: 18 February 2020; Accepted: 07 April 2020;
Published: 28 April 2020.
Edited by:Philip R. Hardwidge, Kansas State University, United States
Reviewed by:Sabrina Mühlen, University of Münster, Germany
Christian Rueter, University Hospital Münster, Germany
Copyright © 2020 Pan, Luo and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shan Li, firstname.lastname@example.org
†These authors have contributed equally to this work