# LeishMANIAb: a comparative resource for Leishmania proteins

Gábor E Tusnády, András Zeke, Zsófia E Kálmán, Marie Fatoux, Sylvie

Ricard-Blum, Toby J Gibson, Laszlo Dobson

► To cite this version:

Gábor E Tusnády, András Zeke, Zsófia E Kálmán, Marie Fatoux, Sylvie Ricard-Blum, et al.. Leish-MANIAb: a comparative resource for Leishmania proteins. Database - The journal of Biological Databases and Curation, 2023, 2023, pp.2023:baad074. 10.1093/database/baad074. hal-04998307

HAL Id: hal-04998307

https://hal.science/hal-04998307v1

Submitted on 20 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# LeishMANIAdb: a comparative resource for *Leishmania* proteins

Gábor E. Tusnády<fnref n="1" />,<fnref n="2" />, András Zeke<fnref n="1" />, Zsófia E. Kálmán<fnref n="3" />, Marie Fatoux<fnref n="4" />,<fnref n="5" />, Sylvie Ricard-Blum<fnref n="1" />,<fnref n="4" />,
Toby J. Gibson<fnref n="6" /> and Laszlo Dobson<fnref n="1" />,<fnref n="6" />,*

Citation details: Tusnády, G.E., Zeke, A., Kálmán, Z.E. et al. LeishMANIAdb: a comparative resource for *Leishmania* proteins. *Database* (2024) Vol. 2024: article ID baad074; DOI: https://doi.org/10.1093/database/baad074

## Abstract

Leishmaniasis is a detrimental disease causing serious changes in quality of life and some forms can lead to death. The disease is spread by the parasite *Leishmania* transmitted by sandfly vectors and their primary hosts are vertebrates including humans. The pathogen penetrates host cells and secretes proteins (the secretome) to repurpose cells for pathogen growth and to alter cell signaling via host-pathogen protein-protein interactions). Here, we present LeishMANIAdb, a database specifically designed to investigate how *Leishmania* virulence factors may interfere with host proteins. Since the secretomes of different *Leishmania* species are only partially characterized, we collated various experimental evidence and used computational predictions to identify *Leishmania* secreted proteins to generate a user-friendly unified web resource allowing users to access all information available on experimental and predicted secretomes. In addition, we manually annotated host-pathogen interactions of 211 proteins and the localization/function of 3764 transmembrane (TM) proteins of different *Leishmania* species. We also enriched all proteins with automatic structural and functional predictions that can provide new insights in the molecular mechanisms of infection. Our database may provide novel insights into *Leishmania* host-pathogen interactions and help to identify new therapeutic targets for this neglected disease.

Database URL: https://leishmaniadb.ttk.hu/

## Introduction

Leishmaniasis is a neglected tropical disease causing severe symptoms, affecting around 1 million new people yearly, with annual deaths estimated to be around 60 000 (1). Although over 90% of cases occur in poor regions south of the Equator, due to climatic changes it also appears in new areas, and it has already shown up in Mediterranean European countries (2) and Texas, USA (3). To this date, no approved human vaccine is available and treatment is most effective at an early stage of the infection. *Leishmania* parasites are unicellular, flagellated trypanosomatids, belonging to the class Kinetoplastea. Upon infection, the amastigote stage pathogen (with reduced flagella) is engulfed by phagocytes, where it ends up in a stable parasitophorous vacuole that protects it (4). *Leishmania* cells then proliferate unhindered within host cells until egress and spreading to nearby phagocytes (5). The parasite secretes proteins that enter various parts of the cell (6). The secreted virulence factors can then interfere with cell signaling

by interacting with the host proteins: they increase glycolytic metabolism (7), perturb microbicidal pathways (8), escape the innate immune response, and repurpose macrophages for parasite replication (9) by disturbing cellular protein-protein interactions (PPIs). Interestingly, these mechanisms are somewhat unique to *Leishmania* among trypanosomes, which are mostly extracellular pathogens and do not enter host cells. In contrast, *Leishmania* secretes proteins which are critical for host cell subjugation, but how they enter the cytoplasm of host cells is still poorly understood.

In many distant, unrelated intracellular pathogens, ranging from viruses and bacteria to unicellular eukaryotes, the host targeted interactions are often mediated via Short Linear Motifs (SLiMs) (10). SLiMs are flexible protein segments composed of a restricted number of residues (between 3 and 10) that usually bind to structured protein domains. Their short length and structural flexibility enable them to bind to a wide range of domains. Cellular SLiMs typically bind their

<footnote n="1">¹Protein Bioinformatics Research Group, Institute of Enzymology, Research Centre for Natural Sciences, Magyar Tudósok körútja 2,
Budapest 1117, Hungary</footnote>

<footnote n="2">²Department of Bioinformatics, Semmelweis University, Tűzoltó u. 7, Budapest 1094, Hungary</footnote>

<footnote n="3">³Faculty of Information Technology and Bionics, Pázmány Péter Catholic University, Práter u. 50/A, Budapest 1083, Hungary</footnote>

<footnote n="4">⁴ICBMS UMR CNRS 5246, University Lyon 1, Rue Victor Grignard, Villeurbanne 69622, France</footnote>

<footnote n="5">⁵UMR CNRS 5086, University Lyon 1, 7 Passage du Vercors, Lyon 69367, France</footnote>

<footnote n="6">⁶Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstraße 1, Heidelberg 69117, Germany</footnote>

<footnote>*Corresponding author: Tel: +36 1 382 6727; E-mail: dobson.laszlo@ttk.hu</footnote>
targets with low micromolar affinity. These weak and tran-
sient interactions enable SLiMs to work in cooperative reg-
ulatory systems (11). Pathogens mimic host SLiMs to inter-
act with host cell proteins (10). Pathogen SLiMs often bind
with higher affinities than the cellular ones, outcompeting
the native interactions, permanently re-wiring the host reg-
ulation network. A few modulatory SLiMs have already
been discovered in eukaryotic pathogens, such as the Tox-
oplasma gondii MapK docking motif (12) and the stage-
specific (promastigote-amastigote) phosphorylation motifs
from Leishmania (13). In addition, several putative SLiMs
were recently detected in Leishmania, such as heparin-binding
sequences or RGD integrin-binding motifs though their func-
tion have not yet been confirmed yet (14).

