Next Article in Journal
Human Cysteine Cathepsins Degrade Immunoglobulin G In Vitro in a Predictable Manner
Next Article in Special Issue
Inhibition of Dephosphorylation of Dolichyl Diphosphate Alters the Synthesis of Dolichol and Hinders Protein N-Glycosylation and Morphological Transitions in Candida albicans
Previous Article in Journal
4-Mercaptobenzoic Acid Labeled Gold-Silver-Alloy-Embedded Silica Nanoparticles as an Internal Standard Containing Nanostructures for Sensitive Quantitative Thiram Detection
Previous Article in Special Issue
Exploring the N-Glycosylation Profile of Glycoprotein B from Human Cytomegalovirus Expressed in CHO and Nicotiana tabacum BY-2 Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Membrane Topological Model of Glycosyltransferases of the GT-C Superfamily

by
Andreia Albuquerque-Wendt
,
Hermann J. Hütte
,
Falk F. R. Buettner
,
Françoise H. Routier
and
Hans Bakker
*
Institute of Clinical Biochemistry, Hannover Medical School, 30625 Hannover, Germany
*
Author to whom correspondence should be addressed.
Current address: Sir William Dunn School of Pathology, University of Oxford, South Parks Road, Oxford OX1 3RE, UK.
Int. J. Mol. Sci. 2019, 20(19), 4842; https://doi.org/10.3390/ijms20194842
Submission received: 20 August 2019 / Revised: 24 September 2019 / Accepted: 25 September 2019 / Published: 29 September 2019

Abstract

:
Glycosyltransferases that use polyisoprenol-linked donor substrates are categorized in the GT-C superfamily. In eukaryotes, they act in the endoplasmic reticulum (ER) lumen and are involved in N-glycosylation, glypiation, O-mannosylation, and C-mannosylation of proteins. We generated a membrane topology model of C-mannosyltransferases (DPY19 family) that concurred perfectly with the 13 transmembrane domains (TMDs) observed in oligosaccharyltransferases (STT3 family) structures. A multiple alignment of family members from diverse organisms highlighted the presence of only a few conserved amino acids between DPY19s and STT3s. Most of these residues were shown to be essential for DPY19 function and are positioned in luminal loops that showed high conservation within the DPY19 family. Multiple alignments of other eukaryotic GT-C families underlined the presence of similar conserved motifs in luminal loops, in all enzymes of the superfamily. Most GT-C enzymes are proposed to have an uneven number of TDMs with 11 (POMT, TMTC, ALG9, ALG12, PIGB, PIGV, and PIGZ) or 13 (DPY19, STT3, and ALG10) membrane-spanning helices. In contrast, PIGM, ALG3, ALG6, and ALG8 have 12 or 14 TMDs and display a C-terminal dilysine ER-retrieval motif oriented towards the cytoplasm. We propose that all members of the GT-C superfamily are evolutionary related enzymes with preserved membrane topology.

Graphical Abstract

1. Introduction

