Skip to contents

Check if the Glycans have the Given Motif(s)

have_motif.Rd

These functions check if the given glycans have the given motif(s).

  • have_motif() checks a single motif against multiple glycans

  • have_motifs() checks multiple motifs against multiple glycans

Technically speaking, they perform subgraph isomorphism tests to determine if the motif(s) are subgraphs of the glycans. Monosaccharides, linkages, and substituents are all considered.

Usage

have_motif(glycans, motif, alignment = NULL, ignore_linkages = FALSE)

have_motifs(glycans, motifs, alignments = NULL, ignore_linkages = FALSE)

Arguments

glycans

A 'glyrepr_structure' object, or an IUPAC-condensed structure string vector.

motif

A 'glyrepr_structure' object, an IUPAC-condensed structure string, or a known motif name (use available_motifs() to see all available motifs).

alignment

A character string. Possible values are "substructure", "core", "terminal" and "whole". If not provided, the value will be decided based on the motif argument. If motif is a motif name, the alignment in the database will be used. Otherwise, "substructure" will be used.

ignore_linkages

A logical value. If TRUE, linkages will be ignored in the comparison.

motifs

A character vector of motif names, IUPAC-condensed structure strings, or a 'glyrepr_structure' object.

alignments

A character vector specifying alignment types for each motif. Can be a single value (applied to all motifs) or a vector of the same length as motifs.

Value

  • have_motif(): A logical vector indicating if each glycan has the motif.

  • have_motifs(): A logical matrix where rows correspond to glycans and columns correspond to motifs. Row names contain glycan identifiers and column names contain motif identifiers.

Monosaccharide type

They can have different monosaccharide types ("concrete" or "generic", see glyrepr::get_mono_type()). The matching rules are:

  • When the motif is "generic", glycans are converted to "generic" type for comparison, allowing both concrete and generic glycans to match generic motifs.

  • When the motif is "concrete", glycans are used as-is, so only concrete glycans with matching monosaccharide names will match, while generic glycans will not match.