Numerous studies investigated *Leishmania* secretomes. Most of them expose promastigotes to a heat shock and pH change (attempting to emulate the conditions that promote promastigote-to-amastigote stage transition) and then analyze the *Leishmania* conditioned medium by proteomics to identify secreted proteins (15), and measure their protein abundance or by transcriptomics to detect mRNA levels (16). While high-throughput experiments inherently suffer from a certain level of noise, experiments on individual proteins may be more reliable—in the case of *Leishmania* the vast majority focuses on leishmanolysin (GP63), a surface-anchored protease important for pathogenesis (17, 18). Furthermore, data were collected on different *Leishmania* species/strains identified via names and identifiers varying from one source to another, making a unified overview challenging. Another key step toward understanding the infection mechanism would be the identification of *Leishmania* surface proteins that can mediate the attachment of the pathogen to the host cell. Some surfaceome experiments were carried out on *Leishmania*-related species, and human host proteins binding to the surface of 24 strains of intact *Leishmania* have been identified (19). Besides the characterization of *Leishmania* secretomes, the identification of host–*Leishmania* PPIs is needed to narrow down virulence factors perturbing the host cell regulation to modules interfering with host proteins. SLiMs have low information content and simply scanning for matches to them in *Leishmania* secretomes may yield many false positives. Their structural and functional context, such as accessibility, conservation and localization, are all key elements to successfully identify those that may have a role in rewiring the host cell regulation. Notably, SLiMs also play a key role in maintaining housekeeping processes in *Leishmania*. Therefore, to find candidate SLiMs that may alter the host regulation, we need to discriminate SLiMs of proteins that reach the host cytoplasm or nucleus but limited information about these proteins are available. Currently, the only publicly available database dealing with *Leishmania* proteins is TriTrypDB (20), which is part of the VEuPathDB (21). TriTrypDB is a functional genomic resource for Trypanosomatidae, offering proteomic datasets; however, it does not focus on protein structure, protein motif search and interactions.

We developed LeishMANIAbd to expedite *Leishmania*
research by unifying scattered information from the litera-
ture in a user-friendly way and to extend available resources
by adding protein level information. We collected high-
throughput experiments and interaction studies on individ-
ual proteins and used various prediction methods to enrich
proteins with structural information.

Results

Selection of *Leishmania* proteomes and homology
mapping of various kinetoplastid proteins

We selected five *Leishmania* species (reference proteomes: *L. braziliensis*, *L. donovani*, *L. infantum*, *L. major* and *L. mexicana*), thirteen *Leishmania* strains (*Lbraziliensis*MHOMBR75M2903, *Lbraziliensis*MHOMBR75M2904, *Lbraziliensis*MHOMBR75M2904_2019, *LdonovaniBPK-282A1*, *LdonovaniCL-SL*, *LdonovaniHU3*, *LdonovaniLV9*, *LinfantumJPCM5*, *LmajorFriedlin*, *LmajorFriedlin2021*, *LmajorLV39c5*, *LmajorSD75.1* and *Lmexicana*MHOMGT2001U1103) and six related species (reference proteomes: *Bodo saltans*, *Leptomonas seymouri*, *Trypansoma brucei*, *Trypansoma cruzi*, *Trypanosoma rangeli* and *Trypanosoma theileri*) as an outgroup (22). *Leishmania* proteins were also cross-referenced to *TriTrypDB* (20). Around 30% of the cross-referenced proteins have different sequences deposited into these resources, and in most cases, the difference is due to the predicted position of the initiator methionine. For data compatibility, we always use the UniProt sequence version, but the conflicts are highlighted in *LeishMANIAbd*. When selecting the species we looked for those that have at least two strains deposited into *TriTrypDB*, and where cross-references and strain information were present and could be assigned with the least errors (Supplementary Table 1). We also performed a similarity search between these proteins and linked close homologs (only close kinetoplastid hits, see Methods) so annotations and predictions can be easily compared between them. All manual annotations and experimental data from different sources were mapped to these proteins. The 13 *Leishmania* strain proteomes were downloaded from *TriTrypDB*. Altogether *LeishMANIAbd* contains 40 537 searchable *Leishmania* proteins from reference proteomes, 108 766 proteins from different strains and 68 924 other kinetoplastid proteins to strengthen predictions.

Manual annotation of host-pathogen PPIs and TM
protein localization

We manually curated hundreds of proteins, using two
strategies.

The first type of annotation was the collection of host-pathogen PPI experiments on individual proteins, with the majority of them involving leishmanolysin (GP63). We collected 29 papers reporting 82 *Leishmania* PPIs with different hosts. Although experiments were mapped back to specific proteins, the results are also displayed on close homologs (with a note that the experimental data is derived from a different protein) resulting in 211 proteins that contain PPI data. Interactions were reported using the Minimum Information required for reporting a Molecular Interaction experiment MIMIx (23) community standard description.

The second type of manually curated data was the localiza-
tion and functional annotation of TM proteins. The aim was
to find surface proteins that may facilitate the infection, but
we annotated hundreds of other TM proteins with their local-
izations too. For this task, we used close homologous protein
groups. Altogether 342 protein families were annotated and
these annotations were shared between 3764 proteins (which
is 45.11% percent of the predicted TM proteomes and 9.28%
of all proteins of the 5 species combined).
## The definition of *Leishmania* secretome and protein localization is still incomplete

*Leishmania* not only exploits host-secretory pathways to distribute effectors, but also utilizes an unusual mechanism to deliver proteins to the cytosol of infected cells by releasing exosomes into the parasitophorous vesicle, which might fuse with the vesicular membrane to release their protein content (24). Therefore, computational methods based on signal peptides and localization predictions are not sufficient to predict *Leishmania* secretomes. To overcome this limitation, we also used high-throughput experiments (15, 25–28) to increase the coverage of *Leishmania* secretomes. Strikingly, the number of proteins in these secretomes varies to a large extent. Other datasets include proteins found in glycosomes (29), stage-dependent (promastigote/amastigote) phosphoproteomics (13), housekeeping gene localizations (30), exosome content (24), protein and mRNA abundance data (16, 31). When we mapped back all secretome and abundance experiments to *Leishmania infantum* (from close homologous proteins of other *Leishmania* species), the number of identified proteins ranges from 10 to 2000 (Figure 1A), and even when experimental conditions were similar they yielded highly different amounts of proteins. For example, pioneer secretome studies only provided a few hundred hits, while the latest ones are more inclusive with thousands of hits.