Glycosyltransferases may be classified and named according to the reaction that they catalyze and their substrate specificity. In addition, a classification based on sequence similarities has been adopted and currently organizes glycosyltransferases in 106 distinct CAZy (Carbohydrate Active enZYmes) GT families (http://www.cazy.org/GlycosylTransferases.html) [1]. The majority of these enzymes use nucleotide sugars as donor substrates, whereas others require a polyisoprenoid lipid carrier for presentation of the monosaccharide or oligosaccharide. In eukaryotes, the glycosyltransferases that use dolichol-phosphate (Dol-P)- or dolichol-diphosphate (Dol-PP)-linked substrates act in the endoplasmic reticulum (ER) and are involved in N-glycosylation, glypiation, O-mannosylation, and C-mannosylation of proteins (Figure 1).
Biosynthesis of N-glycans involves the assembly of a Dol-PP-linked precursor that typically contains three glucose, nine mannose, and two N-acetylglucosamine residues. Four of the nine mannoses, as well as the three glucoses, are added on the luminal side of the ER membrane by the mannosyltransferases ALG3, ALG9, and ALG12, and glucosyltransferases ALG6, ALG8, and ALG10. Dolichol-phosphomannose (Dol-P-Man) and dolichol-phosphoglucose (Dol-P-Glc) are the donor substrates for these reactions [2]. The oligosaccharide moiety of Glc3Man9GlcNAc2-PP-Dol is then transferred from the lipid anchor to protein by the oligosaccharyltransferase (OST) complex [3]. STT3A or STT3B, two alternative catalytic subunits of the complex, are themselves glycosyltransferases using a dolichol-linked donor but transfer an oligosaccharide instead of a monosaccharide.
Protein glypiation involves the preassembly of a glycosylphosphatidylinositol (GPI) anchor and subsequent transfer to the protein. The four mannosyltransferases, PIGM, PIGV, PIGB, and PIGZ, catalyze GPI-anchor assembly using Dol-P-Man on the luminal side of the ER membrane [4]. Finally, protein-O-mannosyltransferases [5,6] and C-mannosyltransferases [7,8] use Dol-P-Man to catalyze the transfer of mannose to the hydroxyl group of serine or threonine and to the indole carbon 2 of tryptophan to form O-mannose and C-mannose, respectively. Like N-glycosylation, protein O-mannosylation is found in eukaryotes, bacteria, and archea [9], whereas C-mannosylation is restricted to metazoa and alveolata [8,10]. The human genome encodes six O-mannosyltransferases, classified into two protein families (POMT1 and POMT2 [11]; and TMTC1, TMTC2, TMTC3, and TMTC4 [5]) and four homologs of the C-mannosyltransferase of C. elegans DPY-19 (dumpy-19) [8] (DPY19L1 (DPY-19like1), DPY19L2, DPY19L3, and DPY19L4). Recently, DPY19L1 and DPY19L3 were shown to be active enzymes with different specificities, whereas no activity could yet be demonstrated for DPY19L2 and DPY19L4 [12].
All the above-mentioned enzymes are multitransmembrane proteins with, based on previous data, 7 to 13 predicted or experimentally defined transmembrane domains (TMDs). They have been categorized as enzymes of the GT-C superfamily, primarily because of their multitransmembrane nature and to distinguish them from enzymes with well-defined GT-A and GT-B structural folds [13,14]. Besides the common polyisoprenol phosphate-bound substrates, they possess a catalytic mechanism leading to an inversion of the transferred monosaccharide anomery (inverting glycosyltransferases). However, they exhibit low or nondetectable sequence identity and have thus been assigned to distinct CAZy GT families (Table 1).
Experimental information about the structure and membrane topology of integral membrane glycosyltransferases is still sparse [15]. Previous studies underlined similarity in the overall predicted topology of ALG (asparagine-linked glycosylation) and PIG (phosphatidylinositol glycan anchor) proteins [14,16] and described the presence of a D×D or DD motif, in which aspartic acid (D) can also be replaced by a glutamic acid (E). This motif is located shortly after the first TMD in the predicted first luminal loop [14,17,18] and has been shown to be essential for enzymatic activity in various members of the GT-C superfamily [19,20,21,22,23,24].
Using protein tagging, the catalytic subunit of yeast Stt3 and protein O-mannosyltransferase Pmt1 were suggested to display 11 and seven TMDs, respectively, whereas prediction programs suggested the presence of 13 and 11 TMDs [25,26,27,28]. The structures of bacterial and archaeal OSTs named PglB (protein glycosylation B) and AglB (archaeal glycosylation B) and of the yeast STT3 demonstrated that these enzymes have 13 membrane helices [21,24,29], and a recent crystallization study of the O-mannosyltransferase revealed that this enzyme has 11 TMDs [30]. However, the defined topology of PglB/AglB proteins differs from the 13 TMD models obtained by algorithms, which predicted a TMD in the first external loop (EL1) (localized in the periplasm in bacteria) and failed to detect the fifth transmembrane helix. In contrast, experimental topology mapping of the yeast STT3 correctly predicted the first luminal loop despite its hydrophobicity, but missed two other domains [26].
Previously, we described clear similarities between the OST catalytic subunit and C. elegans C-mannosyltransferase DPY-19 that suggests an evolutionary link between these enzymes [8]. In this work, we generated a membrane topology model of DPY19 family proteins based on TMDs prediction and information provided by the crystal structure of PglB and AglB, and analyzed the conservation of amino acids in the predicted loops. We strengthen this topology model by demonstrating the importance of several conserved luminally oriented amino acid residues for C-mannosylation of proteins. We further propose a consistent membrane topology model for all eukaryotic members of the GT-C superfamily.

2. Results and Discussion

2.1. Membrane Topology Model of DPY19 Proteins

Membrane topology prediction of members of the GT-C family is rather inconsistent. Enzymes of different species often show distinct patterns of TMD probability. To partly overcome this variation, we used the program PolyPhobius (http://phobius.sbc.su.se/) [25]. This algorithm can be supplied with a multiple sequence alignment and predicts TMDs of the first sequence in the alignment (in our case, the human enzyme STT3A or DPY19L1) but bases the prediction on the multiple sequence alignment. At least 25 eukaryotic STT3 and DPY19 homologs were used to generate multiple sequence alignments. These sequences were selected so that the evolutionary most distinct eukaryotic species are represented. TMD probability of human STT3A and DPY19L1 was compared (Figure 2). The topology model of STT3A has been annotated using information provided by the crystal structures of PglB/AglB and yeast STT3 [21,24,29]. Residues 38 to 114 of STT3A were thus interpreted as the first luminal loop despite the high hydrophobicity of part of this region. Likewise, the large hydrophobic region encompassing residues 200 to 230, predicted as one domain by the program, was marked as TMDs 5 and 6. Annotation of the membrane-spanning domains of DPY19L1 was based on comparison with STT3A and corresponds to the 13 domains previously predicted for the C. elegans enzyme [8]. In spite of low sequence identity between STT3 and DPY19 enzymes, the transmembrane prediction profile of DPY19L1 concurred perfectly with the STT3A model (Figure 2). In DPY19 proteins, the hydrophobicity in the first loop is lower than in STT3 and the 5th and 6th membrane domain are distinguishable, showing a perfect agreement between in silico prediction and the annotation based on the homology to OSTs.
In order to identify amino acids residues potentially involved in C-mannosyltransferase activity, we analyzed the conservation of amino acids in DPY19 proteins using the multiple sequence alignment generated for transmembrane prediction. We produced sequence logos (weblogo.berkeley.edu) [31] of the multiple alignments (Figure S1A,B), in which we maintained the numbering of human STT3A and DPY19L1. Comparison of the logo graphs and the annotated TMDs revealed conservation of amino acids in luminal domains (after uneven TMDs), but not in cytoplasmic domains (Figure S1A,B). In Campylobacter lari PglB, four external (or extracytoplasmic) loops (EL) contain catalytic residues D56 (EL1), R147, D154 and D156 (in EL2), E319 (in EL5), and R375 (in EL6) [21]. These active site amino acids are conserved in archaea AglB, with the exception of D156 that is replaced by a histidine [29], and in eukaryotic STT3s. In human STT3A, the amino acids D49, R160, D167, E169, E351, and R405 have been proposed to be catalytic residues [21,24,29].
DPY19 proteins contain an ExE motif that is located shortly after the first TMD and precedes the hydrophobic stretch of EL1, like the ExD motif of STT3s (Figure 2B). A conserved Rx6DNE or Rx6LRE motif was identified in EL2 of eukaryotic STT3s and DPY19s, respectively. As shown below, the Rx6LRE motif is essential for DPY19 function but was proposed to be located on the cytoplasmic site in a recent topology model of DPY19L3 [32]. This is rather unlikely and probably reflects the limits of the experimental approach using luciferase fusion proteins. The 3rd TMD is entirely embedded in the protein in STT3 and might not function as a transmembrane domain without the presence of the following domains.
The large EL5 of PglB/AglB contains a glutamate (E) interacting with both a divalent cation and the asparagine acceptor of these enzymes. In DPY19 proteins, four to six acidic residues are conserved in EL5, but none seems completely conserved (Figure S1B). Finally, as in STT3, the short EL6 of DPY19s exhibits a highly conserved arginine (R) residue followed by a Proline (P) embedded in transmembrane domain 12 (Figure 2B). In the large C-terminal luminal domain, no conservation between STT3 and DPY19 enzymes was observed. This is unsurprising, since this part of the protein has been shown to be involved in protein acceptor substrate binding in STT3 [33].

2.2. Importance of Conserved Amino Acids for C-Mannosylation

We have previously reported the importance of C-mannosylation for stability and transport of the mouse netrin receptor UNC5 to the cell surface [12]. When a soluble C. elegans UNC-5 construct containing the two thrombospondin repeats was expressed in Drosophila Schneider 2 (S2) cells, the non-C-mannosylated protein fragment was secreted and could be purified from the culture medium when the cells were grown at 21 °C (Figure 3A, lower panel, lane 1). In contrast, when the insect cells were grown at 28 °C, UNC-5 could be obtained from the culture medium only if C. elegans C-mannosyltransferase DPY-19 was co-transfected (Figure 3A, upper panel lanes 1 and 2). This system was used to assess the importance of the above-mentioned conserved residues for C. elegans DPY-19 activity. The amino acids residues E63, E65, R211, E220, E400, and R471 of DPY-19 were individually mutated to alanine, and the resulting enzymes were expressed in S2 cells with UNC-5. Protein secretion at 21 and 28 °C was analyzed by Western blot (Figure 3A). At 28 °C, UNC-5 was secreted when expressed in the presence of DPY-19 carrying the mutations E63A, E65A, or E400A, suggesting that these mutations only had a partial or no effect on C-mannosylation. In contrast, the mutations R211A, E220A, and R471A inactivated DPY-19 since no UNC-5 was secreted when transfected cells were grown at 28 °C. To substantiate these results, UNC-5 purified from medium of cells cultivated at 21 °C was digested with trypsin, and the C-mannosylation degree of the resulting peptides was analyzed by liquid chromatography–mass spectrometry (LC-MS). As previously reported [8], an abundant [M+3H]3+ ion at m/z 884.68, corresponding to the tryptic peptide LDGGWSSWSDWSACSSSCHR substituted with two mannoses was observed when UNC-5 was modified by wild-type DPY-19 (Figure 3B). Triple-charged ions at m/z 776.65 and 830.66 of the same peptide with zero and one mannose, respectively, were of much lower intensity. A very similar profile was observed when UNC-5 was purified from cells expressing the mutant protein E65A and E400A, confirming that these mutations had no significant effect on DPY-19 activity. UNC-5 modified by the mutant proteins E63A yielded reduced di-mannosylated peptide, indicating reduced C-mannosylation activity. The conserved E63 and E65 might be redundant, and mutations at most reduce C-mannosylation activity, as seen for several other GT-C enzymes [19,22]. Cells transfected with an empty vector or with a construct encoding the mutant protein R211A, E220A, or R471A showed no sign of C-mannosylation (Figure 3B). These results indicate that the latter residues, conserved between STT3 and DPY19 proteins, are essential for C-mannosylation. Most likely, these amino acids have the same function in STT3 and DPY19 enzymes in the reaction mechanism, but we cannot exclude an influence on expression or localization of the enzyme.

2.3. Membrane Topology Model of Other GT-C Family Members

Since all eukaryotic glycosyltransferases belonging to the GT-C superfamily use Dol-P- or Dol-PP-linked donor substrates, proceed via an inverting mechanism, and require a divalent cation, we assumed similarities between these proteins and considered the possibility to extend the STT3/DPY19 topology model to other GT-C proteins. Previous studies already underlined similarity in the overall topology of the ALG and PIG proteins [16] and created a membrane topology model with 12 TMDs.
Some sequences of ALG and PIG enzymes can be aligned. ALG9, ALG12, PIGB, and PIGZ are grouped together in the GT22 family. ALG3 and PIGM show clear homology with each other as well, although they are placed in different families (GT58 and GT50, respectively). The glucosyltransferases ALG6 and ALG8 (GT57) are very similar but differ from the 3rd glucosyltransferase (ALG 10, GT59). PIGV (GT76) also forms a separate family. Besides the presence of a D×D or similar motif [14,16], the whole ALG and PIG families of enzymes cannot directly be aligned. The same is true for the O-mannosyltransferases, POMT1 and -2 (GT39) and the recently identified TMTC (GT105) family [5].
To identify conserved motifs and putative membrane topology, we selected the 25 to 30 most diverse eukaryotic sequences of each GT-C glycosyltransferase family and used these to generate transmembrane probability graphs and sequence logos, as described above for DPY19. The transmembrane prediction graph of the POMT family (Figure 4) clearly showed common characteristics with those of STT3 and DPY19 (Figure 2). As in these two enzymes, there is a larger luminal loop after the first TMD, followed by an area with TMDs with small loops. Before the last four TMDs, there is again a larger loop. A conserved Rx4D, similar to the motif found in STT3 and DPY19 proteins, was found between predicted TMDs 3 and 4. As in STT3, domains 5 and 6 were considered as two TDMs separated by a short luminal loop, which displays a highly conserved lysine (K). In POMTs, the conserved Rx6P motif situated before the last two transmembrane domains, described above for STT3s and DPY19s, is present as well, however, with eight amino acid distance (Figure 4).
The topology prediction of POMT can perfectly be aligned with STT3 and DPY19, with the conserved motifs at identical positions, although two TMDs less are predicted in the middle part (Figure 4). Based on the conservation in the first luminal loop, conserved catalytic mechanisms between O-mannosyltransferases and other members of the GT-C superfamily were already suggested by Lommel et al. [22]. Using experiments with hybrid proteins, the yeast O-mannosyltransferase Pmt1 was originally proposed to have seven TMDs [28]. The crystal structure of yeast O-mannosyltransferase Pmt1 and Pmt2 that was just reported [30] is, however, in complete agreement with our proposed 11 TMD model for POMT. It also highlights the role of the conserved lysine (K) in EL3 (K204 in POMT1, Figure 4) and R649 in EL5 in Dol-P-Man binding and the structural resemblance with STT3.
TMTC proteins also have a similar TMD prediction (Figure 4), resembling POMT in the number of predicted domains, but with a long C-terminal luminal domain like in STT3 and DPY19. A multiple alignment of diverse members of the TMTC family again reveals conserved areas resembling those of POMT (Figure 4 and Figure S1D).

2.4. GT-C Glycosyltransferases Involved in N-Glycan and GPI Precursor Biosynthesis

The enzymes analyzed so far use proteins as the acceptor substrate. The remaining members of the GT-C superfamily transfer lipid-linked sugars to a lipid-linked acceptor. These enzymes, however, show again a comparable TMD prediction (Figure 5). With few exceptions, four conserved amino acid patterns can be found in the multiple alignments of all enzymes (Figure 6 and Figure S1E–N). The GT22 enzymes and PIGV are predicted to have 11 TMDs (Figure 5), corresponding to those in the O-mannosyltransferases (Figure 4). The conserved patterns in ALG10 align with STT3 and DPY19 (Figure 2), suggesting an identical 13 TMD model. ALG3 and PIGM show 12 TMDs. The 12th domain is, based on the positions of the conserved domains, a C-terminal addition compared to the enzymes with 11 TMDs and displays a dilysine motif [16]. The 12 TMDs model places this motif on the cytoplasmic side of the membrane where it can function as an ER-retrieval signal [34]. ALG6 and ALG8 display three additional TMDs at the C-terminus, thus comprising 14 TMDs, and contain a dilysine motif as well.
The concept that conserved residues should be located on the luminal side of the membrane, where the catalyzed reaction takes place, helped to create topology models of all enzymes of the GT-C superfamily. Four conserved amino acid patterns can be observed in luminal loop 1, 2, 3, and 6(or 5) in all enzymes (Figure 2,Figure 4 and Figure 6). These are not identical among the different GT-C enzymes, but show similar characteristics. The first D×D motif, in which D can be E and x may be absent, is present in all enzymes. The second conserved pattern is more variable, but mostly contains an R and/or D/E/N or Q. The motif of the 3rd luminal loop comprises mostly one conserved positively charged R or K, but is replaced by N or Q in STT3 and DPY19 (Figure 2). ALG9 is an outlier, in which the motifs of the 2nd and 3rd luminal loop are absent. Possibly, these potential catalytic residues are compensated for in other luminal loops.
The most striking conservation is the presence of an R or K in the luminal loop of the cluster of four TMDs, which are at the C-terminus in most enzymes. This conserved positively charged amino acid, predicted to be involved in donor substrate binding [21,30], is always followed by a conserved P that form an Rx4–14P motif. In addition, the large area between TMD 9/10 (or 7/8) has conserved amino acids in all families (Figure S1) that are potentially involved in catalysis, like in STT3. This area is too large and variable to identify patterns common to the whole superfamily, but the presence of conserved amino acids within families is a strong argument for luminal localization. Finally, in most enzymes the C-terminus is oriented toward the ER lumen and also contains conserved amino acids within the families.
The types of conserved amino acids in the GT-C family are those anticipated to be involved in catalysis. To ascertain that these are not the types of pattern always found on luminal sides of membrane proteins, we carried out the same analysis with an ER transmembrane protein that catalyzes a reaction on the cytoplasmic side. The 10 TMD protein dolichyl-phosphate N-acetylglucosaminephosphotransferase (DPAGT1), which has been crystalized [35], indeed shows major conservation in loops on the cytoplasmic side of the membrane (Figure S1O).

3. Materials and Methods

3.1. Bioinformatics

Using the human protein sequences of members of the GT-C family (Figure 1), a PSI-BLAST (Position-Specific Iterative Basic Local Alignment Search Tool) [36] search was performed in the reference protein database, restricted to eukaryotic species. At least 25 most diverse sequences were selected for each protein. For each protein, paralogs (STT3B; DPY19L2, 3, and 4; POMT2; and TMTC2, 3, 4, and various fungal O-mannosyltransferasess) were included. These sequences were used for multiple alignments using the program Praline (ibi.vu.nl/programs/pralinewww) that considers transmembrane prediction in its alignment [37]. These multiple alignments were subsequently used to predict the membrane topology using PolyPhobius (phobius.sbc.su.se/poly.html), which allows to plot membrane topology probability of a protein, but based on a multiple alignment [25], allowing a more reliable prediction.
The same multiple alignments were used to generate a protein sequence logo (weblogo.berkeley.edu/) [31]. For the logos, sequences of the human enzymes (in case of paralogs, the lowest number) were taken as basis, meaning that inserted amino acids present in some other sequences were deleted from the alignment before the logo was generated, to avoid gaps and to maintain the numbering of the human enzyme.

3.2. Plasmid Constructs for S2 Cell Expression

All C. elegans DPY-19 mutants were created using a previously described plasmid: pIB-DPY-19 [8]. DPY-19-E63A, DPY-19-E65A, DPY-19-R211A, DPY-19-E220A, DPY-19-E400A, and DPY-19-R471A were made by amplifying two PCR fragments extended by recognition sequences for the type IIs restriction enzyme BbsI [38] containing the mutated site, which were, after digestion, ligated back in the same vector. The Netrin receptor UNC-5 construct pMT-UNC-5 is described in [8]. As an internal transfection and secretion control, a plasmid expressing EGF repeats 16–20 of Drosophila melanogaster Notch (pMT-EGF), which is described in [39], was used.

3.3. Protein Expression in S2 Cells

Drosophila Schneider (S2) cells were grown in Insect-XPRESS Protein-free Medium (Lonza Basel, Switzerland) and kept at a density between 3 × 106 and 2 × 107 cells/mL at 24 °C under 30 rpm rotation in 75 cm2 nonadherent T-flasks. Transient expression was carried out as described [8] in S2 cells in 24-well plates. An amount of 0.5 µg of pIB-DPY-19 (wild-type or mutants) or empty vector was combined with 0.5 µg of pMT-UNC-5 and 0.5 µg of pMT-EGF. After Six hours, expression was induced with 0.2 mM CuSO4, and cells were subsequently kept at 21 °C under rotation or at 28 °C for 4 days. To produce protein for mass spectrometry (MS) analysis, the protocol was scaled up to 10 mL in 75 cm2 T-flasks.

3.4. Western Blotting

For western blot analysis, 50 µL of cell suspension was centrifuged at 300× g for 5 min to pellet the cells, followed by a second centrifugation of the supernatant at 4000× g. Ten µL of 5× Laemmli with 5% (v/v) β-mercaptoethanol was added to 40 µL of the supernatant and incubated at 95 °C for 10 min. The entire sample was loaded on a 15% (w/v) SDS-PAGE gel. Blotted nitrocellulose membranes were incubated with mouse V5 antibody (Acris Antibodies GmbH, Herford, Germany) and IRDye 800 goat antimouse secondary antibody (LI-COR BioSciences GmbH, Bad Homburg, Germany) and scanned on an LI-COR Odyssey Infrared Scanner.

3.5. Protein Purification

For nickel affinity chromatography of UNC-5, 10 mL of medium were centrifuged as described above and filtered through a 0.2 µm membrane (Waters GmbH, Eschborn, Germany). Binding buffer was added to generate 50 mL of 500 mM NaCl, 20 mM Tris-Hcl (pH8), and 20 mM imidazole (buffer A). The solution was loaded on a 1 mL HisTrap HP column (GE Healthcare Europe GmbH, Freiburg, Germany), washed with 10 mL of buffer, and eluted with a 7 mL gradient from 20 to 500 mM imidazole in buffer A. The peak fractions (A280) were precipitated with 4 volumes of acetone, pooled, and resuspended in 30 µL of 1× Laemmli (without reducing agents). Finally, 1 µL of the preparation was used for western blotting, and the rest of the sample was loaded on a 15% (w/v) SDS-PAGE and stained with Coomassie blue for MS analysis.

3.6. Mass Spectrometry

Protein bands were excised from gels and digested with trypsin at 37 °C in an overnight reaction. Peptides were extracted and analyzed by mass spectrometry, as described [8].

4. Conclusions

In conclusion, we can state that all enzymes of the GT-C superfamily are part of a true homologous family, as was suggested before [14,16,22]. Amino acids involved in STT3 catalysis can mostly be found in other GT-C enzymes. The conservation of charged R/K and acidic D/E amino acid residues suggests common donor substrate recognition and catalysis involving amino acids from different luminal loops. All members of the GT-C family have a similar topological organization (Figure 7). The majority of enzymes present 11 TMDs (POMTs, TMTCs, ALG12, ALG9, PIGB, PIGZ, and PIGV) with the N-terminus oriented toward the cytoplasmic face of the ER membrane and the C-terminus facing the lumen. The first and fourth luminal loops are relatively large and connect the first and four last TMDs to six central TMDs (1-6-4 organization). In the 13 TMDs model adopted by STT3s, DPY19s, and ALG10, two additional central membrane segments are observed, resulting in a 1-8-4 organization. Finally, one (ALG3, PIGM) or three (ALG6 and ALG8) extra C-terminal TMDs give rise to 12 or 14 TMDs (1-6-5 or 1-6-7 organization). In these enzymes, the dilysine-containing C-terminus is oriented toward the cytoplasm, where it can act as retrieval signal.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/20/19/4842/s1, Figure S1: Full multiple alignments presented as logos of all GT-C families.

Author Contributions

A.A.-W., H.J.H., and F.F.R.B. carried out experiments; F.H.R. and H.B. designed experiments and wrote the paper.

Funding

This work was supported by the European Community FP7 People: Marie-Curie Actions (FP7-PEOPLE-2013-ITN-608295) and the Deutsche Forschungsgemeinschaft (DFG FOR2509 BA4091/6-1).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

CAZycarbohydrate active enzymes
Doldolichol
ERendoplasmic reticulum
ELextracytoplasmic loop
Glcglucose
GlcNAcN-acetylglucosamine
GT-Cglycosyltransferase family C
GPIglycosylphosphatidylinositol
Manmannose
OSToligosaccharyltransferase
TMBtransmembrane domain

References

  1. Lombard, V.; Golaconda Ramulu, H.; Drula, E.; Coutinho, P.M.; Henrissat, B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014, 42, D490–D495. [Google Scholar] [CrossRef]
  2. Aebi, M. N-linked protein glycosylation in the ER. Biochim. Biophys. Acta 2013, 1833, 2430–2437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kelleher, D.J.; Gilmore, R. An evolving view of the eukaryotic oligosaccharyltransferase. Glycobiology 2006, 16, 47R–62R. [Google Scholar] [CrossRef] [PubMed]
  4. Kinoshita, T. Biosynthesis and deficiencies of glycosylphosphatidylinositol. Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 2014, 90, 130–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Larsen, I.S.B.; Narimatsu, Y.; Joshi, H.J.; Siukstaite, L.; Harrison, O.J.; Brasch, J.; Goodman, K.M.; Hansen, L.; Shapiro, L.; Honig, B.; et al. Discovery of an O-mannosylation pathway selectively serving cadherins and protocadherins. Proc. Natl. Acad. Sci. USA 2017, 114, 11163–11168. [Google Scholar] [CrossRef] [PubMed]
  6. Loibl, M.; Strahl, S. Protein O-mannosylation: What we have learned from baker’s yeast. Biochim. Biophys. Acta 2013, 1833, 2438–2446. [Google Scholar] [CrossRef]
  7. Doucey, M.A.; Hess, D.; Cacan, R.; Hofsteenge, J. Protein C-mannosylation is enzyme-catalysed and uses dolichyl-phosphate-mannose as a precursor. Mol. Biol. Cell 1998, 9, 291–300. [Google Scholar] [CrossRef]
  8. Buettner, F.F.; Ashikov, A.; Tiemann, B.; Lehle, L.; Bakker, H.C. elegans DPY-19 is a C-mannosyltransferase glycosylating thrombospondin repeats. Mol. Cell 2013, 50, 2953–3002. [Google Scholar] [CrossRef] [PubMed]
  9. Lommel, M.; Strahl, S. Protein O-mannosylation: Conserved from bacteria to humans. Glycobiology 2009, 19, 816–828. [Google Scholar] [CrossRef] [PubMed]
  10. Hoppe, C.M.; Albuquerque-Wendt, A.; Bandini, G.; Leon, D.R.; Shcherbakova, A.; Buettner, F.F.R.; Izquierdo, L.; Costello, C.E.; Bakker, H.; Routier, F.H. Apicomplexan C-mannosyltransferases modify thrombospondin Type I-containing adhesins of the TRAP family. Glycobiology 2018, 28, 333–343. [Google Scholar] [CrossRef]
  11. Manya, H.; Chiba, A.; Yoshida, A.; Wang, X.; Chiba, Y.; Jigami, Y.; Margolis, R.U.; Endo, T. Demonstration of mammalian protein O-mannosyltransferase activity: Coexpression of POMT1 and POMT2 required for enzymatic activity. Proc. Natl. Acad. Sci. USA 2004, 101, 500–505. [Google Scholar] [CrossRef] [PubMed]
  12. Shcherbakova, A.; Tiemann, B.; Buettner, F.F.; Bakker, H. Distinct C-mannosylation of netrin receptor thrombospondin type 1 repeats by mammalian DPY19L1 and DPY19L3. Proc. Natl. Acad. Sci. USA 2017, 114, 2574–2579. [Google Scholar] [CrossRef] [PubMed]
  13. Lairson, L.L.; Henrissat, B.; Davies, G.J.; Withers, S.G. Glycosyltransferases: Structures, functions, and mechanisms. Annu. Rev. Biochem. 2008, 77, 521–555. [Google Scholar] [CrossRef]
  14. Liu, J.; Mushegian, A. Three monophyletic superfamilies account for the majority of the known glycosyltransferases. Protein Sci. 2003, 12, 1418–1431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Gloster, T.M. Advances in understanding glycosyltransferases from a structural perspective. Curr. Opin. Struct. Biol. 2014, 28, 131–141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Oriol, R.; Martinez-Duncker, I.; Chantret, I.; Mollicone, R.; Codogno, P. Common origin and evolution of glycosyltransferases using Dol-P-monosaccharides as donor substrate. Mol. Biol. Evol. 2002, 19, 1451–1463. [Google Scholar] [CrossRef] [PubMed]
  17. Coutinho, P.M.; Deleury, E.; Davies, G.J.; Henrissat, B. An evolving hierarchical family classification for glycosyltransferases. J. Mol. Biol. 2003, 328, 307–317. [Google Scholar] [CrossRef]
  18. Kikuchi, N.; Kwon, Y.D.; Gotoh, M.; Narimatsu, H. Comparison of glycosyltransferase families using the profile hidden Markov model. Biochem. Biophys. Res. Commun. 2003, 310, 574–579. [Google Scholar] [CrossRef]
  19. Hese, K.; Otto, C.; Routier, F.H.; Lehle, L. The yeast oligosaccharyltransferase complex can be replaced by STT3 from Leishmania major. Glycobiology 2009, 19, 160–171. [Google Scholar] [CrossRef] [PubMed]
  20. Kang, J.Y.; Hong, Y.; Ashida, H.; Shishioh, N.; Murakami, Y.; Morita, Y.S.; Maeda, Y.; Kinoshita, T. PIG-V involved in transferring the second mannose in glycosylphosphatidylinositol. J. Biol. Chem. 2005, 280, 9489–9497. [Google Scholar] [CrossRef] [PubMed]
  21. Lizak, C.; Gerber, S.; Numao, S.; Aebi, M.; Locher, K.P. X-ray structure of a bacterial oligosaccharyltransferase. Nature 2011, 474, 350–355. [Google Scholar] [CrossRef] [PubMed]
  22. Lommel, M.; Schott, A.; Jank, T.; Hofmann, V.; Strahl, S. A conserved acidic motif is crucial for enzymatic activity of protein O-mannosyltransferases. J. Biol. Chem. 2011, 286, 39768–39775. [Google Scholar] [CrossRef]
  23. Maeda, Y.; Watanabe, R.; Harris, C.L.; Hong, Y.; Ohishi, K.; Kinoshita, K.; Kinoshita, T. PIG-M transfers the first mannose to glycosylphosphatidylinositol on the lumenal side of the ER. EMBO J. 2001, 20, 250–261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wild, R.; Kowal, J.; Eyring, J.; Ngwa, E.M.; Aebi, M.; Locher, K.P. Structure of the yeast oligosaccharyltransferase complex gives insight into eukaryotic N-glycosylation. Science 2018, 359, 545–550. [Google Scholar] [CrossRef] [PubMed]
  25. Kall, L.; Krogh, A.; Sonnhammer, E.L. An HMM posterior decoder for sequence feature prediction that includes homology information. Bioinformatics 2005, 21, i251–i257. [Google Scholar] [CrossRef] [PubMed]
  26. Kim, H.; von Heijne, G.; Nilsson, I. Membrane topology of the STT3 subunit of the oligosaccharyl transferase complex. J. Biol. Chem. 2005, 280, 20261–20267. [Google Scholar] [CrossRef]
  27. Krogh, A.; Larsson, B.; von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef] [PubMed]
  28. Strahl-Bolsinger, S.; Scheinost, A. Transmembrane topology of pmt1p, a member of an evolutionarily conserved family of protein O-mannosyltransferases. J. Biol. Chem. 1999, 274, 9068–9075. [Google Scholar] [CrossRef]
  29. Matsumoto, S.; Shimada, A.; Nyirenda, J.; Igura, M.; Kawano, Y.; Kohda, D. Crystal structures of an archaeal oligosaccharyltransferase provide insights into the catalytic cycle of N-linked protein glycosylation. Proc. Natl. Acad. Sci. USA 2013, 110, 17868–17873. [Google Scholar] [CrossRef] [Green Version]
  30. Bai, L.; Kovach, A.; You, Q.; Kenny, A.; Li, H. Structure of the eukaryotic protein O-mannosyltransferase Pmt1-Pmt2 complex. Nat. Struct. Mol. Biol. 2019, 26, 704–711. [Google Scholar] [CrossRef]
  31. Crooks, G.E.; Hon, G.; Chandonia, J.M.; Brenner, S.E. WebLogo: A sequence logo generator. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef] [PubMed]
  32. Niwa, Y.; Nakano, Y.; Suzuki, T.; Yamagishi, M.; Otani, K.; Dohmae, N.; Simizu, S. Topological analysis of DPY19L3, a human C-mannosyltransferase. FEBS J. 2018, 285, 1162–1174. [Google Scholar] [CrossRef] [PubMed]
  33. Igura, M.; Maita, N.; Kamishikiryo, J.; Yamada, M.; Obita, T.; Maenaka, K.; Kohda, D. Structure-guided identification of a new catalytic motif of oligosaccharyltransferase. EMBO J. 2008, 27, 234–243. [Google Scholar] [CrossRef] [PubMed]
  34. Jackson, M.R.; Nilsson, T.; Peterson, P.A. Identification of a consensus motif for retention of transmembrane proteins in the endoplasmic reticulum. EMBO J. 1990, 9, 3153–3162. [Google Scholar] [CrossRef] [PubMed]
  35. Yoo, J.; Mashalidis, E.H.; Kuk, A.C.Y.; Yamamoto, K.; Kaeser, B.; Ichikawa, S.; Lee, S.Y. GlcNAc-1-P-transferase-tunicamycin complex structure reveals basis for inhibition of N-glycosylation. Nat. Struct. Mol. Biol. 2018, 25, 217–224. [Google Scholar] [CrossRef] [PubMed]
  36. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [PubMed]
  37. Pirovano, W.; Feenstra, K.A.; Heringa, J. PRALINETM: A strategy for improved multiple alignment of transmembrane proteins. Bioinformatics 2008, 24, 492–497. [Google Scholar] [CrossRef]
  38. Beck, R.; Burtscher, H. Introduction of arbitrary sequences into genes by use of class IIs restriction enzymes. Nucleic Acids Res. 1994, 22, 886–887. [Google Scholar] [CrossRef] [Green Version]
  39. Pandey, A.; Li-Kroeger, D.; Sethi, M.K.; Lee, T.V.; Buettner, F.F.; Bakker, H.; Jafar-Nejad, H. Sensitized genetic backgrounds reveal differential roles for EGF repeat xylosyltransferases in Drosophila Notch signaling. Glycobiology 2018, 28, 849–859. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Human enzymes using dolichol-linked donor substrates. The official protein symbols (www.genames.org) of all known human enzymes are indicated. PIG (phosphatidylinositol glycan anchor) enzymes generate the GPI (glycosylphosphatidylinositol)anchor. ALG (asparagine-linked glycosylation) enzymes transfer four mannoses and three glucoses to the lipid-linked oligosaccharide before it is transferred to proteins by the oligosaccharyltransferase complex. STT3 (staurosporine- and temperature-sensitive) is the catalytic subunit of the oligosaccharyltransferase and occurs in two forms. Two POMTs (protein O-mannosyltransferases) and four TMTC (transmembrane and tetratricopeptide repeat containing) proteins transfer O-mannose to protein. Finally, four homologues of the Caenorhaditis elegans C-mannosyltransferase DPY-19 are present in human.