Examples:

  • Man (concrete glycan) vs Hex (generic motif) → TRUE (Man converted to Hex for comparison)

  • Hex (generic glycan) vs Man (concrete motif) → FALSE (names don't match)

  • Man (concrete glycan) vs Man (concrete motif) → TRUE (exact match)

  • Hex (generic glycan) vs Hex (generic motif) → TRUE (exact match)

Linkages

Obscure linkages (e.g. "??-?") are allowed in the motif graph (see glyrepr::possible_linkages()). "?" in a motif graph means "anything could be OK", so it will match any linkage in the glycan graph. However, "?" in a glycan graph will only match "?" in the motif graph. You can set ignore_linkages = TRUE to ignore linkages in the comparison.

Some examples:

  • "b1-?" in motif will match "b1-4" in glycan.

  • "b1-?" in motif will match "b1-?" in glycan.

  • "b1-4" in motif will NOT match "b1-?" in glycan.

  • "a1-?" in motif will NOT match "b1-4" in glycan.

  • "a1-?" in motif will NOT match "a?-4" in glycan.

Both motifs and glycans can have a "half-linkage" at the reducing end, e.g. "GlcNAc(b1-". The half linkage in the motif will be matched to any linkage in the glycan, or the half linkage of the glycan. e.g. Glycan "GlcNAc(b1-4)Gal(a1-" will have both "GlcNAc(b1-" and "Gal(a1-" motifs.

Alignment

According to the GlycoMotif database, a motif can be classified into four alignment types:

  • "substructure": The motif can be anywhere in the glycan. This is the default. See substructure for details.

  • "core": The motif must align with at least one connected substructure (subtree) at the reducing end of the glycan. See glycan core for details.

  • "terminal": The motif must align with at least one connected substructure (subtree) at the nonreducing end of the glycan. See nonreducing end for details.

  • "whole": The motif must align with the entire glycan. See whole-glycan for details.

When using known motifs in the GlycoMotif GlyGen Collection, the best practice is to not provide the alignment argument, and let the function decide the alignment based on the motif name. However, it is still possible to override the default alignments. In this case, the user-provided alignments will be used, but a warning will be issued.

Substituents

Substituents (e.g. "Ac", "SO3") are matched in strict mode. Both single and multiple substituents are supported:

  • Single substituents: "Neu5Ac-9Ac" will only match "Neu5Ac-9Ac" but not "Neu5Ac"

  • Multiple substituents: "Glc3Me6S" (has both 3Me and 6S) will only match motifs that contain both substituents, e.g., "Glc3Me6S", "Glc?Me6S", "Glc3Me?S"

  • "Glc3Me6S" will NOT match "Glc3Me" (missing 6S) or "Glc" (missing both)

For multiple substituents, they are internally stored as comma-separated values (e.g. "3Me,6S") and matched individually. Each substituent in the motif must have a corresponding match in the glycan, and vice versa.

Obscure linkages in motif substituents will match any linkage in glycan substituents:

  • Motif "Neu5Ac-?Ac" will match "Neu5Ac-9Ac" in the glycan

  • Motif "Glc?Me6S" will match "Glc3Me6S" in the glycan (? matches 3)

  • Motif "Glc3Me?S" will match "Glc3Me6S" in the glycan (? matches 6)

Implementation

Under the hood, the function uses igraph::graph.get.subisomorphisms.vf2() to get all possible subgraph isomorphisms between glycan and motif. color vertex attributes are added to the graphs to distinguish monosaccharides. For all possible matches, the function checks the following:

  • Alignment: using alignment_check()

  • Substituents: using substituent_check()

  • Linkages: using linkage_check()

  • Anomer: using anomer_check() The function returns TRUE if any of the matches pass all checks.

See also

count_motif(), count_motifs()

Examples

library(glyparse)
library(glyrepr)

(glycan <- o_glycan_core_2(mono_type = "concrete"))
#> <glycan_structure[1]>
#> [1] Gal(b1-3)[GlcNAc(b1-6)]GalNAc(a1-
#> # Unique structures: 1

# The glycan has the motif "Gal(b1-3)GalNAc"
have_motif(glycan, "Gal(b1-3)GalNAc")
#> [1] TRUE

# But not "Gal(b1-4)GalNAc" (wrong linkage)
have_motif(glycan, "Gal(b1-4)GalNAc")
#> [1] FALSE

# Set `ignore_linkages` to `TRUE` to ignore linkages
have_motif(glycan, "Gal(b1-4)GalNAc", ignore_linkages = TRUE)
#> [1] TRUE

# Different monosaccharide types are allowed
have_motif(glycan, "Hex(b1-3)HexNAc")
#> [1] TRUE

# Obscure linkages in the `motif` graph are allowed
have_motif(glycan, "Gal(b1-?)GalNAc")
#> [1] TRUE

# However, obscure linkages in `glycan` will only match "?" in the `motif` graph
glycan_2 <- parse_iupac_condensed("Gal(b1-?)[GlcNAc(b1-6)]GalNAc")
have_motif(glycan_2, "Gal(b1-3)GalNAc")
#> [1] FALSE
have_motif(glycan_2, "Gal(b1-?)GalNAc")
#> [1] TRUE

# The anomer of the motif will be matched to linkages in the glycan
have_motif(glycan_2, "GlcNAc(b1-")
#> [1] TRUE

# Alignment types
# The default type is "substructure", which means the motif can be anywhere in the glycan.
# Other options include "core", "terminal" and "whole".
glycan_3 <- parse_iupac_condensed("Gal(a1-3)Gal(a1-4)Gal(a1-6)Gal")
motifs <- c(
  "Gal(a1-3)Gal(a1-4)Gal(a1-6)Gal",
  "Gal(a1-3)Gal(a1-4)Gal",
           "Gal(a1-4)Gal(a1-6)Gal",
           "Gal(a1-4)Gal"
)

purrr::map_lgl(motifs, ~ have_motif(glycan_3, .x, alignment = "whole"))
#> [1]  TRUE FALSE FALSE FALSE
purrr::map_lgl(motifs, ~ have_motif(glycan_3, .x, alignment = "core"))
#> [1]  TRUE FALSE  TRUE FALSE
purrr::map_lgl(motifs, ~ have_motif(glycan_3, .x, alignment = "terminal"))
#> [1]  TRUE  TRUE FALSE FALSE
purrr::map_lgl(motifs, ~ have_motif(glycan_3, .x, alignment = "substructure"))
#> [1] TRUE TRUE TRUE TRUE

# Substituents
glycan_4 <- "Neu5Ac9Ac(a2-3)Gal(b1-4)GlcNAc"
glycan_5 <- "Neu5Ac(a2-3)Gal(b1-4)GlcNAc"

have_motif(glycan_4, glycan_5)
#> [1] FALSE
have_motif(glycan_5, glycan_4)
#> [1] FALSE
have_motif(glycan_4, glycan_4)
#> [1] TRUE
have_motif(glycan_5, glycan_5)
#> [1] TRUE

# Multiple substituents
glycan_6 <- "Glc3Me6S"  # has both 3Me and 6S substituents
have_motif(glycan_6, "Glc3Me6S")  # TRUE: exact match
#> [1] TRUE
have_motif(glycan_6, "Glc?Me6S")  # TRUE: obscure linkage ?Me matches 3Me
#> [1] TRUE
have_motif(glycan_6, "Glc3Me?S")  # TRUE: obscure linkage ?S matches 6S
#> [1] TRUE
have_motif(glycan_6, "Glc3Me")    # FALSE: missing 6S substituent
#> [1] FALSE
have_motif(glycan_6, "Glc")       # FALSE: missing all substituents
#> [1] FALSE

# Vectorization with single motif
glycans <- c(glycan, glycan_2, glycan_3)
motif <- "Gal(b1-3)GalNAc"
have_motif(glycans, motif)
#> [1]  TRUE FALSE FALSE

# Multiple motifs with have_motifs()
glycan1 <- o_glycan_core_2(mono_type = "concrete")
glycan2 <- parse_iupac_condensed("Gal(b1-?)[GlcNAc(b1-6)]GalNAc")
glycans <- c(glycan1, glycan2)

motifs <- c("Gal(b1-3)GalNAc", "Gal(b1-4)GalNAc", "GlcNAc(b1-6)GalNAc")
result <- have_motifs(glycans, motifs)
print(result)
#>                                   Gal(b1-3)GalNAc Gal(b1-4)GalNAc
#> Gal(b1-3)[GlcNAc(b1-6)]GalNAc(a1-            TRUE           FALSE
#> GlcNAc(b1-6)[Gal(b1-?)]GalNAc(?1-           FALSE           FALSE
#>                                   GlcNAc(b1-6)GalNAc
#> Gal(b1-3)[GlcNAc(b1-6)]GalNAc(a1-               TRUE
#> GlcNAc(b1-6)[Gal(b1-?)]GalNAc(?1-               TRUE