Gene duplication is often acting on protein families responsible for host-pathogen PPIs; therefore, we also collected proteins that are highly expanded. Notably, as all kinetoplastids have a polycistronic transcription system, the main way to amplify expression of critical proteins is through gene duplication. Thereby, highly expanded gene families can be directly mapped to functions critical for these parasites (32). In this case, we could discriminate between proteins with many paralogs found in other kinetoplastids versus *Leishmania*-exclusive amplified proteins. When we searched for close homologs of *Leishmania infantum* proteins, we found distinct amino acid transporter and cofactor families already expanded in all kinetoplastids including *Leishmania*. In contrast, amastins, leishmanolysin, 3'A2-related proteins, kinase-containing putative receptor proteins (and several uncharacterized proteins) seemed to be highly abundant in *Leishmania* proteomes compared to all kinetoplastids (Figure 1B).

Comparing complete proteomic datasets yielded only a small overlap. We defined (1) *Leishmania* novelty proteins, which are proteins without close homologs in SwissProt, without characterized Pfam domains, and expanded in *Leishmania infantum* (compared with other kinetoplastids); (2) abundant proteins, which are proteins showing increased abundance upon infection; and (3) secreted proteins experimentally identified in at least two secretome experiments. These definitions provided markedly different protein sets, with some overlap between secreted and abundant proteins (611 proteins) and with only 22 proteins contained in all datasets (Figure 1C).

## AlphaFold2 provides an alternative way to estimate structural features

We used different methods to predict the structural features of proteins. Classical sequence-based methods can detect globular domains (33), TM regions (34) and intrinsically disordered regions (IDRs) (35). However, the use of AlphaFold2 (AF2) (36) provides alternative ways to obtain structural

information. In LeishMANIAbd, we used structures available in the AlphaFold database (37) (however, we could not find 3192 proteins (~6% of all *Leishmania* proteins in LeishMANIAbd)). We not only displayed the predicted 3D structure of the proteins, but also information derived from the AF2 models, such as the secondary structures and the position of the lipid bilayer for membrane proteins using the method introduced in the TmAlphaFold database (38). Although AF2 was originally built to predict protein structure, the scientific community quickly realized it is as much (if not more) efficient at predicting protein disorder (38). To analyze IDRs, we displayed predicted local distance difference test (pLDDT) values and relative surface accessibility from AF2. For IDR prediction in TM proteins, we tailored MemDis (39) to incorporate features from AF2 instead of sequence-based predictors (see Methods).

## Short linear motif candidates that may hijack host cell regulation

We scanned *Leishmania* proteins for SLiMs using the regular expressions stored in the Eukaryotic Linear Motif (ELM) resource (40). Scanning SLiMs alone would mostly yield false positive hits, so we developed a scoring system that ranges from 0 to 1, and that takes into account most information we collected. We aimed to develop a scoring system where conservation and accessibility/disorder has a reasonably high weight, while keeping in mind that proteomic experiments and localization information are a good way to narrow down the potentially large number of false positive hits. Unfortunately, due to the lack of data, in the case of *Leishmania*, it is not possible to construct a benchmark set to evaluate motif scores. We can still assume that a good starting point can be when most predictions and proteomic data agree. Considering *Leishmania infantum* alone, we detected over a million putative motifs, from which 1.21% had a score above 0.85, on 343 proteins (Figure 1D).

## The LeishMANIAbd web resource

To visualize all the collected and calculated information, we developed an open-access resource. In LeishMANIAbd, users can search for proteins using their UniProt Accession (AC), Entry name (formerly ID), gene name and protein name. We also provide several protein sets as examples to help users browsing the database. Currently, proteins are sorted based on (1) species: *L. braziliensis*, *L. donovani*, *L. infantum*, *L. major* and *L. mexicana*; (2) manual curation data; (3) experimental data: secreted proteins, protein abundance/mRNA level data, proteins with any kind of experimental data listed above; (4) computationally predicted information: proteins expanded in *Leishmania* (score ≥ 0.8—see Methods, Supplementary Material), transmembrane (TM) proteins, proteins with high disordered content (at least 70% predicted disorder), proteins with high-scoring linear motifs (score ≥ 0.85), and novel kinetoplastid proteins (proteins without SwissProt homologs or Pfam domains). After searching (or selecting a protein set), users can further narrow their selection by choosing any other criterion (Figure 2A).

The entry page for each protein consists of up to 10 sections, which are only visible if they contain data. The 'Quick info' displays the protein name, species, cross-references, and its number of amino acid residues. Data curation appears under the 'Annotations' section. PPI (curated
<figure><img src="image_2.png" /><figcaption>**Figure 1.** LeishMANIAdb content. All data were calculated on *Leishmania infantum*. A: Number of proteins in different proteomic datasets (purple: promastigote secretome, red: amastigote secretome, orange: exosome, green: housekeeping genes, blue: higher protein abundance level upon infection). B: Number of kinetoplastid and Leishmania close homologs. Each dot represents a protein (red: at least 80% of close homologs are in Leishmania, blue: other proteins). Green circles represent distinctive groups. C: Overlap between abundant, secreted and 'Leishmania novelty' proteins (for more detail see text). D: Distribution of all predicted SLiMs with different scores. Red marks candidate motifs above 0.85 cutoff (for more details see text).</figcaption></figure>

at the MIMIx level), localization, and function annotations are mirrored from close homologous proteins. We also display functional annotations for *Leishmania donovani* (and close homologs) by Jardim *et al.* (30) (last update in 2018). The ‘Localization’ section contains high-throughput experiment data—promastigote and amastigote secretion, an exosome experiment and the glycosome. Protein localization, signal peptide and glycosylphosphatidyl (GPI) anchor predictions are also displayed here. Since the reliability of both predictions and experiments may vary, we also display all this data for close homologous proteins, so users can quickly check the robustness of information—by checking if different sources on closely related proteins agree. Furthermore, we also collected Gene Ontology (GO) (41) annotations for cellular compartments. In this case, the specificity of the term (how deep it is on the tree) is shown in the level column. GO annotations are collected for all close homologous proteins too and the number of occurrences of each term is displayed. We highlighted terms that are associated with the