Figure 1. Human enzymes using dolichol-linked donor substrates. The official protein symbols (www.genames.org) of all known human enzymes are indicated. PIG (phosphatidylinositol glycan anchor) enzymes generate the GPI (glycosylphosphatidylinositol)anchor. ALG (asparagine-linked glycosylation) enzymes transfer four mannoses and three glucoses to the lipid-linked oligosaccharide before it is transferred to proteins by the oligosaccharyltransferase complex. STT3 (staurosporine- and temperature-sensitive) is the catalytic subunit of the oligosaccharyltransferase and occurs in two forms. Two POMTs (protein O-mannosyltransferases) and four TMTC (transmembrane and tetratricopeptide repeat containing) proteins transfer O-mannose to protein. Finally, four homologues of the Caenorhaditis elegans C-mannosyltransferase DPY-19 are present in human.
Ijms 20 04842 g001
Figure 2. (A) Transmembrane probability of STT3A and DPY19L1. Transmembrane domain prediction was carried out using the algorithm PolyPhobius with a multiple alignment comprising 25 most diverse members of the two families, as described in experimental procedures. Proposed transmembrane domain (TMD) numbering is based on the structure of bacterial Campylobacter lari oligosaccharyltransferase [21]. (B) Upper panel: Conserved residues involved in catalysis in STT3. Lower panel: Conserved areas (arrows in Figure 2A) in extracytoplasmic (luminal) loops (EL) of the STT3 and DPY19 family. *E in EL5 is part of the catalytic pocket, but due to the large size of this loop and the presence of multiple D (aspartic acid) and E (glutamic acid) residues, no clear equivalent can be observed in DPY19 enzymes (see supplemental Figure S1A,B). The asparagine (N) in EL3 is highly conserved in STT3, but has not been associated with the catalytic pocket.
Figure 2. (A) Transmembrane probability of STT3A and DPY19L1. Transmembrane domain prediction was carried out using the algorithm PolyPhobius with a multiple alignment comprising 25 most diverse members of the two families, as described in experimental procedures. Proposed transmembrane domain (TMD) numbering is based on the structure of bacterial Campylobacter lari oligosaccharyltransferase [21]. (B) Upper panel: Conserved residues involved in catalysis in STT3. Lower panel: Conserved areas (arrows in Figure 2A) in extracytoplasmic (luminal) loops (EL) of the STT3 and DPY19 family. *E in EL5 is part of the catalytic pocket, but due to the large size of this loop and the presence of multiple D (aspartic acid) and E (glutamic acid) residues, no clear equivalent can be observed in DPY19 enzymes (see supplemental Figure S1A,B). The asparagine (N) in EL3 is highly conserved in STT3, but has not been associated with the catalytic pocket.
Ijms 20 04842 g002
Figure 3. Secretion and C-mannosylation of an UNC-5 fragment modified by native and mutated C. elegans DPY-19. (A) Secretion of UNC-5 from S2 cells grown at 21 and 28 °C. A V5-His6-tagged UNC-5 construct encompassing the two TSRs was transfected in S2 cells without (−) or with (+) C. elegans DPY-19 or with C. elegans DPY-19 carrying the mutation E63A, E65A, R211A, E220A, E400A, or R471A. A plasmid expressing V5-tagged epidermial growth factor (EGF)-repeats 16–20 (EGF) of Drosophila melanogaster Notch was used as a transfection and secretion control. Transfected cells were grown at 21 or 28 °C, and secreted proteins were detected on Western blot using an anti-V5 antibody. (B) Extracted ion chromatograms of the unmodified, mono-C-mannosylated, and di-C-mannosylated TSR2 tryptic peptide LDGGWSSWSDWSACSSSCHR. All graphs are shown with equal intensities (13,000 counts per second).
Figure 3. Secretion and C-mannosylation of an UNC-5 fragment modified by native and mutated C. elegans DPY-19. (A) Secretion of UNC-5 from S2 cells grown at 21 and 28 °C. A V5-His6-tagged UNC-5 construct encompassing the two TSRs was transfected in S2 cells without (−) or with (+) C. elegans DPY-19 or with C. elegans DPY-19 carrying the mutation E63A, E65A, R211A, E220A, E400A, or R471A. A plasmid expressing V5-tagged epidermial growth factor (EGF)-repeats 16–20 (EGF) of Drosophila melanogaster Notch was used as a transfection and secretion control. Transfected cells were grown at 21 or 28 °C, and secreted proteins were detected on Western blot using an anti-V5 antibody. (B) Extracted ion chromatograms of the unmodified, mono-C-mannosylated, and di-C-mannosylated TSR2 tryptic peptide LDGGWSSWSDWSACSSSCHR. All graphs are shown with equal intensities (13,000 counts per second).
Ijms 20 04842 g003
Figure 4. Predicted membrane topology and amino acid conservation of human POMT1 and TMTC1. Transmembrane domain prediction was carried out using the algorithm PolyPhobius, as described in experimental procedures. The position of the predicted extracytoplasmic loops EL1, EL2, EL3 and EL5 is indicated by arrows and amino acid conservation in these loops is shown in the sequence logos.
Figure 4. Predicted membrane topology and amino acid conservation of human POMT1 and TMTC1. Transmembrane domain prediction was carried out using the algorithm PolyPhobius, as described in experimental procedures. The position of the predicted extracytoplasmic loops EL1, EL2, EL3 and EL5 is indicated by arrows and amino acid conservation in these loops is shown in the sequence logos.
Ijms 20 04842 g004
Figure 5. Membrane topology prediction of other eukaryotic glycosyltransferases of the GT-C superfamily. Arrows indicate the positions of the conserved areas that are depicted in Figure 6, corresponding to EL1, EL2, EL3, and EL5 (EL6 for ALG10, due to the two additional TMDs). The dilysine motif at the C-terminus of proteins with even numbers of predicted TMDs is shown.
Figure 5. Membrane topology prediction of other eukaryotic glycosyltransferases of the GT-C superfamily. Arrows indicate the positions of the conserved areas that are depicted in Figure 6, corresponding to EL1, EL2, EL3, and EL5 (EL6 for ALG10, due to the two additional TMDs). The dilysine motif at the C-terminus of proteins with even numbers of predicted TMDs is shown.
Ijms 20 04842 g005
Figure 6. Conserved domains in extracytoplasmic loops of GT-C enzymes. Conserved domains in EL1, EL2, EL3, and EL5 (EL6 for ALG10) of indicated GT-C members, presented as logos. The positions correspond to arrows in Figure 5.
Figure 6. Conserved domains in extracytoplasmic loops of GT-C enzymes. Conserved domains in EL1, EL2, EL3, and EL5 (EL6 for ALG10) of indicated GT-C members, presented as logos. The positions correspond to arrows in Figure 5.
Ijms 20 04842 g006
Figure 7. General membrane topology model of glycosyltransferase-C members. Location of the four conserved extracytoplasmic loops in the GT-C families are indicated by color-coded balls representing acidic amino acids. The 11 TMDs (numbered in black) present in all GT-C members are depicted. Inserted TMDs (in broken lines) increase the number of TMDs to 13 (blue) or to 12/14 with the C-terminus in the cytoplasm (red).
Figure 7. General membrane topology model of glycosyltransferase-C members. Location of the four conserved extracytoplasmic loops in the GT-C families are indicated by color-coded balls representing acidic amino acids. The 11 TMDs (numbered in black) present in all GT-C members are depicted. Inserted TMDs (in broken lines) increase the number of TMDs to 13 (blue) or to 12/14 with the C-terminus in the cytoplasm (red).
Ijms 20 04842 g007
Table 1. Assignment of human glycosyltransferases that uses dolichol linked substrate to CAZy glycosyltransferase families.
Table 1. Assignment of human glycosyltransferases that uses dolichol linked substrate to CAZy glycosyltransferase families.
EnzymesGT Family
ALG9, ALG12, PIGB, PIGZGT22
POMT1, POMT2GT39
PIGMGT50
ALG6, ALG8GT57
ALG3GT58
ALG10GT59
STT3A, STT3BGT66
PIGVGT76
DPY19L1, DPY19L2, DPY19L3, DPY19L4
TMTC1, TMTC2, TMTC3, TMTC4
GT98
GT105