inspected protein itself that is displayed on the page. The ‘Abundance’ module can display the mRNA and protein level experiments: static/single point (upregulated or not) or time-course experiments (e.g. mRNA and protein levels available for 7 timesteps across 120 h (16), Figure 2B). In the “Expansion” section, the number of close homologs are displayed by species with a color-code to identify *Leishmania* intracellular and extracellular/free-living relatives (Figure 2C). The ‘Sequence features’ displays various information (Figure 2D). At the top, the gapless multiple sequence alignment (MSA) of proteins from the reference proteomes is visible. In this alignment gaps from the entry protein were removed (the original alignment with all strains can be downloaded) so other protein features could be visualized. Protein disorder, secondary structures, TM topology prediction, domains, signal peptides and GPI anchors, and stage-dependent phosphorylation are also displayed. Predicted SLiMs are shown with a color-coded score (see Methods, Supplementary Material). In the ‘Structure’ section, the AF2 predicted structure is available (with
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Layout of LeishMANIAbd. A: The search/browse result menu. B: Expression section of the entry page. C: Expansion section. D: Sequence features section. E: Putative motif mimicry section. F: Function section.</figcaption></figure>

the position of the membrane domain for TM proteins). The ‘Putative motif mimicry’ section is the table format version of SLiMs from the “Sequence features” module (Figure 2E). The ‘Function’ section contains GO Molecular Function and Biological Process terms. As done for the Cellular Component, terms are transferred from close homologous proteins, and they can be sorted based on their specificity (how deep they are located in the tree) and occurrences considering close homologs (Figure 2F). Finally, each result from a BLAST search against SwissProt and kinetoplastids relatives are listed in the ‘Homologs’ table (not only close homologs, see Methods).

For each protein the full MSA, high-throughput experiments, annotation and predicted sequential features can be downloaded from the bottom of the page. Batch download is also available to download the full database or different protein sets.

# Discussion

## Reliability of data

In LeishMANIAbd, we aimed to collect high-throughput experimental data, PPI data on individual proteins, predictions and localization information based on close homologs. We noticed that the amount of data from MS experiments differs highly, and therefore likely the quality also varies. Not unexpected from high-throughput techniques applied to less-studied organisms, this can be attributed to the quality of sample preparation, MS experimental techniques and most

likely to the sequences in the background databases. One striking finding was that the secretory datasets contain a large number of proteins that are likely to take part in the housekeeping processes of *Leishmania* cells, such as cytoskeletal proteins, nuclear histones and metabolic enzymes. Exosomes are known to contain a relatively high amount of “background” proteins leaking from the cytosol of cells. Another explanation is that several housekeeping genes (such as intracellular chaperones and enzymes) are moonlighting proteins, they are generally constitutively expressed and have high levels of expression, while they are fulfilling other functions outside the cells (42). Due to the lack of comparative studies, we cannot assess the enrichment ratios of secreted molecules, to see if there is selective exosomal packaging of a well-defined subset of *Leishmania* proteins. However, *Leishmania* exosomal-like secretion also differs from the typical exosomal sorting seen in other eukaryotic organisms because budding primarily initiates at the cell membrane, and not inside multivesicular bodies (endosomes). Therefore, it is equally possible that in Leishmaniids, the budding is non-selective for its cytoplasmic cargos. Instead, it would be initiated by cell surface receptors and primarily serve as a defense mechanism against membrane-attached host complement and other immune complexes, removing them before they could damage the parasite membrane. Currently, testing of the latter hypothesis is impossible, since only soluble components, but not the integral membrane proteins of *Leishmania* exosomes have been studied in depth in the above cited studies.
From a computational point of view, predicting any features on *Leishmania* proteins might be highly challenging, as methods established were mostly trained on sequences that show little or no similarity to *Leishmania* proteins. The five *Leishmania* reference proteomes contain 10 267 uncharacterized proteins combined, which is ~25% of LeishMANIAbdb. TmAlphaFold provides an objective quality measurement option for α-helical membrane proteins. When we compared the TM proteome of *Homo sapiens* and *Leishmania infantum*, we noticed that the ratio of good and excellent quality structures was much lower in *Leishmania*, probably caused by the different coverage of kinetoplastid and human structures deposited into the PDB (Figure 3A).

## Case studies

LeishMANIAbdb can be utilized for different purposes and can be a good starting point for various analyses. We selected three examples that highlight some use cases of the resource.

Using the Browse menu, after selecting a category, users can further narrow down their search for proteins selecting additional categories to refine the results. For instance, if users are looking for *Leishmania* SLiMs that may alter or rewire host cell regulation network, they can look for proteins that were experimentally proven to be secreted, and then select proteins with disordered regions because SLiMs are mostly located in IDRs. ‘Kinetoplastid novelty’ selection ensures that the protein and its domains are not present in organisms belonging to other lineages, while *Leishmania* novelty/expansions select proteins that are new or highly expanded in *Leishmania* species. Last, by selecting high-scoring motifs, users get a list of proteins where the motif is most likely to be functional (Figure 3B shows the Venn diagram of the selection). These proteins may be an interesting starting point for further analyses.

When performing systematic searches to identify possible parasite hits of integrin ligand motifs (that only function in the host, as kinetoplastids have no integrins), we identified a striking set of examples in a family of poorly-known *Leishmania* genes called 3′A2 related ORFs. This kinetoplastid-specific family of genes is actually expanded in *Leishmania* species together with the canonically secreted A2 proteins, which are known pathogenicity factors (43). While the actual sequences of these proteins are poorly conserved and very little is known about their subcellular location, the *Leishmania* versions have at least one TM region and a C-terminal cytoplasmic tail, with an N-terminal signal peptide (or possibly another TM segment). Nevertheless, in the predicted, largely disordered extracellular segment we observed multiple, short, conserved stretches that may have amyloidogenic properties (high Val, Ala and Gly content, upon visual inspection), presumably capable of oligomerization and amphiphilic interaction with membranes (Figure 3C). A highly conserved cysteine residue preceding the first amyloidogenic sequence might help the homodimerization by forming a disulfide bridge with neighboring 3′A2 related protein. Strikingly, in *Leishmania infantum* and *Leishmania donovani* (both species capable of causing visceral leishmaniasis), the N-terminus of these proteins carries canonical RGD (Arg-Gly-Asp) sequences, immediately after the putative signal peptide cleavage site. In addition, *Leishmania donovani* and *Leishmania infantum* proteins contain an NGR motif where asparagine deamidation might yield an isoDGR motif. If these proteins are

expressed on the cell surface, they might bind to host integrins in an oligomeric state, and might even attack the host membrane as if it were a β-barrel pore-forming toxin. However, much more experiments are needed to test any of these hypotheses.

Amastins are a large family of kinetoplastid-specific membrane proteins that belong to the broader claudin-like superfamily, implicated in the maintenance of parasitophorous vacuoles (44). Accordingly, the majority of amastins have four tightly packed TM segments, with cytosolic tail regions. Similarly to their vertebrate counterparts that form tight cell-cell junctions by complex oligomerization processes, amastins might also engage in a variety of interactions with internal as well as external, host proteins (44). Although their exact function is not known, among the 221 identified amastins with 4 TM regions we looked for SLiMs that occur in multiple proteins. Predicted SLiMs (within disordered regions) were packed in their cytoplasmic tail regions (Figure 3D). Since these regions face inward the parasite, we further narrowed hits based on their binding domain to be present in *Leishmania*. We identified multiple potential phosphorylation sites and PPI motifs, such as SH3 ligands (*Leishmania* species do encode SH3 domain proteins) as well as vesicular trafficking signals. The tail region of amastins seems to be highly variable, likely acting as a hotspot in the pathogen-host arms race.

## Comparison with other resources

In the past decades, several databases were built to investigate *Leishmania*; however, they most unfortunately often offline and no longer updated. LeishCyc (45) focused on biochemical pathways. LeishDB (46) included coding genes and non-RNAs and provided new annotation to them. The cysteine protease database in *Leishmania* species (47) was designed to find data related to cysteine protease and LeishBase (48) was a structural database. There are a few active databases: Leish-ExP (http://www.hpppi.iicb.res.in/Leish-ex) (which has not so far been published in a peer-reviewed journal) contains proteins exclusively present in *Leishmania*. Leish-ExP incorporates localization tools, includes GO annotations and calculates physico-chemical properties of proteins. LmSmdB (49) focuses on metabolic and biosynthetic pathways. TriTrypDB (20) is a kinetoplastid database that is part of the VEuPathDB resource (21). These databases contain a lot of experimental data and various tools to analyze eukaryotic pathogens, but they are mostly focused on genomic data—although proteomic datasets, and some protein prediction algorithms are also incorporated.

There are also a handful of databases that include information on host-pathogen interactions: HPIDB (50), PHIDIAS (51) and PHI-base (52) contain information about PPIs between the host and pathogen, while ImitateDB (53) specifically focuses on motif mimicry. These resources contain no or very little data about *Leishmania*.

In LeishMANIAbdb, our main goal was to include protein information relevant to the infection and to complement previously established and still available resources. We included several proteomic datasets and enriched experimental information with state-of-the-art prediction tools. Still, the most powerful way to explore uncharted proteomes is to inspect MSAs and check for conserved residues and
<figure><img src="image_4.png" /><figcaption>**Figure 3.** A: Distribution of membrane protein quality levels of AlphaFold structure in Homo sapiens and Leishmania infantum. B: Upset diagram of proteins that are (1) secreted, (2) novel kinetoplastid, (3) expanded (or new) in Leishmania, (4) disordered, (5) contain candidate SLiMs. C: left: Multiple Sequence Alignment of 3'A2 related proteins (alignment is available under UniProt AC: E9AGZ3). Amyloidogenic regions, conserved cysteine and integrin-binding motifs are highlighted; right: proposed topology of 3'A2 related proteins. D: Frequent SLiMs in the cytoplasmic tail regions of amastins (the numbers denote the unique/total occurrences).</figcaption></figure>

regions—LeishMANIAbdb contains precalculated alignments for all proteins. We also added hundreds of annotations to thousands of proteins, including localization and interaction information. While several databases seem to be shut down after a couple of years, our laboratory hosts several resources and we routinely update them. We plan to do so with LeishMANIAbdb as well as to expand its repertoire to host-Leishmania interactions involving glycans and glycolipids, which play major roles in the infection.

## Methods

### Resources

Protein sequences were retrieved from UniProtKB (release 2022_05) (22) and from TriTrypDB, release_59 (20) based on the UniProt cross-references (*L. braziliensis*, *L. donovani*, *L. infantum*, *L. major*, *L. mexicana*, *Bodo saltans*, *Lep-tomonas seymouri*, *Trypanosoma brucei*, *Trypanosoma cruzi*, *Trypanosoma rangeli* and *Trypanosoma theileri*). Homologs
in other kinetoplastids and in SwissProt were searched with BLAST using E-value: <10⁻⁵; sequence identity >20%; and coverage >50%. In the 'Homologs' section, all results are displayed, and this information was used to calculate the "out-group score" component of the motif score (see below and in Supplementary Material). In most other sections (and calculation), we only used homologous proteins until the first non-kinetoplastid SwissProt hit considering sequence identity (termed as 'close homologs', Supplementary Figure S1). Further similar kinetoplastid proteins were therefore considered as a different homology group, and this way huge superfamilies with a common ancestor from other species are split into smaller families. This consideration seemed to be useful for calculating expanded proteins. Furthermore, we downloaded strains belonging to the five selected *Leishmania* species from TriTrypDB. In this case, a more stringent condition was used in BLAST, by setting E-value: 10⁻⁵; sequence identity > 80%; coverage > 80%. All kinetoplastid species and strains were used to calculate motif conservation.

We prepared three different types of MSAs using ClustalΩ (54): (1) 'non-redundant' MSA using close homologous proteins from kinetoplastid reference proteomes; (2) the same MSA but with gaps removed from the 'reference' protein that is currently displayed on the webpage; and (3) a more redundant MSA using close homologous kinetoplastid proteins in all species and strains (used to calculate motif conservation).

High-throughput experiments were first mapped to the corresponding protein using the identifier provided in the original paper, then mirrored to close *Leishmania* close homologs if their sequence were identical.

IDRs were predicted using IUPred3 with default settings (long) (35) and using the AF2 models' pLDDT and accessibility values—the latter was calculated by DSSP 3.1.5 with default settings (32), normalized using maximum values calculated as in Tien et al. (55), the exposed value threshold defined as suggested by Rost et al. (56). In the case of TM proteins, IDRs were also predicted by MemDis 1.0, sensitive settings (39). In this in-house modified version, the Position-Specific Scoring Matrices) were generated using kinetoplastid sequence library instead of SwissProt sequences. Furthermore secondary structure and accessibility were calculated using AlphaFold2 instead of sequence based predictors trained on distant proteins that was used in the original version. TM protein topology was predicted by CCTOP (34), however to minimize sporadic erroneous predictions, after an initial prediction we performed a constrained iteration where the topologies of close homologous proteins were used as a constraint. Using this approach, closely related proteins will likely have the same topology. Secondary structure elements derived from AF2 structures are also displayed. Pfam domains were identified using InterPro (33). Protein localization was displayed from the following tools: GO (41), DeepLoc (57) and SignalP6.0 (58). NetGPI (59) was used to predict GPI-anchors (all prediction results are displayed, therefore in case of a contradiction it is up to the user to evaluate the results). For these predictions, we used default settings and thresholds.

To detect SLiMs that may alter or rewire host cell regulation, we used the regular expressions from ELM (40) on all *Leishmania* sequences. We defined different contextual filters and merged them into a single score to rank motifs (for more details, see Supplementary Material): (1) Disordered: The score is the average of the IUPred3, AF2-based pLDDT

and accessibility values. These disordered scores were first transformed so they range from 0 to 1, with 0.5 being the threshold, before calculating their mean; (2) Conservation of the motif was checked among close homologs with some permission for slight misalignment, and penalizing motifs that are present across all kinetoplastids—notably, in this case proteins from different *Leishmania* strains were also considered; (3) Localization: we used a simplified (intracellular/extracellular) distinction. Motif localization was determined using ELM GO annotations, secretion information and CCTOP prediction, while the binding domain localization was determined from TOPDOM (60). We looked for motif-domain pairs where they both have the same simplified (in/out) localization; (4) mRNA level: using transcriptomic experiments about expression data; (5) protein level: from experiments about protein abundance; (6) Secretion score based on secretome experiments; (7) Expansion score: reflecting how much the protein is expanded in *Leishmania* species (strains not included) compared to all kinetoplastids; (8) Outgroups score favoring proteins without homologs in SwissProt.

Structure data reflects structure data deposited in the PDB (61) before 26.03.2023, AlphaFold database (v3) and the TmAlphaFold database (v1). To generate Figure 3A we normalized the following graph: https://tmalphafold.ttk.hu/stati stics?parameters=MSZjMzllMDA0ME5EQXINREU1TVRJ ek1DWXINVFEzTkRnek5qVTJKakkwTVRVNU1Ua3hNR GdtTXpJM05qZz0%3D. All other data was downloaded in October, 2022 from the source databases.

## Manual curation

We manually curated hundreds of proteins, using two strategies. First, we searched PubMed and Google scholar for 'Leishmania' host-PPI and manually processed the results. Each protein in the experiments was mapped to the corresponding UniProt entry. Then, we mapped interaction data to the five *Leishmania* proteomes. When the experiment was performed on a protein from a different species, we mirrored it to the closest homology group in LeishMANIAbd, and we also indicated on the webpage that the experiment is from a different protein. All interactions were reported according to the community standard MIMIX level (23).

For annotating function/localization, we used the BLAST and alignment results from SwissProt (release 2022_05) (22) and kinetoplastid species, and information on annotated SwissProt entries if they were found. We also used publications when they were available (references are added to the website and to the downloaded files). We also used high-throughput studies on surface proteins (62, 63), and if a homologous protein was measured on the surface we took into account this information. We used prediction/annotation tools (Gene Ontology, DeepLoc, SignalP and NetGPI). These annotations were not made on proteins one by one, but rather for larger sets of kinetoplastid proteins belonging to the same family. We manually processed the entries using this approach, taking distant homologues, domain architectures and conservation patterns into consideration.

## Website design

The LeishMANIAbd website is written in PHP (v8.0) using the Laravel (v9.19) framework. All downloaded, predicted, or calculated data are stored in a local MySQL (v8.0) database.
To visualize sequence features over amino acid sequences, we developed a javascript package using React (18.2), while 3D structures are visualized using the original (for non-TM proteins) or a locally modified version of Mol\* (64) for TM proteins that can display the membrane bilayer (the modified version is available at: https://git.enzim.ttk.hu/web/TmMolStar). The modified version can display the membrane as two planes around the investigated TM protein using the results of TMDET 2.0 (65).

## Supplementary material

Supplementary material is available at Database online.

## Funding

This project has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie Grant Agreement No 101028908; Ministry of Innovation and Technology of Hungary from the National Research, Development and Innovation Fund [K132522].

## Conflict of interest

None declared.

## Acknowledgements

The authors would like to thank the help of the Gibson team at the European Molecular Biology Laboratory, especially Jesus Alvarado-Valverde for helpful suggestions and discussions and Philippe Esterre (Pasteur Network, Institut Pasteur, Paris) for his help in selecting a number of Leishmania interactions to be curated.

## References

1. Torres-Guerrero,E., Quintanilla-Cedillo,M.R., Ruiz-Esmenjaud,J. et al. (2017) Leishmaniasis: a review. *F1000Res.*, **6**, 1–15.

2. Gianechci,E. and Montomoli,E. (2020) The enemy at home: leishmaniasis in the Mediterranean basin, Italy on the focus. *Expert Rev. Anti. Infect. Ther.*, **18**, 563–577.

3. McIlwee,B.E., Weis,S.E. and Hosler,G.A. (2018) Incidence of endemic human cutaneous leishmaniasis in the United States. *JAMA Dermatol.*, **154**, 1032–1039.

4. Arango Duque,G., Jardim,A., Gagnon,É. et al. (2019) The host cell secretory pathway mediates the export of Leishmania virulence factors out of the parasitophorous vacuole. *PLoS Pathog.*, **15**, 1–25.

5. Real,F., Florentino,P.T.V., Reis,L.C. et al. (2014) Cell-to-cell transfer of *Leishmania amazonensis* amastigotes is mediated by immunomodulatory LAMP-rich parasitophorous extrusions. *Cell. Microbiol.*, **16**, 1549–1564.

6. Atayde,V.D., Aslan,H., Townsend,S. et al. (2015) Exosome secretion by the parasitic protozoan leishmania within the sand fly midgut. *Cell Rep.*, **13**, 957–967.

7. Ohms,M., Ferreira,C., Busch,H. et al. (2021) Enhanced glycolysis is required for antileishmanial functions of neutrophils upon infection with. *Front. Immunol.*, **12**, 1–14.

8. Matheoud,D., Moradin,N., Bellemare-Pelletier,A. et al. (2013) Leishmania evades host immunity by inhibiting antigen cross-presentation through direct cleavage of the SNARE VAMP8. *Cell Host Microbe*, **14**, 15–25.

9. Atayde,V.D., Hassani,K., da Silva Lira Filho,A. et al. (2016) Leishmania exosomes and other virulence factors: impact on innate immune response and macrophage functions. *Cell. Immunol.*, **309**, 7–18.

10. Davey,N.E., Travé,G. and Gibson,T.J. (2011) How viruses hijack cell regulation. How viruses hijack cell regulation. *Trends Biochem. Sci.*, **(2011)**, **36**: 159–169.

11. Van Roey,K., Uyar,B., Weatheritt,R.J. et al. (2014) Short linear motifs: ubiquitous and functionally diverse protein interaction modules directing cell regulation. *Chem. Rev.*, **114**, 6733–6778.

12. Pellegrini,E., Palencia,A., Braun,L. et al. (2017) Structural basis for the subversion of MAP kinase signaling by an intrinsically disordered parasite secreted agonist. *Structure*, **25**, 16–26.

13. Tsigankov,P., Gherardini,P.F., Helmer-Citterich,M. et al. (2013) Phosphoproteomic analysis of differentiating leishmania parasites reveals a unique stage-specific phosphorylation motif. *J. Proteome Res.*, **12**, 3405–3412.

14. Peysselon,F., Launay,G., Lisacek,F. et al. (2013) Comparative analysis of Leishmania exoproteomes: implication for host-pathogen interactions. *Biochim Biophys Acta – Proteins Proteom*, **1834**: 2653–2662.

15. Cuervo,P., De Jesus,J.B., Saboia-Vahia,L. et al. (2009) Proteomic characterization of the released/secreted proteins of *Leishmania* (*Viannia*): *braziliensis* promastigotes. *J. Proteomics*, **73**, 79–92.

16. Lahav,T., Sivam,D., Volpin,H. et al. (2011) Multiple levels of gene regulation mediate differentiation of the intracellular pathogen Leishmania. *FASEB J.*, **25**, 515–525.

17. Guay-Vincent,-M.-M., Matte,C., Berthiaume,A.-M. et al. (2022) Revisiting Leishmania GP63 host cell targets reveals a limited spectrum of substrates. *PLoS Pathog.*, **18**, 1–17.

18. Gregory,D.J., Godbout,M., Contreras,I. et al. (2008) A novel form of NF-kappaB is induced by Leishmania infection: involvement in macrophage gene expression. *Eur. J. Immunol.*, **38**, 1071–1081.

19. Fatoux-Ardore,M., Peysselon,F., Weiss,A. et al. (2014) Large-scale investigation of Leishmania interaction networks with host extracellular matrix by surface plasmon resonance imaging. *Infect. Immun.*, **82**, 594–606.

20. Shanmugasundram,A., Starns,D., Böhme,U. et al. (2023) TriT-rypDB: an integrated functional genomics resource for kinetoplastida. *PLoS Negl. Trop. Dis.*, **17**, 1–20.

21. Amos,B., Aurrecoechea,C., Barba,M. et al. (2022) VEuPathDB: the eukaryotic pathogen, vector and host bioinformatics resource center. *Nucleic Acids Res.*, **50**, D898–D911.

22. UniProt Consortium. (2023) UniProt: the universal protein knowledgebase in 2023. *Nucleic Acids Res.*, **51**, D523–D531.

23. Orchard,S., Salwinski,L., Kerrien,S. et al. (2007) The minimum information required for reporting a molecular interaction experiment (MIMIx). *Nat. Biotechnol.*, **25**, 894–898.

24. Silverman,J.M., Clos,J., De'oliveira,C.C. et al. (2010) An exosome-based secretion pathway is responsible for protein export from Leishmania and communication with macrophages. *J. Cell Sci.*, **123**, 842–52.

25. Hassani,K., Antoniak,E., Jardim,A. et al. (2011) Temperature-induced protein secretion by *Leishmania mexicana* modulates macrophage signalling and function. *PLoS One*, **6**, 1–10.

26. Forrest,D.M., Batista,M., Marchini,F.K. et al. (2020) Proteomic analysis of exosomes derived from procyclic and metacyclic-like cultured *Leishmania infantum chagasi*. *J. Proteomics*, **227**, 1–10.

27. Silverman,J.M., Chan,S.K., Robinson,D.P. et al. (2008) Proteomic analysis of the secretome of Leishmania donovani. *Genome Biol.*, **9**, 1–21.

28. Pissarra,J., Pagniez,J., Petitdidier,E. et al. (2022) Proteomic analysis of the promastigote secretome of seven leishmania species. *J. Proteome Res.*, **21**, 30–48.

29. Jamdhade, Pawar,H., Chavan,S. et al. (2015) Comprehensive proteomics analysis of glycosomes from *Leishmania donovani*. *OMICS*, **19**: 157–170.
30. Jardim,A., Hardie,D.B., Boitz,J. et al. (2018) Proteomic profiling of Leishmania donovani promastigote subcellular organelles. *J. Proteome Res.*, **17**, 1194–1215.

31. Pescher,P., Blisnick,T., Bastin,P. et al. (2011) Quantitative proteome profiling informs on phenotypic traits that adapt Leishmania donovani for axenic and intracellular proliferation. *Cell. Microbiol.*, **13**, 978–91.

32. Jackson,A.P., Otto,T.D., Aslett,M. et al. (2016) Kinetoplas- tid phylogenomics reveals the evolutionary innovations associated with the origins of parasitism. *Curr. Biol.*, **26**, 161–172.

33. Paysan-Lafosse,T., Blum,M., Chuguransky,S. et al. (2023) InterPro in 2022. *Nucleic Acids Res.*, **51**, D418–D427.

34. Dobson,L., Reményi,I. and Tusnády,G.E. (2015) CCTOP: a con- sensus constrained TOPology prediction web server. *Nucleic Acids Res.*, **43**, W408–12.

35. Erdős,G., Pajkos,M. and Dosztányi,Z. (2021) IUPred3: predic- tion of protein disorder enhanced with unambiguous experimental annotation and visualization of evolutionary conservation. *Nucleic Acids Res.*, **49**, W297–W303.

36. Jumper,J., Evans,R., Pritzel,A. et al. (2021) Highly accurate protein structure prediction with AlphaFold. *Nature*, **596**, 583–589.

37. Varadi,M., Anyango,S., Deshpande,M. et al. (2022) AlphaFold protein structure database: massively expanding the structural coverage of protein-sequence space with high-accuracy models. *Nucleic Acids Res.*, **50**, D439–D444.

38. Dobson,L., Szekeres,L.I., Gerdán,C. et al. (2023) TmAlphaFold database: membrane localization and evaluation of AlphaFold2 predicted alpha-helical transmembrane protein structures. *Nucleic Acids Res.*, **51**, D517–D522.

39. Dobson,L. and Tusnády,G.E. (2021) MemDis: predicting disor- dered regions in transmembrane proteins. *Int. J. Mol. Sci.*, **22**, 1–7.

40. Kumar,M., Michael,S., Alvarado-Valverde,J. et al. (2022) The eukaryotic linear motif resource: 2022 release. *Nucleic Acids Res.*, **50**, D497–D508.

41. The Gene Ontology Consortium. (2019) The gene ontology resource: 20 years and still GOing strong. *Nucleic Acids Res.*, **47**, D330–D338.

42. Jeffery,C.J. (2018) Protein moonlighting: what is it, and why is it important? *Philos. Trans. R. Soc. Lond. B Biol. Sci.*, **373**, 1–8.

43. Zhang,W.W. and Matlashewski,G. (2001) Characterization of the A2-A2rel gene cluster in *Leishmania donovani*: involvement of A2 in visceralization during infection. *Mol. Microbiol.*, **39**, 935–948.

44. de Paiva,R.M.C., Grazielle-Silva,V., Cardoso,M.S. et al. (2015) Amastin knockdown in *Leishmania braziliensis* affects parasite- macrophage interaction and results in impaired viability of intra- cellular amastigotes. *PLoS Pathog.*, **11**, 1–24.

45. Saunders,E.C., MacRae,J.I., Naderer,T. et al. (2012) LeishCyc: a guide to building a metabolic pathway database and visualization of metabolomic data. *Methods Mol. Biol.*, **881**, 505–529.

46. Torres,F., Arias-Carrasco,R., Caris-Maldonado,J.C. et al. (2017) LeishDB: a database of coding gene annotation and non- coding RNAs in *Leishmania braziliensis*. *Database*, **2017**, 1–11.

47. Rana,S., Dikhit,M.R., Rani,M. et al. (2012) CPDB: cysteine pro- tease annotation database in *Leishmania* species. *Integr. Biol.*, **4**, 1351–1357.

48. Kaur,S., Patel,H., Sharma,V. et al. (2009) LeishBase: *Leishmania major* structural database. *IJIB*, **7**, 63–68.

49. Patel,P., Mandlik,V. and Singh,S. (2016) LmSmdB: an integrated database for metabolic and gene regulatory network in *Leishmania major* and *Schistosoma mansoni*. *Genom. Data*, **7**, 115–118.

50. Ammari,M.G., Gresham,C.R., McCarthy,F.M. et al. (2016) HPIDB 2.0: a curated database for host-pathogen interactions. *Database*, **2016**, 1–9.

51. Xiang,Z., Tian,Y. and He,Y. (2007) PHIDIAS: a pathogen–host interaction data integration and analysis system. *Genome Biol.*, **8**, 1–15.

52. Urban,M., Cuzick,A., Seager,J. et al. (2022) PHI-base in 2022: a multi-species phenotype database for pathogen–host interactions. *Nucleic Acids Res.*, **50**, D837–D847.

53. Tayal,S., Bhatia,V., Mehrotra,T. et al. (2022) ImitateDB: A database for domain and motif mimicry incorporating host and pathogen protein interactions. *Amino Acids*, **54**, 923–934.

54. Sievers,F., Wilm,A., Dineen,D. et al. (2011) Fast, scalable gener- ation of high-quality protein multiple sequence alignments using Clustal Omega. *Mol. Syst. Biol.*, **7**, 1–6.

55. Tien,M.Z., Meyer,A.G., Sydykova,D.K. et al. (2013) Maximum allowed solvent accessibilities of residues in proteins. *PLoS ONE*, **8**, 1–8.

56. Rost,B. and Sander,C. (1994) Conservation and prediction of sol- vent accessibility in protein families. *Proteins*, **20**, 216–226.

57. Almagro Armenteros,J.J., Sønderby,C.K., Sønderby,S.K. et al. (2017) DeepLoc: prediction of protein subcellular localization using deep learning. *Bioinformatics*, **33**, 3387–3395.

58. Teufel,F., Almagro Armenteros,J.J., Johansen,A.R. et al. (2022) Sig- nalP 6.0 predicts all five types of signal peptides using protein language models. *Nat. Biotechnol.*, **40**, 1023–1025.

59. Gíslason,M.H., Nielsen,H., Armenteros,J.J.A. et al. (2021) Predic- tion of GPI-Anchored proteins with pointer neural networks. *Curr. Res. Biotechnol.*, **3**, 6–13.

60. Varga,J., Dobson,L. and Tusnády,G.E. (2016) TOPDOM: database of conservatively located domains and motifs in proteins. *Bioinfor- matics*, **32**, 2725–2726.

61. Varadi,M., Anyango,S., Appasamy,S.D. et al. (2022) PDBe and PDBe-KB: Providing high-quality, up-to-date and inte- grated resources of macromolecular structures to support basic and applied research and education. *Protein Sci.*, **31**, 1–10.

62. Langó,T., Róna,G., Hunyadi-Gulyás,É. et al. (2017) Identification of extracellular segments by mass spectrometry improves topology prediction of transmembrane proteins. *Sci. Rep.*, **7**, 1–9.

63. Bausch-Fluck,D., Goldmann,U., Müller,S. et al. (2018) The in silico human surfaceome. *Proc. Natl. Acad. Sci.*, **115**, E10988–E10997.

64. Sehnal,D., Bittrich,S., Deshpande,M. et al. (2021) Mol* Viewer: modern web app for 3D visualization and analysis of large biomolecular structures. *Nucleic Acids Res.*, **49**, W431–W437.

65. Tusnády,G.E., Dosztányi,Z. and Simon,I. (2005) TMDET: web server for detecting transmembrane regions of proteins by using their 3D coordinates. *Bioinformatics*, **21**, 1276–1277.