Share and Cite

MDPI and ACS Style

Albuquerque-Wendt, A.; Hütte, H.J.; Buettner, F.F.R.; Routier, F.H.; Bakker, H. Membrane Topological Model of Glycosyltransferases of the GT-C Superfamily. Int. J. Mol. Sci. 2019, 20, 4842. https://doi.org/10.3390/ijms20194842

AMA Style

Albuquerque-Wendt A, Hütte HJ, Buettner FFR, Routier FH, Bakker H. Membrane Topological Model of Glycosyltransferases of the GT-C Superfamily. International Journal of Molecular Sciences. 2019; 20(19):4842. https://doi.org/10.3390/ijms20194842

Chicago/Turabian Style

Albuquerque-Wendt, Andreia, Hermann J. Hütte, Falk F. R. Buettner, Françoise H. Routier, and Hans Bakker. 2019. "Membrane Topological Model of Glycosyltransferases of the GT-C Superfamily" International Journal of Molecular Sciences 20, no. 19: 4842. https://doi.org/10.3390/ijms20194842

APA Style

Albuquerque-Wendt, A., Hütte, H. J., Buettner, F. F. R., Routier, F. H., & Bakker, H. (2019). Membrane Topological Model of Glycosyltransferases of the GT-C Superfamily. International Journal of Molecular Sciences, 20(19), 4842. https://doi.org/10.3390/ijms20194842

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop