# Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution

Laurine Levillayer, Camille Brighelli, Caroline Demeret, Anavaj Sakuntabhai, Jean-François Bureau

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

## To cite this version:

Laurine Levillayer, Camille Brighelli, Caroline Demeret, Anavaj Sakuntabhai, Jean-François Bureau. Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution. PLoS ONE, 2024, 19 (3), pp.e0296230. 10.1371/journal.pone.0296230. pasteur-04506017

HAL Id: pasteur-04506017

https://pasteur.hal.science/pasteur-04506017

Submitted on 15 Mar 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>
# Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution

Laurine Levillayer, Camille Brighelli, Caroline Demeret, Anavaj Sakuntabhai, Jean-François Bureau

► To cite this version:

Laurine Levillayer, Camille Brighelli, Caroline Demeret, Anavaj Sakuntabhai, Jean-François Bureau. Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution. PLoS ONE, 2024, 19 (3), pp.e0296230. 10.1371/journal.pone.0296230. pasteur-04506017

HAL Id: pasteur-04506017

https://pasteur.hal.science/pasteur-04506017v1

Submitted on 15 Mar 2024

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
RESEARCH ARTICLE

# Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution

Laurine Levillayer<sup>1</sup>, Camille Brighelli<sup>1</sup>, Caroline Demeret<sup>2</sup>, Anavaj Sakuntabhai<sup>1</sup>, Jean-François Bureau<sup>1*</sup>

1 Institut Pasteur, Institut National de Recherche pour l'Agriculture, Université de Paris-Cité, CNRS UMR 2000, l'Alimentation et l'Environnement (INRAE) USC 1510, Unité Écologie et Émergence des Pathogènes Transmis par les Arthropodes (EEPTA), Paris, France, 2 Institut Pasteur, Université de Paris-Cité, Laboratoire Interactomique, ARN et Immunité - Interactomics, RNA and Immunity, Paris, France

* jfb@pasteur.fr

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

## Abstract

SRC kinase associated phosphoprotein 1 (SKAP1), an adaptor for protein assembly, plays an important role in the immune system such as stabilizing immune synapses. Understanding how these functions are controlled at the level of the protein-protein interactions is necessary to describe these processes and to develop therapeutics. Here, we dissected the SKAP1 modular organization to recognize SRC kinases and compared it to that of its paralog SRC kinase associated phosphoprotein 2 (SKAP2). Different conserved motifs common to either both proteins or specific to SKAP2 were found using this comparison. Two modules harboring different binding properties between SKAP1 and SKAP2 were identified: one composed of two conserved motifs located in the second interdomain interacting at least with the SH2 domain of SRC kinases and a second one composed of the DIM domain modulated by the SH3 domain and the activation of SRC kinases. This work suggests a convergent evolution of the binding properties of some SRC kinases interacting specifically with either SKAP1 or SKAP2.

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

### OPEN ACCESS

**Citation:** Levillayer L, Brighelli C, Demeret C, Sakuntabhai A, Bureau J-Franois (2024) Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution. PLoS ONE 19(3): e0296230. https://doi.org/10.1371/journal.pone.0296230

**Editor:** Laszlo Buday, Hungarian Academy of Sciences, HUNGARY

**Received:** December 7, 2023

**Accepted:** February 28, 2024

**Published:** March 14, 2024

**Copyright:** © 2024 Levillayer et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** The author(s) received no specific funding for this work.

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

SRC kinase associated phosphoprotein family is composed of two members, SKAP1 and SKAP2. Due to their localization to human chromosomes 7p15.2 for SKAP2 and 17q21.32 for SKAP1 at one extremity to respectively HoxA and HoxB cluster, their duplication occurs during evolution of chordates to vertebrates [1]. At least since gnathostome radiation, SKAP1 and SKAP2 are highly conserved. SKAP1 is mainly expressed in the immune system particularly in T-cells [2,3] in contrast to SKAP2 that is widely expressed [4] even if at higher level in some cells of the immune system such as macrophages and neutrophils [5-8]. SKAP2 also plays an important role during the development of the Central Nervous System [9,10]. The structure of
RESEARCH ARTICLE

# Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution

Laurine Levillayer<sup>1</sup>, Camille Brighelli<sup>1</sup>, Caroline Demeret<sup>2</sup>, Anavaj Sakuntabhai<sup>1</sup>, Jean-François Bureau<sup>1*</sup>

1 Institut Pasteur, Institut National de Recherche pour l'Agriculture, Université de Paris-Cité, CNRS UMR 2000, l'Alimentation et l'Environnement (INRAE) USC 1510, Unité Écologie et Émergence des Pathogènes Transmis par les Arthropodes (EEPTA), Paris, France, 2 Institut Pasteur, Université de Paris-Cité, Laboratoire Interactomique, ARN et Immunité - Interactomics, RNA and Immunity, Paris, France

* jfb@pasteur.fr

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

## Abstract

SRC kinase associated phosphoprotein 1 (SKAP1), an adaptor for protein assembly, plays an important role in the immune system such as stabilizing immune synapses. Understanding how these functions are controlled at the level of the protein-protein interactions is necessary to describe these processes and to develop therapeutics. Here, we dissected the SKAP1 modular organization to recognize SRC kinases and compared it to that of its paralog SRC kinase associated phosphoprotein 2 (SKAP2). Different conserved motifs common to either both proteins or specific to SKAP2 were found using this comparison. Two modules harboring different binding properties between SKAP1 and SKAP2 were identified: one composed of two conserved motifs located in the second interdomain interacting at least with the SH2 domain of SRC kinases and a second one composed of the DIM domain modulated by the SH3 domain and the activation of SRC kinases. This work suggests a convergent evolution of the binding properties of some SRC kinases interacting specifically with either SKAP1 or SKAP2.

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

### OPEN ACCESS

**Citation:** Levillayer L, Brighelli C, Demeret C, Sakuntabhai A, Bureau J-Franois (2024) Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution. PLoS ONE 19(3): e0296230. https://doi.org/10.1371/journal.pone.0296230

**Editor:** Laszlo Buday, Hungarian Academy of Sciences, HUNGARY

**Received:** December 7, 2023

**Accepted:** February 28, 2024

**Published:** March 14, 2024

**Copyright:** © 2024 Levillayer et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** The author(s) received no specific funding for this work.

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

SRC kinase associated phosphoprotein family is composed of two members, SKAP1 and SKAP2. Due to their localization to human chromosomes 7p15.2 for SKAP2 and 17q21.32 for SKAP1 at one extremity to respectively HoxA and HoxB cluster, their duplication occurs during evolution of chordates to vertebrates [1]. At least since gnathostome radiation, SKAP1 and SKAP2 are highly conserved. SKAP1 is mainly expressed in the immune system particularly in T-cells [2,3] in contrast to SKAP2 that is widely expressed [4] even if at higher level in some cells of the immune system such as macrophages and neutrophils [5-8]. SKAP2 also plays an important role during the development of the Central Nervous System [9,10]. The structure of
both proteins is composed of a dimerization domain (DIM) followed by a first interdomain (Link1), a Pleckstrin homology domain (PH), a second interdomain (Link2), and a SRC homology 3 domain (SH3) (Fig 1). The homologous or heterologous association of two DIM domains by hydrophobic interactions formed a four-helix bundle structure that has been solved for SKAP2 [11,12]. Their PH domain interacts with membranes via two phosphate-derivatives of phosphatidyl inositol, PI(3,4,5)P₃ and either PI(4,5)P₂ for SKAP1 [13] or PI(3,4)P₂ for SKAP2 [12]. Interestingly, the DIM and the PH domains of the same chain interact together with an important role of D120 and D129, respectively for SKAP1 and SKAP2. The interdomain Link2 contains two conserved tyrosine phosphorylated motifs that have been suspected to bind to the SH2 domain of SRC kinases [14]. Both SH3 domains bind proline-rich motif of PRAM1 [15], FYB1 [4,16], FYB2 and FAM102A [17,18], in which W333 and W336 play a necessary role for SKAP1 and SKAP2, respectively.

At the T-cell immune synapse during the inside-out signaling, SKAP1 stabilizes SLP-76/LCP2 microclusters by interacting with FYB1 [19,20]. Interestingly, the interaction of either SKAP1 or SKAP2 with FYB1 protects them from degradation [21]. SKAP1 also formed a complex with LAT and GRB2 induced by the integrin LFA-1 and PTK2A/PTK2B, which terminate T-cell conjugate formation [20]. During this process, SKAP1 is recruited to lipid raft, in which the SRC kinase FYN phosphorylates its tyrosine 271, a necessary event for its interaction with GRB2 [14]. Note also that Skap1 interacts by its DIM domain with Rapl and RIAM, regulators of respectively α and β-integrin [19,22,23]. SKAP1 was also reported to interact with Polo-like kinase 1 for an optimal control of cell-cycling [24]. A recent review [25] discuss these data in more details. The FYB1-SKAP1 complex has also been detected in mast cells but its function is unknown [26]. In conclusion, both SKAP1 and SKAP2 play an important role as adaptors for the same proteins such as FYB1, FAM102A, and SRC kinases. Our recent work revealed that the binding of SRC kinases by SKAP2 is dictated by its complex modular architecture, inviting to study the SKAP1 one [27].

The SRC kinase family is composed of the SrcA subfamily with FGR, FYN, SRC, and YES, the SrcB subfamily with BLK, HCK, LCK, and LYN and finally the FRK-related family with FRK, SRMS, and PTK6 [28,29]. The composition of the FRK-related family is still being debated. SRC kinases [30], which phosphorylate tyrosine residues, are composed of three domains, a SRC homology 3 (SH3) domain, a SRC homology 2 (SH2) domain, and a tyrosine kinase catalytic domain (TYRKC) (Fig 1D). A fine allosteric control of the kinase activity is modulated from the passage of a closed to open conformation [31–33] (Fig 1C). The closed, inactive conformation is provided by intramolecular interactions between (1) the SH3 domain and a proline-rich motif located in the hinge between the SH2 and the catalytic kinase domains and (2) the SH2 domain and a phosphorylated tyrosine motif in the C-terminus of the kinase. In the open active conformation, these intramolecular interactions are lost, and a transphosphorylation of a tyrosine located in the active loop of the catalytic site locks the system [34]. The localization of most SRC kinases in different membranous compartments is finely regulated by N-terminal myristylation and palmitoylation [35,36]. However, neither SRMS nor FRK have an N-terminal myristylation site.

Evolution of protein-protein interactions has been a subject of rich debates on the gain and loss of interactors particularly after a duplication event [37–42]. Alternative splicing also plays an important role to explain tissue-specificity of an interaction and of monogenic disorders [43–45]. Also, evolution information was frequently used to predict interaction [38]. Large protein and domain families have been particularly studied [46]. Kinases with an SRC module constitute one family in which such analysis was performed [47–49]. This kind of analysis was also conducted on domains such as the SH2 and the SH3 domains [50–53]. One interesting case was protein-protein interactions between viral proteins and host ones [54]. Most of the
both proteins is composed of a dimerization domain (DIM) followed by a first interdomain (Link1), a Pleckstrin homology domain (PH), a second interdomain (Link2), and a SRC homology 3 domain (SH3) (Fig 1). The homologous or heterologous association of two DIM domains by hydrophobic interactions formed a four-helix bundle structure that has been solved for SKAP2 [11,12]. Their PH domain interacts with membranes via two phosphate-derivatives of phosphatidyl inositol, PI(3,4,5)P₃ and either PI(4,5)P₂ for SKAP1 [13] or PI(3,4)P₂ for SKAP2 [12]. Interestingly, the DIM and the PH domains of the same chain interact together with an important role of D120 and D129, respectively for SKAP1 and SKAP2. The interdomain Link2 contains two conserved tyrosine phosphorylated motifs that have been suspected to bind to the SH2 domain of SRC kinases [14]. Both SH3 domains bind proline-rich motif of PRAM1 [15], FYB1 [4,16], FYB2 and FAM102A [17,18], in which W333 and W336 play a necessary role for SKAP1 and SKAP2, respectively.

At the T-cell immune synapse during the inside-out signaling, SKAP1 stabilizes SLP-76/LCP2 microclusters by interacting with FYB1 [19,20]. Interestingly, the interaction of either SKAP1 or SKAP2 with FYB1 protects them from degradation [21]. SKAP1 also formed a complex with LAT and GRB2 induced by the integrin LFA-1 and PTK2A/PTK2B, which terminate T-cell conjugate formation [20]. During this process, SKAP1 is recruited to lipid raft, in which the SRC kinase FYN phosphorylates its tyrosine 271, a necessary event for its interaction with GRB2 [14]. Note also that Skap1 interacts by its DIM domain with Rapl and RIAM, regulators of respectively α and β-integrin [19,22,23]. SKAP1 was also reported to interact with Polo-like kinase 1 for an optimal control of cell-cycling [24]. A recent review [25] discuss these data in more details. The FYB1-SKAP1 complex has also been detected in mast cells but its function is unknown [26]. In conclusion, both SKAP1 and SKAP2 play an important role as adaptors for the same proteins such as FYB1, FAM102A, and SRC kinases. Our recent work revealed that the binding of SRC kinases by SKAP2 is dictated by its complex modular architecture, inviting to study the SKAP1 one [27].

The SRC kinase family is composed of the SrcA subfamily with FGR, FYN, SRC, and YES, the SrcB subfamily with BLK, HCK, LCK, and LYN and finally the FRK-related family with FRK, SRMS, and PTK6 [28,29]. The composition of the FRK-related family is still being debated. SRC kinases [30], which phosphorylate tyrosine residues, are composed of three domains, a SRC homology 3 (SH3) domain, a SRC homology 2 (SH2) domain, and a tyrosine kinase catalytic domain (TYRKC) (Fig 1D). A fine allosteric control of the kinase activity is modulated from the passage of a closed to open conformation [31–33] (Fig 1C). The closed, inactive conformation is provided by intramolecular interactions between (1) the SH3 domain and a proline-rich motif located in the hinge between the SH2 and the catalytic kinase domains and (2) the SH2 domain and a phosphorylated tyrosine motif in the C-terminus of the kinase. In the open active conformation, these intramolecular interactions are lost, and a transphosphorylation of a tyrosine located in the active loop of the catalytic site locks the system [34]. The localization of most SRC kinases in different membranous compartments is finely regulated by N-terminal myristylation and palmitoylation [35,36]. However, neither SRMS nor FRK have an N-terminal myristylation site.

Evolution of protein-protein interactions has been a subject of rich debates on the gain and loss of interactors particularly after a duplication event [37–42]. Alternative splicing also plays an important role to explain tissue-specificity of an interaction and of monogenic disorders [43–45]. Also, evolution information was frequently used to predict interaction [38]. Large protein and domain families have been particularly studied [46]. Kinases with an SRC module constitute one family in which such analysis was performed [47–49]. This kind of analysis was also conducted on domains such as the SH2 and the SH3 domains [50–53]. One interesting case was protein-protein interactions between viral proteins and host ones [54]. Most of the
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Map of SKAP1, SKAP2, and HCK mutants. For each protein, the delineation of domain boundaries is indicated according to Uniprot database. Wild type and replacement amino acid identities at mutated position involved in domain or motif inactivation are indicated in red and blue respectively. (A) Schematic of SKAP1 organization with the three domains, domain of dimerization (DIM), Pleckstrin homology domain (PH), SRC homology 3 domain (SH3), and the two interdomains, Link1 and Link2 with deletion mutants. (B) Schematic of SKAP2 organization with the three domains, DIM, PH, and SH3. (C) Schematic illustrating inactive (bottom) and active (top) forms of HCK. (D) Schematic of HCK protein organization with its three domains: SRC homology 3 domain (SH3), SRC homology 2 domain (SH2) and tyrosine kinase catalytic domain (TYRC).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g001

works were done at the global level of the interactome but some of them also analyzed the level of protein-protein interfaces [38,55,56]. Knowledge at this level is necessary to develop protein-protein inhibitors [57,58]. Specific cellular expression occurs among some SRC kinases as for SKAP1 and SKAP2 allowing to test possible convergent evolution. The FYN kinase is particularly interesting because one isoform, FYNT, is mainly expressed in T lymphocytes as SKAP1, and a second one, FYNB, is expressed in the brain as SKAP2 [59]. Both isoforms differ by their exon 7 sequence, which were absent to the FYN isoform that we used in our previous work [60,61]. Interestingly, mis-expression of these two isoforms in T lymphocytes as in the brain has been associated to an increasing risk for some diseases [62]. LCK kinase is also mainly expressed in T lymphocytes in contrast to HCK kinase, which is expressed in the monocytes/macrophages lineage and neutrophils as SKAP2.

In the present paper, we described SKAP1 architecture recognizing SRC kinases, and compared it to that of its paralog SKAP2 [27]. Comparison of similar deletion mutants and mutations between these two adaptors allowed to precisely define conserved and divergent functions. We identified a motif in Link1 interdomain, specific to SKAP2, which affects the interactions with SRC kinases. The conservation of the interdomain Link2 function between both proteins allowed the detection of two common conserved motifs, which play non-redundant and possibly cooperative roles in the binding to the SH2 domain of SRC kinases. These two conserved motifs and at least the SH2 domain of SRC kinases constitute the first interaction unit identified; the second one is composed of the DIM domain, modulated by its SH3 domain and the activated SRC kinases. These two modules in SKAP1 and SKAP2 harbor different binding properties to SRC kinases. Interestingly, LCK and FYNT kinases expressed in T lymphocytes as SKAP1 show different binding properties than HCK kinase expressed in the monocytes/macrophages lineage and neutrophils, and FYNB kinase expressed mainly in the brain. SKAP2 is also
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Map of SKAP1, SKAP2, and HCK mutants. For each protein, the delineation of domain boundaries is indicated according to Uniprot database. Wild type and replacement amino acid identities at mutated position involved in domain or motif inactivation are indicated in red and blue respectively. (A) Schematic of SKAP1 organization with the three domains, domain of dimerization (DIM), Pleckstrin homology domain (PH), SRC homology 3 domain (SH3), and the two interdomains, Link1 and Link2 with deletion mutants. (B) Schematic of SKAP2 organization with the three domains, DIM, PH, and SH3. (C) Schematic illustrating inactive (bottom) and active (top) forms of HCK. (D) Schematic of HCK protein organization with its three domains: SRC homology 3 domain (SH3), SRC homology 2 domain (SH2) and tyrosine kinase catalytic domain (TYRC).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g001

works were done at the global level of the interactome but some of them also analyzed the level of protein-protein interfaces [38,55,56]. Knowledge at this level is necessary to develop protein-protein inhibitors [57,58]. Specific cellular expression occurs among some SRC kinases as for SKAP1 and SKAP2 allowing to test possible convergent evolution. The FYN kinase is particularly interesting because one isoform, FYNT, is mainly expressed in T lymphocytes as SKAP1, and a second one, FYNB, is expressed in the brain as SKAP2 [59]. Both isoforms differ by their exon 7 sequence, which were absent to the FYN isoform that we used in our previous work [60,61]. Interestingly, mis-expression of these two isoforms in T lymphocytes as in the brain has been associated to an increasing risk for some diseases [62]. LCK kinase is also mainly expressed in T lymphocytes in contrast to HCK kinase, which is expressed in the monocytes/macrophages lineage and neutrophils as SKAP2.

In the present paper, we described SKAP1 architecture recognizing SRC kinases, and compared it to that of its paralog SKAP2 [27]. Comparison of similar deletion mutants and mutations between these two adaptors allowed to precisely define conserved and divergent functions. We identified a motif in Link1 interdomain, specific to SKAP2, which affects the interactions with SRC kinases. The conservation of the interdomain Link2 function between both proteins allowed the detection of two common conserved motifs, which play non-redundant and possibly cooperative roles in the binding to the SH2 domain of SRC kinases. These two conserved motifs and at least the SH2 domain of SRC kinases constitute the first interaction unit identified; the second one is composed of the DIM domain, modulated by its SH3 domain and the activated SRC kinases. These two modules in SKAP1 and SKAP2 harbor different binding properties to SRC kinases. Interestingly, LCK and FYNT kinases expressed in T lymphocytes as SKAP1 show different binding properties than HCK kinase expressed in the monocytes/macrophages lineage and neutrophils, and FYNB kinase expressed mainly in the brain. SKAP2 is also
mainly expressed in the monocytes/macrophages lineage, in neutrophils and in the brain suggesting a possible coevolution between these two adaptors and SRC kinases.

## Materials and methods

### Plasmids, mutagenesis, BP and LR cloning, and luciferase complementation assay

Previously used plasmids and a more detailed methodology will be found elsewhere [17,27]. The ORF of FYNT (NM153047.4), FYNB (NM002037.5) flanked by two gateway sites have been synthesized and cloned in PUC57 plasmid (GenScript, Hong Kong). These two clones were individually transferred by Gateway recombinational cloning into pDONR207 vector. The ORF of the new SKAP2 deletion mutants flanked by two gateway sites were amplified by PCR, cloned into pCR-II TOPO vector using TOPO TA cloning kit (Thermo Fisher Scientific, Waltham) and transferred into pDONR207 vector. Mutagenesis of SKAP1, SKAP2, and HCK ORFs were performed using QuickChange Lightening Site-Directed Mutagenesis kit (Agilent Technologies, Les Ulis) according to the manufacturer's protocol. S1 Table in S2 File shows the new primers used for all these purposes. The sequence of other primers will be found in previous publications [17,27]. After each mutagenesis and PCR amplification, the insert was completely sequenced at Eurofin services using LightRun Barcodes (Eurofins Genomics, Ebersberg, Germany). Sequence alignments were performed using DNA Strider [63]. The ORF of SKAP1, and SKAP2 mutants were transferred into pSPICA-N2 vector and those of SRC kinases and HCK mutants into pSPICA-N1 and pSPICA-C1 ones. The pSPICA vectors are mammalian expressing vectors designed for luciferase complementation assay [64]. They expressed *Gaussia princeps* Luciferase fragment 1 (amino acid residues from 18 to 109) or 2 (amino acid residues from 110 to 185) at the N or C-terminal part of the fusion protein. The luciferase complementation assay was performed as described in human HEK293T cells [17,64]. Briefly the first day, 30 000 HEK293T cells are cultured in 100 µl of DMEM supplemented with 10% Fetal Bovine Serum and antibiotics per well of microplate 96 wells (Greiner, Kremsmünster). One day later, cells were transfected with 100 ng GPCA plasmid pair using polyethylenimine (PEI) method. The third day, culture medium was discarded, cells were washed once with 150 µl of Phosphate Buffer Saline without calcium and magnesium (PBS) and incubated for 25 min in 40 µl of lysis buffer. Luminescence monitoring was performed after addition of native Coelenterazine on a Centro XS³ LB 960 microplate luminometer (Berthold Technologies, Thoiry). Transfections using PEI were performed at least in triplicate. Protein-protein interactions were monitored by measuring interaction-mediated normalized luminescence ratio (NLR) [64]. Briefly, NLR is the ratio between the interaction signal and a control one, the sum of signals of each component of the interaction with the complementary empty hemi-luciferase expressing vector.

### Experimental design and statistical rationale

Approximation of the mean and the variance of a ratio were performed using Taylor expansions with a null covariance between the numerator and the denominator. Two-sample z test and two-way ANOVA after log-normal transformation were used to compare these ratios. Differential interaction scatterplots and PPI-mutation plots were performed as described in [17]. Three levels of significance were used for the P-value lower than 0.05 (*), 0.01 (**), and 0.001 (***) . Bonferroni correction was used for multiple testing with an initial threshold p = 0.05. To compare different experiments, we used modified NLR for each SRC kinases as the ratio between NLR of the SKAP1 mutant and NLR of wild-type SKAP1 protein, which was named
mainly expressed in the monocytes/macrophages lineage, in neutrophils and in the brain suggesting a possible coevolution between these two adaptors and SRC kinases.

## Materials and methods

### Plasmids, mutagenesis, BP and LR cloning, and luciferase complementation assay

Previously used plasmids and a more detailed methodology will be found elsewhere [17,27]. The ORF of *FYNT* (NM153047.4), *FYNB* (NM002037.5) flanked by two gateway sites have been synthesized and cloned in PUC57 plasmid (GenScript, Hong Kong). These two clones were individually transferred by Gateway recombinational cloning into pDONR207 vector. The ORF of the new SKAP2 deletion mutants flanked by two gateway sites were amplified by PCR, cloned into pCR-II TOPO vector using TOPO TA cloning kit (Thermo Fisher Scientific, Waltham) and transferred into pDONR207 vector. Mutagenesis of SKAP1, SKAP2, and HCK ORFs were performed using QuickChange Lightening Site-Directed Mutagenesis kit (Agilent Technologies, Les Ulis) according to the manufacturer's protocol. S1 Table in S2 File shows the new primers used for all these purposes. The sequence of other primers will be found in previous publications [17,27]. After each mutagenesis and PCR amplification, the insert was completely sequenced at Eurofin services using LightRun Barcodes (Eurofins Genomics, Ebersberg, Germany). Sequence alignments were performed using DNA Strider [63]. The ORF of SKAP1, and SKAP2 mutants were transferred into pSPICA-N2 vector and those of SRC kinases and HCK mutants into pSPICA-N1 and pSPICA-C1 ones. The pSPICA vectors are mammalian expressing vectors designed for luciferase complementation assay [64]. They expressed *Gaussia princeps* Luciferase fragment 1 (amino acid residues from 18 to 109) or 2 (amino acid residues from 110 to 185) at the N or C-terminal part of the fusion protein. The luciferase complementation assay was performed as described in human HEK293T cells [17,64]. Briefly the first day, 30 000 HEK293T cells are cultured in 100 µl of DMEM supplemented with 10% Fetal Bovine Serum and antibiotics per well of microplate 96 wells (Greiner, Kremsmünster). One day later, cells were transfected with 100 ng GPCA plasmid pair using polyethylenimine (PEI) method. The third day, culture medium was discarded, cells were washed once with 150 µl of Phosphate Buffer Saline without calcium and magnesium (PBS) and incubated for 25 min in 40 µl of lysis buffer. Luminescence monitoring was performed after addition of native Coelenterazine on a Centro XS³ LB 960 microplate luminometer (Berthold Technologies, Thoiry). Transfections using PEI were performed at least in triplicate. Protein-protein interactions were monitored by measuring interaction-mediated normalized luminescence ratio (NLR) [64]. Briefly, NLR is the ratio between the interaction signal and a control one, the sum of signals of each component of the interaction with the complementary empty hemi-luciferase expressing vector.

### Experimental design and statistical rationale

Approximation of the mean and the variance of a ratio were performed using Taylor expansions with a null covariance between the numerator and the denominator. Two-sample z test and two-way ANOVA after log-normal transformation were used to compare these ratios. Differential interaction scatterplots and PPI-mutation plots were performed as described in [17]. Three levels of significance were used for the P-value lower than 0.05 (*), 0.01 (**), and 0.001 (***) . Bonferroni correction was used for multiple testing with an initial threshold p = 0.05. To compare different experiments, we used modified NLR for each SRC kinases as the ratio between NLR of the SKAP1 mutant and NLR of wild-type SKAP1 protein, which was named
wild-type normalized NLR. Experiments were repeated three times. For studying the effect of the different domains and mutations, we used the ratio between NLR of the SKAP1 mutant and NLR of the same SKAP2 mutant, which was named SKAP2 normalized NLR. Study of each mutant experiment was repeated at least two times. Statistical analyses were performed using Stata/IC version 10.1 (Stata Corporation, College Station, USA). Comparison of mean to theoretical value (i.e. Mean ratio to 1 and its log-normal transformation to 0) and between two means was performed by using Student t test and unpaired Student t test, respectively. Comparison of mean to theoretical value was also studied using binomial distribution after calculating the number of samples upper/lower to theoretical value. Comparison of two means was also performed by Kruskal-Wallis test. Parametric and non-parametric tests giving similar results, only results of the parametric are shown in Table 1. We also decided to present preferentially ratio data that are easier to interpret than those after log-normal transformation when both analyses gave similar results. The only exception was the two-way ANOVA, which was performed after a log-normal transformation of the modified NLR and using analytic weights. For data modelling, the full model contains six explanatory variables and the data on the N1- and C1-fused SRC kinases were analyzed separately. Six different explicative variables were generated from the initial experiments on the different mutants. An explicative model was generated from this full model by a step-by-step procedure that removed at each step the lowest non-significant variable until no non-significant variables were present. A principal component analysis (PCA) from the present mutant data and those of SKAP2 was performed after a log-transformation using the R package Factominer [65] and Factoshiny according to the author specification to determine which factors explained the variability amongst C1-fused SRC kinases. More details on the methodology, complete results, and interpretation are shown in supplemental data.

## Conservation of the Casein kinase 1 motif in the interdomain Link1 of SKAP2

The N-terminal fragment of human SKAP2 protein structure was predicted from I-TASSER server [66] using two crystal structures of a N-terminal fragment of murine SKAP2 containing both the helical dimerization domain and the PH domain (PDB ID: 2OTX and 1U5E). Clustalw program was used to align respectively 38 protein sequences of SKAP2 from the Coelacanth to mammals (NCBI database, S2 Table in S2 File). ConSurf server [67,68] was used to visualize evolutionary conservation in this macromolecule from its predict structure and aligned sequences using a maximum likelihood procedure. UCSF Chimera version 1.13.1 [69] was used to visualize the conservation of these structures using a script from ConSurf server.

## Results

In this paper, we study the SKAP1 organization, which recognizes the SRC kinases. The goal of this work was 1) to define the SKAP1 modular architecture of this interaction; 2) to specify the composition and binding properties of its two modules by looking for common or specific motifs to SKAP1 and SKAP2 architecture; 3) to compare the two architectures in the context of evolution. For that purpose, we used a luciferase complementation assay in HEK293T cells where the C-terminal hemi-gaussia luciferase (Gluc2) is fused at the N terminus of SKAP1 (N2-SKAP1) and the N-terminal hemi-gaussia luciferase (Gluc1) is fused to either the N terminus (N1) or C terminus (C1) of the SRC kinases. The two orientations of hemi-gaussian luciferase fused to SRC kinases were studied because we previously showed that they affect their subcellular localization and activation but also the interactions with SKAP2 [27]. By similarity of what we found for SKAP2 adaptor [27], our assumption was that there is a general SKAP1
wild-type normalized NLR. Experiments were repeated three times. For studying the effect of the different domains and mutations, we used the ratio between NLR of the SKAP1 mutant and NLR of the same SKAP2 mutant, which was named SKAP2 normalized NLR. Study of each mutant experiment was repeated at least two times. Statistical analyses were performed using Stata/IC version 10.1 (Stata Corporation, College Station, USA). Comparison of mean to theoretical value (i.e. Mean ratio to 1 and its log-normal transformation to 0) and between two means was performed by using Student t test and unpaired Student t test, respectively. Comparison of mean to theoretical value was also studied using binomial distribution after calculating the number of samples upper/lower to theoretical value. Comparison of two means was also performed by Kruskal-Wallis test. Parametric and non-parametric tests giving similar results, only results of the parametric are shown in Table 1. We also decided to present preferentially ratio data that are easier to interpret than those after log-normal transformation when both analyses gave similar results. The only exception was the two-way ANOVA, which was performed after a log-normal transformation of the modified NLR and using analytic weights. For data modelling, the full model contains six explanatory variables and the data on the N1- and C1-fused SRC kinases were analyzed separately. Six different explicative variables were generated from the initial experiments on the different mutants. An explicative model was generated from this full model by a step-by-step procedure that removed at each step the lowest non-significant variable until no non-significant variables were present. A principal component analysis (PCA) from the present mutant data and those of SKAP2 was performed after a log-transformation using the R package Factominer [65] and Factoshiny according to the author specification to determine which factors explained the variability amongst C1-fused SRC kinases. More details on the methodology, complete results, and interpretation are shown in supplemental data.

## Conservation of the Casein kinase 1 motif in the interdomain Link1 of SKAP2

The N-terminal fragment of human SKAP2 protein structure was predicted from I-TASSER server [66] using two crystal structures of a N-terminal fragment of murine SKAP2 containing both the helical dimerization domain and the PH domain (PDB ID: 2OTX and 1U5E). Clustalw program was used to align respectively 38 protein sequences of SKAP2 from the Coelacanth to mammals (NCBI database, S2 Table in S2 File). ConSurf server [67,68] was used to visualize evolutionary conservation in this macromolecule from its predict structure and aligned sequences using a maximum likelihood procedure. UCSF Chimera version 1.13.1 [69] was used to visualize the conservation of these structures using a script from ConSurf server.

## Results

In this paper, we study the SKAP1 organization, which recognizes the SRC kinases. The goal of this work was 1) to define the SKAP1 modular architecture of this interaction; 2) to specify the composition and binding properties of its two modules by looking for common or specific motifs to SKAP1 and SKAP2 architecture; 3) to compare the two architectures in the context of evolution. For that purpose, we used a luciferase complementation assay in HEK293T cells where the C-terminal hemi-gaussia luciferase (Gluc2) is fused at the N terminus of SKAP1 (N2-SKAP1) and the N-terminal hemi-gaussia luciferase (Gluc1) is fused to either the N terminus (N1) or C terminus (C1) of the SRC kinases. The two orientations of hemi-gaussian luciferase fused to SRC kinases were studied because we previously showed that they affect their subcellular localization and activation but also the interactions with SKAP2 [27]. By similarity of what we found for SKAP2 adaptor [27], our assumption was that there is a general SKAP1
<table><caption>Table 1. Statistical analysis of the interaction between SKAP1 mutants and SRC kinases.</caption><thead><tr><th>Mutant</th><th>Nb Exp (n)</th><th>N1 effect*</th><th>C1 effect*</th><th>ΔC1-N1*</th></tr></thead><tbody><tr><td>D120K</td><td>3 (27)</td><td><strong>1.17 ± 0.05</strong> **</td><td><strong>1.49 ± 0.08</strong> ***</td><td><strong>0.32 ± 0.09</strong> **</td></tr><tr><td>DIM-</td><td>3 (27)</td><td>0.86 ± 0.06 *</td><td>1.10 ± 0.14</td><td>0.23 ± 0.15</td></tr><tr><td>ΔLink1PH</td><td>3 (27)</td><td><strong>1.57 ± 0.11</strong> ***</td><td>0.96 ± 0.07</td><td><strong>-0.60 ± 0.12</strong> ***</td></tr><tr><td>ΔPH</td><td>3 (27)</td><td><strong>1.35 ± 0.09</strong> ***</td><td>0.87 ± 0.07 *</td><td><strong>-0.48 ± 0.11</strong> ***</td></tr><tr><td>ΔPHLink2</td><td>3 (27)</td><td><strong>0.40 ± 0.04</strong> ***</td><td><strong>0.30 ± 0.03</strong> ***</td><td>-0.09 ± 0.05 *</td></tr><tr><td>ΔSH3</td><td>3 (27)</td><td>1.51 ± 0.20 **</td><td><strong>2.53 ± 0.37</strong> ***</td><td>1.02 ± 0.42 **</td></tr><tr><td>DomDIM</td><td>3 (27)</td><td><strong>0.11 ± 0.01</strong> ***</td><td><strong>0.32 ± 0.04</strong> ***</td><td><strong>0.22 ± 0.04</strong> ***</td></tr><tr><td>Link2SH3</td><td>3 (27)</td><td><strong>0.29 ± 0.02</strong> ***</td><td><strong>0.51 ± 0.05</strong> ***</td><td><strong>0.22 ± 0.05</strong> ***</td></tr><tr><td>W333K</td><td>3 (27)</td><td>1.16 ± 0.05 **</td><td><strong>1.24 ± 0.06</strong> ***</td><td>0.08 ± 0.08</td></tr></tbody></table>

architecture, which defines its interaction with all the SRC kinases, and was analyzed using the mean of SRC kinase PPI signals. We measured each PPI signal by using the normalized luciferase ratio [64] and the mean value between SRC kinases of each mutant to detect the general architecture. In practice, to assess the effect of SKAP1 deletions or mutations, the NLR generated by their binding to an SRC kinase has been divided by the binding of wildtype SKAP1 to the same kinase. This generated a wild-type normalized NLR enabling to standardize the interaction of mutated SKAP1 for each SRC kinase amongst the different experiments. The SKAP1 interaction with nine SRC kinases, BLK, FGR, FRK, FYN, HCK, LCK, LYN, SRMS, and YES, was studied with 12 SKAP1 mutants, nine of which are deletion mutants (Fig 1A) and the others contain nonsynonymous mutations. To compare architecture between SKAP1 and SKAP2, we also studied similar SKAP2 mutations for their interaction with the nine SRC kinases. Fig 1B shows the localization of all SKAP2 nonsynonymous mutations. Lastly, we studied the specificity of some SRC kinase interactions with SKAP1 and SKAP2 by the analysis of two nonsynonymous mutations on HCK (Fig 1D) or LCK background and of the three FYN isoforms, and more systematically by using a Principal Component Analysis (PCA) on C1-fused SRC kinases and either SKAP1 or SKAP2 mutant data.

## SKAP1 modular architecture interacting with SRC kinases

We first studied the interaction of nine SKAP1 mutants with nine SRC kinases, which has been able to previously define the SKAP2 architecture [27]. Fig 2 shows one experiment for each SKAP1 mutant. The interaction of SRC kinases systematically decreases for ΔPHLink2 (deletion from aa 107 to 326), DomDIM1 (sequence from aa 1 to 99), and Link2SH3 (sequence from aa 213 to 358) SKAP1 mutants, independently of the localization of the hemi-gaussia fusion. This strongly suggests that the domain DIM and the interdomain Link2 are necessary and sufficient for the interaction with the SRC kinases. The Link2SH3 and the DomDIM mutants have a similar loss of signal indicating that neither the interdomain Link2 nor the DIM domain are sufficient. The ΔPH (deletion from aa 107 to 210), the ΔLink1PH (deletion from aa 67 to 210), and the ΔSH3 (sequence from aa 1 to 291) mutants do not exhibit such decrease of the PPI signal supporting that they are not necessary. The W333K SKAP1 mutant, invalidated for the binding to the proline-rich binding motif, behave similarly as the ΔSH3 mutant, showing a slight increase in their binding capacity. The D120K and DIM* mutations have a slight effect on the PPI respectively increasing the signal for C1-fused SRC kinases and decreasing it for N1-fused ones.

<footnote>* Mean of normalized NLR ± SEM p = 0.05 (*); p = 0.01 (**); p = 0.001 (***); bold: p = 0.05 after Bonferroni correction.</footnote>
<table><caption>Table 1. Statistical analysis of the interaction between SKAP1 mutants and SRC kinases.</caption><thead><tr><th>Mutant</th><th>Nb Exp (n)</th><th>N1 effect*</th><th>C1 effect*</th><th>ΔC1-N1*</th></tr></thead><tbody><tr><td>D120K</td><td>3 (27)</td><td><strong>1.17 ± 0.05</strong> **</td><td><strong>1.49 ± 0.08</strong> ***</td><td><strong>0.32 ± 0.09</strong> **</td></tr><tr><td>DIM-</td><td>3 (27)</td><td>0.86 ± 0.06 *</td><td>1.10 ± 0.14</td><td>0.23 ± 0.15</td></tr><tr><td>ΔLink1PH</td><td>3 (27)</td><td><strong>1.57 ± 0.11</strong> ***</td><td>0.96 ± 0.07</td><td><strong>-0.60 ± 0.12</strong> ***</td></tr><tr><td>ΔPH</td><td>3 (27)</td><td><strong>1.35 ± 0.09</strong> ***</td><td>0.87 ± 0.07 *</td><td><strong>-0.48 ± 0.11</strong> ***</td></tr><tr><td>ΔPHLink2</td><td>3 (27)</td><td><strong>0.40 ± 0.04</strong> ***</td><td><strong>0.30 ± 0.03</strong> ***</td><td>-0.09 ± 0.05 *</td></tr><tr><td>ΔSH3</td><td>3 (27)</td><td>1.51 ± 0.20 **</td><td><strong>2.53 ± 0.37</strong> ***</td><td>1.02 ± 0.42 **</td></tr><tr><td>DomDIM</td><td>3 (27)</td><td><strong>0.11 ± 0.01</strong> ***</td><td><strong>0.32 ± 0.04</strong> ***</td><td><strong>0.22 ± 0.04</strong> ***</td></tr><tr><td>Link2SH3</td><td>3 (27)</td><td><strong>0.29 ± 0.02</strong> ***</td><td><strong>0.51 ± 0.05</strong> ***</td><td><strong>0.22 ± 0.05</strong> ***</td></tr><tr><td>W333K</td><td>3 (27)</td><td>1.16 ± 0.05 **</td><td><strong>1.24 ± 0.06</strong> ***</td><td>0.08 ± 0.08</td></tr></tbody></table>

architecture, which defines its interaction with all the SRC kinases, and was analyzed using the mean of SRC kinase PPI signals. We measured each PPI signal by using the normalized luciferase ratio [64] and the mean value between SRC kinases of each mutant to detect the general architecture. In practice, to assess the effect of SKAP1 deletions or mutations, the NLR generated by their binding to an SRC kinase has been divided by the binding of wildtype SKAP1 to the same kinase. This generated a wild-type normalized NLR enabling to standardize the interaction of mutated SKAP1 for each SRC kinase amongst the different experiments. The SKAP1 interaction with nine SRC kinases, BLK, FGR, FRK, FYN, HCK, LCK, LYN, SRMS, and YES, was studied with 12 SKAP1 mutants, nine of which are deletion mutants (Fig 1A) and the others contain nonsynonymous mutations. To compare architecture between SKAP1 and SKAP2, we also studied similar SKAP2 mutations for their interaction with the nine SRC kinases. Fig 1B shows the localization of all SKAP2 nonsynonymous mutations. Lastly, we studied the specificity of some SRC kinase interactions with SKAP1 and SKAP2 by the analysis of two nonsynonymous mutations on HCK (Fig 1D) or LCK background and of the three FYN isoforms, and more systematically by using a Principal Component Analysis (PCA) on C1-fused SRC kinases and either SKAP1 or SKAP2 mutant data.

## SKAP1 modular architecture interacting with SRC kinases

We first studied the interaction of nine SKAP1 mutants with nine SRC kinases, which has been able to previously define the SKAP2 architecture [27]. Fig 2 shows one experiment for each SKAP1 mutant. The interaction of SRC kinases systematically decreases for ΔPHLink2 (deletion from aa 107 to 326), DomDIM1 (sequence from aa 1 to 99), and Link2SH3 (sequence from aa 213 to 358) SKAP1 mutants, independently of the localization of the hemi-gaussia fusion. This strongly suggests that the domain DIM and the interdomain Link2 are necessary and sufficient for the interaction with the SRC kinases. The Link2SH3 and the DomDIM mutants have a similar loss of signal indicating that neither the interdomain Link2 nor the DIM domain are sufficient. The ΔPH (deletion from aa 107 to 210), the ΔLink1PH (deletion from aa 67 to 210), and the ΔSH3 (sequence from aa 1 to 291) mutants do not exhibit such decrease of the PPI signal supporting that they are not necessary. The W333K SKAP1 mutant, invalidated for the binding to the proline-rich binding motif, behave similarly as the ΔSH3 mutant, showing a slight increase in their binding capacity. The D120K and DIM* mutations have a slight effect on the PPI respectively increasing the signal for C1-fused SRC kinases and decreasing it for N1-fused ones.

<footnote>* Mean of normalized NLR ± SEM p = 0.05 (*); p = 0.01 (**); p = 0.001 (***); bold: p = 0.05 after Bonferroni correction.</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Study of the interaction between SKAP1 mutants and N1-fused or C1-fused SRC kinases. Graphs show wild type normalized NLR of nine SKAP1 mutants with N1-fused (blue) and C1-fused (red) SRC kinases in a logarithmic scale. Wild type normalized NLR is the ratio of the NLR between a SKAP1 mutant and one SRC kinase to the corresponding NLR between wild type SKAP1 and the same SRC kinase fused with hemi-gaussia at the same position. For each graph, the normalized NLR mean ± SEM among SRC kinases is shown as well as its level of significance (p = 0.05 (*), p = 0.01 (**), p = 0.001 (***)).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g002

In conclusion, our data suggest a SKAP1 binding architecture independent of its SH3 domain and of the status of the SRC kinases (i.e the localization of the hemi-luciferase), which clearly differs from the one we identified for SKAP2 using the same strategy [27]. S3 Table in S2 File summarizes our knowledge on the role of the different domains and interdomains of SKAP1 and SKAP2 in the interaction with SRC kinases.

## Statistical validation of the model of SKAP1 architecture

Previous hypotheses were confirmed using univariate analysis combining three experiments. For each SKAP1 mutant, three mean-value effects were tested on N1-fused SRC kinases and on C1-fused ones comparing to theoretical mean of no effect equal to one and on the difference between the signal of the C1-fused SRC kinases and that of N1-fused ones. Results were considered statistically significant only after Bonferroni's correction (Table 1).

Comparison of the deletion mutants shows that (1) only three mutants, ΔPHLink2, Dom-DIM, and Link2SH3, have a significant decrease of the PPI signal; (2) other mutants except DIM one have at least one significant increase of the PPI signal; (3) two mutants, ΔLink1PH and ΔPH have a significant increase of the N1-fused signal with a significant interaction; (4) two mutants, W333K and ΔSH3, have a significant increase of the C1-fused signal; (5) D120K mutant has a significant increase of N1-fused and C1-fused signal significantly greater for C1-fused one. These results validate the model of SKAP1 architecture in which the domain DIM and the interdomain Link2 are necessary and sufficient for binding to SRC kinases independently of the location of hemi-gaussia fusion. In addition, the similar interaction pattern of ΔLink1PH and ΔPH mutants suggests that the PH domain decreases probably indirectly the interaction with N1-fused SRC kinases. The increase of binding to SRC kinases for D120K mutant supports either a role of the intramolecular interaction between the DIM and the PH
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Study of the interaction between SKAP1 mutants and N1-fused or C1-fused SRC kinases. Graphs show wild type normalized NLR of nine SKAP1 mutants with N1-fused (blue) and C1-fused (red) SRC kinases in a logarithmic scale. Wild type normalized NLR is the ratio of the NLR between a SKAP1 mutant and one SRC kinase to the corresponding NLR between wild type SKAP1 and the same SRC kinase fused with hemi-gaussia at the same position. For each graph, the normalized NLR mean ± SEM among SRC kinases is shown as well as its level of significance (p = 0.05 (*), p = 0.01 (**), p = 0.001 (***)).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g002

In conclusion, our data suggest a SKAP1 binding architecture independent of its SH3 domain and of the status of the SRC kinases (i.e the localization of the hemi-luciferase), which clearly differs from the one we identified for SKAP2 using the same strategy [27]. S3 Table in S2 File summarizes our knowledge on the role of the different domains and interdomains of SKAP1 and SKAP2 in the interaction with SRC kinases.

## Statistical validation of the model of SKAP1 architecture

Previous hypotheses were confirmed using univariate analysis combining three experiments. For each SKAP1 mutant, three mean-value effects were tested on N1-fused SRC kinases and on C1-fused ones comparing to theoretical mean of no effect equal to one and on the difference between the signal of the C1-fused SRC kinases and that of N1-fused ones. Results were considered statistically significant only after Bonferroni's correction (Table 1).

Comparison of the deletion mutants shows that (1) only three mutants, ΔPHLink2, Dom-DIM, and Link2SH3, have a significant decrease of the PPI signal; (2) other mutants except DIM one have at least one significant increase of the PPI signal; (3) two mutants, ΔLink1PH and ΔPH have a significant increase of the N1-fused signal with a significant interaction; (4) two mutants, W333K and ΔSH3, have a significant increase of the C1-fused signal; (5) D120K mutant has a significant increase of N1-fused and C1-fused signal significantly greater for C1-fused one. These results validate the model of SKAP1 architecture in which the domain DIM and the interdomain Link2 are necessary and sufficient for binding to SRC kinases independently of the location of hemi-gaussia fusion. In addition, the similar interaction pattern of ΔLink1PH and ΔPH mutants suggests that the PH domain decreases probably indirectly the interaction with N1-fused SRC kinases. The increase of binding to SRC kinases for D120K mutant supports either a role of the intramolecular interaction between the DIM and the PH
domains [12] or the appearance of a new binding site for lipids by increasing the interaction of
the PH domain to membranes [12,13].

A statistical modeling was performed taking account of both literature curation and our
results. Six variables (see Table 2) were generated, describing function possibly affected by
mutations and a status for each variable was assigned for each mutant. However due to overfit-
ting, these six variables cannot be analyzed together, and we decided to suppress the variable
DIM⁻. Note that the loss of PH domain combined ΔLink1PH and ΔPH mutant data. A two-
way ANOVA analysis was carried out after lognormal transformation of the dependent vari-
able on these five variables with either N1-fused or C1-fused SRC kinase data. This transfor-
mation was chosen since most data are linearized after log-log transformation in differential
interaction scatterplot (see S1 Fig in S1 File). An explicative model was generated from a full
model by a step-by-step procedure that removed at each step the lowest nonsignificant variable
until no nonsignificant variable was present (Table 2).

For C1-fused SRC kinase data, four variables were significant, the SKAP1 module consist-
ing of the DIM and Link2 regions necessary and sufficient to bind to SRC kinases, D120K,
ΔSH3, and W333K. For N1-fused SRC kinase data, the five variables, the SKAP1 module,
D120K, ΔSH3, W333K, and the loss of PH domain were significant. This ANOVA thus cor-
roborates the SKAP1 interaction module as well as the effects of the D120K, W333K, and
ΔSH3 mutants as defined in the univariate analysis and supports that the SRC kinase status
affects only slightly their binding capacity to SKAP1. The differences between the results
obtained from the study of N1-fused and C1-fused SRC kinases support the investigation of
both localizations of the fusion protein. They also support the idea that both SRC kinase locali-
zation and activation affect these PPIs.

Comparing SKAP1 and SKAP2 architecture

To analyze the difference between the two architectures more finely, we compared the effect of
each mutant and wild type form between SKAP1 and SKAP2. In practice, the NLR generated
by the binding of a SKAP1 mutant or wild type SKAP1 to an SRC kinase has been divided by
the NLR of same SKAP2 mutant or wild type SKAP2 to the same kinase. This ratio was named
SKAP2 normalized NLR. Fig 3 shows the result of one experiment for each of the ten pairwise
proteins. For each of the ten data, the mean-value of NLR for the nine SRC kinases was com-
pared to the theorical value of no effect (equal to 1) depending on the localization of the hemi-
gaussia fusion. First, we compared the pattern of significant differences between each mutant
and the wildtype protein to classify the data. As previously described [17], SKAP1 has a signifi-
cant decrease of the PPI signal for C1-fused SRC kinases compared to SKAP2. Five mutants,
D120K, DIM⁻, ΔPH, ΔLink2PH, and W333K, have a similar pattern suggesting that these dele-
tions or mutations are not responsible for the difference between SKAP1 and SKAP2. Note for
the DIM⁻ mutant that its mean of NLR for C1-fused SRC kinase is greater than that for wild-
type SKAP1 (mC1 = 0.75 versus 0.39). For the ΔSH3 and ΔLink1PH mutants, SKAP1 PPI sig-
nals to N1-fused SRC kinases are stronger than the SKAP2 ones. Link2SH3 and DomDIM
mutants have a slightly different pattern with no significant difference of signal for N1-fused
and C1-fused SRC kinase. From the differential interaction scatterplot between SKAP1 versus
SKAP2, the position of each barycenter and the corresponding regression line were calculated
(S2 Fig in S1 File). The geometric mean of the NLRs from the SKAP1 and SKAP2 data defines
the coordinates of the barycenter along the y-axis and x-axis, respectively (S3 Fig in S1 File). If
the position of this barycenter varies on a line parallel to the no-effect line with a slope equal to
1, the difference between the NLR mean-values of SKAP1 and SKAP2 will remain constant.
All other changes in barycenter position will affect this difference and support that the
domains [12] or the appearance of a new binding site for lipids by increasing the interaction of the PH domain to membranes [12,13].

A statistical modeling was performed taking account of both literature curation and our results. Six variables (see Table 2) were generated, describing function possibly affected by mutations and a status for each variable was assigned for each mutant. However due to overfitting, these six variables cannot be analyzed together, and we decided to suppress the variable DIM⁻. Note that the loss of PH domain combined ΔLink1PH and ΔPH mutant data. A two-way ANOVA analysis was carried out after lognormal transformation of the dependent variable on these five variables with either N1-fused or C1-fused SRC kinase data. This transformation was chosen since most data are linearized after log-log transformation in differential interaction scatterplot (see S1 Fig in S1 File). An explicative model was generated from a full model by a step-by-step procedure that removed at each step the lowest nonsignificant variable until no nonsignificant variable was present (Table 2).

For C1-fused SRC kinase data, four variables were significant, the SKAP1 module consisting of the DIM and Link2 regions necessary and sufficient to bind to SRC kinases, D120K, ΔSH3, and W333K. For N1-fused SRC kinase data, the five variables, the SKAP1 module, D120K, ΔSH3, W333K, and the loss of PH domain were significant. This ANOVA thus corroborates the SKAP1 interaction module as well as the effects of the D120K, W333K, and ΔSH3 mutants as defined in the univariate analysis and supports that the SRC kinase status affects only slightly their binding capacity to SKAP1. The differences between the results obtained from the study of N1-fused and C1-fused SRC kinases support the investigation of both localizations of the fusion protein. They also support the idea that both SRC kinase localization and activation affect these PPIs.

## Comparing SKAP1 and SKAP2 architecture

To analyze the difference between the two architectures more finely, we compared the effect of each mutant and wild type form between SKAP1 and SKAP2. In practice, the NLR generated by the binding of a SKAP1 mutant or wild type SKAP1 to an SRC kinase has been divided by the NLR of same SKAP2 mutant or wild type SKAP2 to the same kinase. This ratio was named SKAP2 normalized NLR. Fig 3 shows the result of one experiment for each of the ten pairwise proteins. For each of the ten data, the mean-value of NLR for the nine SRC kinases was compared to the theoretical value of no effect (equal to 1) depending on the localization of the hemigaussia fusion. First, we compared the pattern of significant differences between each mutant and the wildtype protein to classify the data. As previously described [17], SKAP1 has a significant decrease of the PPI signal for C1-fused SRC kinases compared to SKAP2. Five mutants, D120K, DIM⁻, ΔPH, ΔLink2PH, and W333K, have a similar pattern suggesting that these deletions or mutations are not responsible for the difference between SKAP1 and SKAP2. Note for the DIM⁻ mutant that its mean of NLR for C1-fused SRC kinase is greater than that for wild-type SKAP1 ($m_{C1} = 0.75$ versus 0.39). For the ΔSH3 and ΔLink1PH mutants, SKAP1 PPI signals to N1-fused SRC kinases are stronger than the SKAP2 ones. Link2SH3 and DomDIM mutants have a slightly different pattern with no significant difference of signal for N1-fused and C1-fused SRC kinase. From the differential interaction scatterplot between SKAP1 versus SKAP2, the position of each barycenter and the corresponding regression line were calculated (S2 Fig in S1 File). The geometric mean of the NLRs from the SKAP1 and SKAP2 data defines the coordinates of the barycenter along the y-axis and x-axis, respectively (S3 Fig in S1 File). If the position of this barycenter varies on a line parallel to the no-effect line with a slope equal to 1, the difference between the NLR mean-values of SKAP1 and SKAP2 will remain constant. All other changes in barycenter position will affect this difference and support that the
<table><caption>Table 2. Statistical modeling of the interaction between SKAP1 and SRC kinases.</caption><thead><tr><th rowspan="2">Variable</th><th colspan="5">N1-Loc</th><th colspan="5">C1-Loc</th></tr><tr><th>Partial SS*</th><th>df*</th><th>MS*</th><th>F (sense)*</th><th>P</th><th>Partial SS</th><th>df</th><th>MS</th><th>F (sense)</th><th>P</th></tr></thead><tbody><tr><td>Explanatory model</td><td>146.81</td><td>5</td><td>29.36</td><td>130.05</td><td>7.20 10<sup>-66</sup></td><td>142.21</td><td>4</td><td>35.55</td><td>136.10</td><td>2.62 10<sup>-60</sup></td></tr><tr><td>Loss of PH domain</td><td>5.26</td><td>1</td><td>5.26</td><td>23.75 (+)</td><td>2.01 10<sup>-6</sup></td><td></td><td></td><td></td><td></td><td>NS</td></tr><tr><td>D120K</td><td>2.08</td><td>1</td><td>2.08</td><td>9.21 (+)</td><td>0.0027</td><td>4.44</td><td>1</td><td>4.44</td><td>16.99 (+)</td><td>0.0001</td></tr><tr><td>ΔSH3</td><td>1.41</td><td>1</td><td>1.41</td><td>6.23 (+)</td><td>0.0132</td><td>33.19</td><td>1</td><td>33.19</td><td>127.07(+)</td><td>6.74 10<sup>-24</sup></td></tr><tr><td>W333K</td><td>1.61</td><td>1</td><td>1.61</td><td>7.14 (+)</td><td>0.0081</td><td>1.35</td><td>1</td><td>1.35</td><td>5.17 (+)</td><td>0.0239</td></tr><tr><td>SKAP1 module</td><td>28.46</td><td>1</td><td>28.46</td><td>126.03</td><td>9.89 10<sup>-24</sup></td><td>37.02</td><td>1</td><td>37.02</td><td>141.72</td><td>6.04 10<sup>-26</sup></td></tr><tr><td>Residual</td><td>53.51</td><td>237</td><td>0.24</td><td></td><td></td><td>62.17</td><td>238</td><td>0.26</td><td></td><td></td></tr></tbody></table>

* Partial SS: Partial Sum of Squares; df: Degree of fredoom; MS: Mean of Squares.

* F (sense): Fisher statistic (+: increase).

* NS non-significant (deleted variable during the procedure of minimization).

https://doi.org/10.1371/journal.pone.0296230.t002

corresponding mutant has different interactions with SRC kinases than its wild type protein. The position of both barycenters are similar to that of wildtype protein for D120K, ΔPH and W333K mutants. This position of both baricenters moves mainly according to the line of no effect for ΔPHLink2 mutant and probably also for DomDIM one (see C1-fused data and below) suggesting no difference between SKAP1 and SKAP2 data. The position of barycenters for both N1-fused and C1-fused data moves mainly according to the x-axis supporting an effect of ΔLink1PH and ΔSH3 mutants. Only for C1-fused data, a similar x-axis movement is detected for DomDIM and Link2SH3 mutants even if we could not exclude a similar effect for DomDIM mutant. In conclusion of these two analyses, ΔLink1PH, ΔSH3, and Link2SH3

<figure><img src="image_7.png" /><figcaption>**Fig 3.** Comparison of the interaction with SRC kinases between SKAP1 and SKAP2 mutants. Graphs show SKAP2 normalized NLR of wild type SKAP1 and nine SKAP1 mutants with N1-fused (blue) and C1-fused (red) SRC kinases in a logarithmic scale. SKAP2 normalized NLR is the ratio of the NLR between a SKAP1 mutant and one SRC kinase to the corresponding NLR between the same SKAP2 mutant and the same SRC kinase fused with hemi-gaussia at the same position. For each graph, the normalized NRL mean ± SEM among SRC kinases is shown as its level of significance (*p* = 0.05 (*), *p* = 0.01 (**), *p* = 0.001 (**)).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g003
<table><caption>Table 2. Statistical modeling of the interaction between SKAP1 and SRC kinases.</caption><thead><tr><th rowspan="2">Variable</th><th colspan="3">N1-Loc</th><th colspan="3">C1-Loc</th></tr><tr><th>Partial SS*</th><th>df*</th><th>MS*</th><th>F (sense)*</th><th>P</th><th>Partial SS</th><th>df</th><th>MS</th><th>F (sense)</th><th>P</th></tr></thead><tbody><tr><td>Explanatory model</td><td>146.81</td><td>5</td><td>29.36</td><td>130.05</td><td>7.20 10<sup>-66</sup></td><td>142.21</td><td>4</td><td>35.55</td><td>136.10</td><td>2.62 10<sup>-60</sup></td></tr><tr><td>Loss of PH domain</td><td>5.26</td><td>1</td><td>5.26</td><td>23.75 (+)</td><td>2.01 10<sup>-6</sup></td><td></td><td></td><td></td><td></td><td>NS</td></tr><tr><td>D120K</td><td>2.08</td><td>1</td><td>2.08</td><td>9.21 (+)</td><td>0.0027</td><td>4.44</td><td>1</td><td>4.44</td><td>16.99 (+)</td><td>0.0001</td></tr><tr><td>ΔSH3</td><td>1.41</td><td>1</td><td>1.41</td><td>6.23 (+)</td><td>0.0132</td><td>33.19</td><td>1</td><td>33.19</td><td>127.07(+)</td><td>6.74 10<sup>-24</sup></td></tr><tr><td>W333K</td><td>1.61</td><td>1</td><td>1.61</td><td>7.14 (+)</td><td>0.0081</td><td>1.35</td><td>1</td><td>1.35</td><td>5.17 (+)</td><td>0.0239</td></tr><tr><td>SKAP1 module</td><td>28.46</td><td>1</td><td>28.46</td><td>126.03</td><td>9.89 10<sup>-24</sup></td><td>37.02</td><td>1</td><td>37.02</td><td>141.72</td><td>6.04 10<sup>-26</sup></td></tr><tr><td>Residual</td><td>53.51</td><td>237</td><td>0.24</td><td></td><td></td><td>62.17</td><td>238</td><td>0.26</td><td></td><td></td></tr></tbody></table>

* Partial SS: Partial Sum of Squares; df: Degree of fredoom; MS: Mean of Squares.

* F (sense): Fisher statistic (+: increase).

* NS non-significant (deleted variable during the procedure of minimization).

https://doi.org/10.1371/journal.pone.0296230.t002

corresponding mutant has different interactions with SRC kinases than its wild type protein. The position of both barycenters are similar to that of wildtype protein for D120K, ΔPH and W333K mutants. This position of both baricenters moves mainly according to the line of no effect for ΔPHLink2 mutant and probably also for DomDIM one (see C1-fused data and below) suggesting no difference between SKAP1 and SKAP2 data. The position of barycenters for both N1-fused and C1-fused data moves mainly according to the x-axis supporting an effect of ΔLink1PH and ΔSH3 mutants. Only for C1-fused data, a similar x-axis movement is detected for DomDIM and Link2SH3 mutants even if we could not exclude a similar effect for DomDIM mutant. In conclusion of these two analyses, ΔLink1PH, ΔSH3, and Link2SH3

<figure><img src="image_6.png" /><figcaption>**Fig 3.** Comparison of the interaction with SRC kinases between SKAP1 and SKAP2 mutants. Graphs show SKAP2 normalized NLR of wild type SKAP1 and nine SKAP1 mutants with N1-fused (blue) and C1-fused (red) SRC kinases in a logarithmic scale. SKAP2 normalized NLR is the ratio of the NLR between a SKAP1 mutant and one SRC kinase to the corresponding NLR between the same SKAP2 mutant and the same SRC kinase fused with hemi-gaussia at the same position. For each graph, the normalized NRL mean ± SEM among SRC kinases is shown as its level of significance (*p* = 0.05 (*), *p* = 0.01 (**), *p* = 0.001 (**)).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g003
mutants show clear differences between SKAP1 and SKAP2 data. Conflictual results have been obtained for DomDIM and DIM⁻ mutants.

A statistical validation was performed by comparing the mean-value of SKAP2 normalized NLR for each mutant with that of wildtype protein for N1-fused or C1-fused SRC kinases. Table 3 shows the results. It confirms that ΔLink1PH and ΔSH3 mutants have a significantly greater mean-value of SKAP2 normalized NLR than that for wildtype protein independently of the localization of the hemi-gaussa fusion. Link2SH3 and DIM⁻ mutants have a significantly greater mean-value of SKAP2 normalized NLR than that for wildtype protein only for C1 data. These results showed that these four mutations have different roles between SKAP1 and SKAP2 since SKAP1 mutant and wildtype SKAP1 NLRs are normalized according the corresponding NLR of SKAP2 mutant and wildtype SKAP2.

In conclusion, these results strongly suggest that the SH3 domain and the interdomain Link1 have different functions between SKAP1 and SKAP2 proteins. These results support also that dimerization affects differently the interaction to C1-fused SRC kinases between SKAP1 and SKAP2 proteins. Furthermore, the DIM domain has also a major role to explain the differences detected between SKAP1 and SKAP2 proteins for their binding capacity to SRC kinases depending on the localization of the hemi-gaussa fusion. This result explains why the Link2SH3 mutant, which lost the DIM domain, has a significantly greater PPI signal for C1-fused SRC kinases.

## Role of a SKAP2-specific motif located in the interdomain Link1

Our previous results strongly suggest that the interdomain Link1 affects specifically the interaction with SRC kinases for SKAP2. Search for a motif specific to SKAP2 and conserved in its phylogeny identified a motif that has been defined as a site of phosphorylation for Casein kinase 1 [70]. This site has been localized in a predicted α-helix according to I-Tasser (Fig 4A), which is absent in SKAP1 such specific to SKAP2. The two serine residues of this site were mutated to either proline or alanine to destroy or not the predicted α-helix. Fig 4B shows that both mutants increased the PPI signal of SRC kinases independently of the localization of hemi-gaussa fusion. This result supports the view that the phosphorylation of this serine affects the interaction with SRC kinases and not the putative secondary structure. This motif is probably at least one of the elements, which explains the specific role of the interdomain Link1 of SKAP2.

<table><caption>Table 3. Statistical analysis comparing effect of each mutation between SKAP1 and SKAP2.</caption><thead><tr><th>Mutant</th><th>Nb Exp (n)</th><th>ΔN1-Mutant-N1-SKAP1*</th><th>ΔC1-Mutant-C1-SKAP1*</th></tr></thead><tbody><tr><td>D120K</td><td>2 (18)</td><td>-0.05 ± 0.17</td><td>-0.01 ± 0.07</td></tr><tr><td>DIM-</td><td>2 (18)</td><td>0.38 ± 0.18</td><td><strong>0.36 ± 0.06</strong></td></tr><tr><td>ΔLink1PH</td><td>2 (18)</td><td><strong>2.00 ± 0.34</strong></td><td><strong>0.57 ± 0.10</strong></td></tr><tr><td>ΔPH</td><td>1 (9)</td><td>-0.24 ± 0.15</td><td>-0.10 ± 0.06</td></tr><tr><td>ΔPHLink2</td><td>2 (18)</td><td>-0.22 ± 0.14</td><td>-0.04 ± 0.06</td></tr><tr><td>ΔSH3</td><td>2 (18)</td><td><strong>7.39 ± 1.48</strong></td><td><strong>1.07 ± 0.22</strong></td></tr><tr><td>DomDIM</td><td>2 (18)</td><td>0.37 ± 0.18</td><td>0.23 ± 0.16</td></tr><tr><td>Link2SH3</td><td>2 (18)</td><td>0.03 ± 0.14</td><td><strong>0.60 ± 0.10</strong></td></tr><tr><td>W333K</td><td>2 (18)</td><td>0.24 ± 0.17</td><td>0.16 ± 0.07</td></tr></tbody></table>

<footnote>* Difference of means of normalized NLR ± SEM. Bold: p = 0.05 after Bonferroni correction.</footnote>
mutants show clear differences between SKAP1 and SKAP2 data. Conflictual results have been obtained for DomDIM and DIM⁻ mutants.

A statistical validation was performed by comparing the mean-value of SKAP2 normalized NLR for each mutant with that of wildtype protein for N1-fused or C1-fused SRC kinases. Table 3 shows the results. It confirms that ΔLink1PH and ΔSH3 mutants have a significantly greater mean-value of SKAP2 normalized NLR than that for wildtype protein independently of the localization of the hemi-gaussa fusion. Link2SH3 and DIM⁻ mutants have a significantly greater mean-value of SKAP2 normalized NLR than that for wildtype protein only for C1 data. These results showed that these four mutations have different roles between SKAP1 and SKAP2 since SKAP1 mutant and wildtype SKAP1 NLRs are normalized according the corresponding NLR of SKAP2 mutant and wildtype SKAP2.

In conclusion, these results strongly suggest that the SH3 domain and the interdomain Link1 have different functions between SKAP1 and SKAP2 proteins. These results support also that dimerization affects differently the interaction to C1-fused SRC kinases between SKAP1 and SKAP2 proteins. Furthermore, the DIM domain has also a major role to explain the differences detected between SKAP1 and SKAP2 proteins for their binding capacity to SRC kinases depending on the localization of the hemi-gaussa fusion. This result explains why the Link2SH3 mutant, which lost the DIM domain, has a significantly greater PPI signal for C1-fused SRC kinases.

## Role of a SKAP2-specific motif located in the interdomain Link1

Our previous results strongly suggest that the interdomain Link1 affects specifically the interaction with SRC kinases for SKAP2. Search for a motif specific to SKAP2 and conserved in its phylogeny identified a motif that has been defined as a site of phosphorylation for Casein kinase 1 [70]. This site has been localized in a predicted α-helix according to I-Tasser (Fig 4A), which is absent in SKAP1 such specific to SKAP2. The two serine residues of this site were mutated to either proline or alanine to destroy or not the predicted α-helix. Fig 4B shows that both mutants increased the PPI signal of SRC kinases independently of the localization of hemi-gaussa fusion. This result supports the view that the phosphorylation of this serine affects the interaction with SRC kinases and not the putative secondary structure. This motif is probably at least one of the elements, which explains the specific role of the interdomain Link1 of SKAP2.

<table><caption>Table 3. Statistical analysis comparing effect of each mutation between SKAP1 and SKAP2.</caption><thead><tr><th>Mutant</th><th>Nb Exp (n)</th><th>ΔN1-Mutant-N1-SKAP1*</th><th>ΔC1-Mutant-C1-SKAP1*</th></tr></thead><tbody><tr><td>D120K</td><td>2 (18)</td><td>-0.05 ± 0.17</td><td>-0.01 ± 0.07</td></tr><tr><td>DIM-</td><td>2 (18)</td><td>0.38 ± 0.18</td><td><strong>0.36 ± 0.06</strong></td></tr><tr><td>ΔLink1PH</td><td>2 (18)</td><td><strong>2.00 ± 0.34</strong></td><td><strong>0.57 ± 0.10</strong></td></tr><tr><td>ΔPH</td><td>1 (9)</td><td>-0.24 ± 0.15</td><td>-0.10 ± 0.06</td></tr><tr><td>ΔPHLink2</td><td>2 (18)</td><td>-0.22 ± 0.14</td><td>-0.04 ± 0.06</td></tr><tr><td>ΔSH3</td><td>2 (18)</td><td><strong>7.39 ± 1.48</strong></td><td><strong>1.07 ± 0.22</strong></td></tr><tr><td>DomDIM</td><td>2 (18)</td><td>0.37 ± 0.18</td><td>0.23 ± 0.16</td></tr><tr><td>Link2SH3</td><td>2 (18)</td><td>0.03 ± 0.14</td><td><strong>0.60 ± 0.10</strong></td></tr><tr><td>W333K</td><td>2 (18)</td><td>0.24 ± 0.17</td><td>0.16 ± 0.07</td></tr></tbody></table>

<footnote>* Difference of means of normalized NLR ± SEM. Bold: p = 0.05 after Bonferroni correction.</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig 4. Role of a conserved site of phosphorylation for Casein kinase 1 in Link1 interdomain.** (A) The N-terminal part of one chain of SKAP2 from the DIM domain to the PH domain has been modelized by I-Tasser server [66]. Conservation has been evaluated by Consurf website [67,68] using a scale from highly variable residues (dark blue) to highly conserved one (dark red) passing through white. The phosphorylation site with its predicted α-helix and the two highly conserved serine residues (arrows in the lower left) is in the Link1 interdomain between the DIM domain with its predicted α-helix (lower right) and the PH domain (upper). (B) Graphs show wild type normalized NLR of two SKAP2 mutants, S87AS90A and S87PS90P, with N1-fused (blue) and C1-fused (red) SRC kinases in a logarithmic scale. Wild type normalized NLR is the ratio of the NLR between a SKAP2 mutant and one SRC kinase to the corresponding NLR between wild type SKAP2 and the same SRC kinase fused with hemi-gaussia at the same position. For each graph, the normalized NRL mean ± SEM among SRC kinases is shown as its level of significance (p = 0.05 (*), p = 0.01 (**), p = 0.001 (***)).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g004

## Role of two conserved motifs located in the interdomain Link2

The interdomain Link2 plays a major role for the binding of SRC kinases on both SKAP1 and SKAP2 proteins supporting that similar motifs exist for the two proteins. Search for motifs conserved in both SKAP1 or SKAP2 phylogeny found two motifs consisting of a binding motif for a SH2 domain followed or forward by a glutamine acid-rich region [70,71]. We decided to test their effect on the binding to SRC kinases by delating one of them or both in the two proteins. Fig 5 shows the result of four experiments for SKAP1 (ΔLink2A: deletion from aa 262 to 274; ΔLink2B: deletion from aa 219 to 235) and six experiments for SKAP2 (ΔLink2A: deletion from aa 261 to 272; ΔLink2B: deletion from aa 228 to 240) mutants. For SKAP2, we have also tested the mutation of three tyrosine to phenylalanine located at position 75, 237 and 261, which have been predicted to be inside a binding motif for a SH2 domain. This mutant, named 3YF, was mutated at the three predicted binding motifs for a SH2 domain. For each of the two proteins, deletion of one motif slightly decreases the binding to SRC kinases and deletion of both motifs has a greater effect than the deletion of each separatly. Note also that the binding capacity to SRC kinases appears similar or slightly greater for the 3YF SKAP2 mutant than for ΔLink2AB SKAP2 mutant. S3 Table in S2 File, which shows the results of the statistical analysis, supports this analysis. Interestingly, the deletion of motif A affects differently the level of the PPI signal than that of motif B depending on SRC kinases. This difference frequently persists in both SKAP1 and SKAP2 proteins. For example, the decrease of binding for HCK, LYN, FYN, and SRMS kinases is greater after deletion of motif A than that of motif B in
<figure><img src="image_7.png" /><figcaption>**Fig 4. Role of a conserved site of phosphorylation for Casein kinase 1 in Link1 interdomain.** (A) The N-terminal part of one chain of SKAP2 from the DIM domain to the PH domain has been modelized by I-Tasser server [66]. Conservation has been evaluated by Consurf website [67,68] using a scale from highly variable residues (dark blue) to highly conserved one (dark red) passing through white. The phosphorylation site with its predicted α-helix and the two highly conserved serine residues (arrows in the lower left) is in the Link1 interdomain between the DIM domain with its predicted α-helix (lower right) and the PH domain (upper). (B) Graphs show wild type normalized NLR of two SKAP2 mutants, S87AS90A and S87PS90P, with N1-fused (blue) and C1-fused (red) SRC kinases in a logarithmic scale. Wild type normalized NLR is the ratio of the NLR between a SKAP2 mutant and one SRC kinase to the corresponding NLR between wild type SKAP2 and the same SRC kinase fused with hemi-gaussia at the same position. For each graph, the normalized NRL mean ± SEM among SRC kinases is shown as its level of significance (p = 0.05 (*), p = 0.01 (**), p = 0.001 (***)).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g004

## Role of two conserved motifs located in the interdomain Link2

The interdomain Link2 plays a major role for the binding of SRC kinases on both SKAP1 and SKAP2 proteins supporting that similar motifs exist for the two proteins. Search for motifs conserved in both SKAP1 or SKAP2 phylogeny found two motifs consisting of a binding motif for a SH2 domain followed or forward by a glutamine acid-rich region [70,71]. We decided to test their effect on the binding to SRC kinases by delating one of them or both in the two proteins. Fig 5 shows the result of four experiments for SKAP1 (ΔLink2A: deletion from aa 262 to 274; ΔLink2B: deletion from aa 219 to 235) and six experiments for SKAP2 (ΔLink2A: deletion from aa 261 to 272; ΔLink2B: deletion from aa 228 to 240) mutants. For SKAP2, we have also tested the mutation of three tyrosine to phenylalanine located at position 75, 237 and 261, which have been predicted to be inside a binding motif for a SH2 domain. This mutant, named 3YF, was mutated at the three predicted binding motifs for a SH2 domain. For each of the two proteins, deletion of one motif slightly decreases the binding to SRC kinases and deletion of both motifs has a greater effect than the deletion of each separatly. Note also that the binding capacity to SRC kinases appears similar or slightly greater for the 3YF SKAP2 mutant than for ΔLink2AB SKAP2 mutant. S3 Table in S2 File, which shows the results of the statistical analysis, supports this analysis. Interestingly, the deletion of motif A affects differently the level of the PPI signal than that of motif B depending on SRC kinases. This difference frequently persists in both SKAP1 and SKAP2 proteins. For example, the decrease of binding for HCK, LYN, FYN, and SRMS kinases is greater after deletion of motif A than that of motif B in
<figure><img src="image_9.png" /><figcaption>**Fig 5.** Study of the deletion of the two conserved motif in the interdomain Link2 for SKAP1 and SKAP2 protein.

Graphs show the mean of wild type normalized NLR of three SKAP1 deletion mutants (4 experiments) and four SKAP2 mutants (7 experiments) with N1-fused (blue) and C1-fused (red) SRC kinases in a logarithmic scale. Wild type normalized NLR is the ratio of the NLR between a mutant and one SRC kinase to the corresponding NLR between wild type protein and the same SRC kinase fused with hemi-gaussia at the same position. ΔLink2A: Deletion of motif A; ΔLink2B: Deletion of motif B; ΔLink2AB: Deletion of both motifs; 3YF: Mutation to phenylalanine of the three tyrosine residues located at position 75, 237, and 261 of SKAP2.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g005

both proteins and the opposite is true for LCK. But the decrease of binding for BLK and FRK kinases is greater after deletion of motif A than that of motif B only for SKAP2.

## Dissociating the binding capacity of the SH2 domain from activation for HCK and LCK kinase

These two conserved motifs contain a binding site for SH2 domain which role for binding to SRC kinases was evaluated. It is not a trivial task because the SH2 domain of SRC kinases also plays a major role for their activation (Fig 1C). To dissociate the activation to the binding capacity of the SH2 domain of the kinase, we decided to study two mutants, either Y501F and R150A for HCK kinase or Y505F and R154A for LCK one. The R150A HCK/ R154A LCK mutant has lost the binding capacity of the SH2 domain and only activation of HCK/LCK kinase occurs (Fig 1C). The Y501F HCK/ Y505F LCK mutant has lost the binding site of its SH2 domain and activation of the HCK/LCK kinase occurs associated to the capacity of its SH2 domain to interact with other proteins. To dissociate a direct binding effect of SKAP1/ SKAP2 on SRC kinases from an indirect effect of other proteins, whose binding modified the interaction between SKAP1/SKAP2 and SRC kinases, we also study their ΔLink2AB mutant. Furthermore, the study of this mutant is necessary to exclude than another motif than the SH2/ ΔLink2AB module affects the signal. The indirect role is detected by a difference of ΔLink2AB binding between the mutant bearing the binding motif to SH2 domain (R150A HCK or R154A LCK mutant) and that of the SH2 binding site (Y501F HCK or Y505F LCK mutant) and the direct role by a decrease of ΔLink2AB binding compare to wild type adaptor with mutant bearing the binding motif to SH2 domain that disappears with mutant of the SH2 binding site. In the remainder of this paragraph, we will use the term activation in reference to these 2 mutations without having validated it experimentally. Fig 6 shows the result of one
<figure><img src="image_8.png" /><figcaption>**Fig 5.** Study of the deletion of the two conserved motif in the interdomain Link2 for SKAP1 and SKAP2 protein.
Graphs show the mean of wild type normalized NLR of three SKAP1 deletion mutants (4 experiments) and four SKAP2 mutants (7 experiments) with N1-fused (blue) and C1-fused (red) SRC kinases in a logarithmic scale. Wild type normalized NLR is the ratio of the NLR between a mutant and one SRC kinase to the corresponding NLR between wild type protein and the same SRC kinase fused with hemi-gaussia at the same position. ΔLink2A: Deletion of motif A; ΔLink2B: Deletion of motif B; ΔLink2AB: Deletion of both motifs; 3YF: Mutation to phenylalanine of the three tyrosine residues located at position 75, 237, and 261 of SKAP2.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g005

both proteins and the opposite is true for LCK. But the decrease of binding for BLK and FRK kinases is greater after deletion of motif A than that of motif B only for SKAP2.

## Dissociating the binding capacity of the SH2 domain from activation for HCK and LCK kinase

These two conserved motifs contain a binding site for SH2 domain which role for binding to SRC kinases was evaluated. It is not a trivial task because the SH2 domain of SRC kinases also plays a major role for their activation (Fig 1C). To dissociate the activation to the binding capacity of the SH2 domain of the kinase, we decided to study two mutants, either Y501F and R150A for HCK kinase or Y505F and R154A for LCK one. The R150A HCK/ R154A LCK mutant has lost the binding capacity of the SH2 domain and only activation of HCK/LCK kinase occurs (Fig 1C). The Y501F HCK/ Y505F LCK mutant has lost the binding site of its SH2 domain and activation of the HCK/LCK kinase occurs associated to the capacity of its SH2 domain to interact with other proteins. To dissociate a direct binding effect of SKAP1/ SKAP2 on SRC kinases from an indirect effect of other proteins, whose binding modified the interaction between SKAP1/SKAP2 and SRC kinases, we also study their ΔLink2AB mutant. Furthermore, the study of this mutant is necessary to exclude than another motif than the SH2/ ΔLink2AB module affects the signal. The indirect role is detected by a difference of ΔLink2AB binding between the mutant bearing the binding motif to SH2 domain (R150A HCK or R154A LCK mutant) and that of the SH2 binding site (Y501F HCK or Y505F LCK mutant) and the direct role by a decrease of ΔLink2AB binding compare to wild type adaptor with mutant bearing the binding motif to SH2 domain that disappears with mutant of the SH2 binding site. In the remainder of this paragraph, we will use the term activation in reference to these 2 mutations without having validated it experimentally. Fig 6 shows the result of one
<figure><img src="image_10.png" /><figcaption>**Fig 6. Dissociating the role of activation and inactivation of its SH2 domain in the interaction with SKAP1 and SKAP2.** Graphs show wild type normalized NLR of two SKAP2 mutants (upper) and two SKAP1 deletion mutants (lower) with either N1-fused and C1-fused wild type HCK and HCK mutants (A) or N1-fused and C1-fused wild type LCK and LCK mutants (B) in a logarithmic scale. Wild type normalized NLR is the ratio of the NLR between a couple of SKAP2 and kinase mutants and the corresponding NLR between wild type SKAP2 and wild type kinase fused with hemi-gaussia at the same position. Similar normalization was performed for SKAP1 data. ΔLink2AB: Deletion of motif A and B in interdomain Link2; 3YF: Mutation to phenylalanine of the three tyrosine residues located at position 75, 237, and 261 of SKAP2; ΔPHLink2: Deletion of the PH domain and interdomain Link2: DIMPH: N-terminal part of SKAP2 containing the DIM and the PH domains; Y501F HCK and Y505F LCK mutations inactivate the intramolecular binding to the SH2 domain in contrast to R150A HCK and R154A LCK mutations, which inactivate the binding capacity of SH2 domain to its phosphorylated motif.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g006

experiment for SKAP2 (upper) or SKAP1 (lower) challenged against HCK (A) or LCK(B) mutants. First, we will analyze HCK data (Fig 6A). For SKAP2 (upper), the activation of the HCK kinase shown by comparing results of wildtype HCK to R150A mutant increases the interaction between SKAP2 and HCK independently of the localization of the hemi-gaussia fusion. Similar level of interaction of R150A HCK mutant with either wildtype SKAP2 or with its ΔLink2AB mutant supports the view that no other motif than the SH2/ΔLink2AB module affects the signal. The role of the SH2 domain obtained by comparing results of Y501F and R150A HCK mutants shows that its effect on the binding to SKAP2 is moderate (N1-fused data) or null (C1-fused data). Study of ΔLink2AB mutant confirms that results and suggests that the moderate effect shown in N1-fused Y501F data is due to an indirect effect. The study of 3YF SKAP2 mutant suggests also that the glutamic acid-rich region plays a role since differences of signal between ΔLink2AB mutant and 3YF one was present for wildtype HCK and N1-fused Y501F HCK mutant. However, we cannot exclude a role of Y75F mutation. Similar analysis for SKAP1 (Fig 6A lower) shows different results. If activation of HCK has a similar effect than for SKAP2, the role of the SH2 domain was completely different. Similar level of interaction of R150A HCK mutant with either wildtype SKAP1 or with its ΔLink2AB mutant supports the view that no other motif than the SH2/ΔLink2AB module affects the signal. The difference of PPI signal between Y501F HCK mutant with either wildtype SKAP1 or ΔLink2AB mutant suggests that the HCK SH2 domain has an effect. Further analysis supports both a direct binding effect on SKAP1 and an indirect one due to others interactors, which modified the interaction between SKAP1 and HCK. Note that the results for ΔPHLink2 mutant is
<figure><img src="image_9.png" /><figcaption>**Fig 6. Dissociating the role of activation and inactivation of its SH2 domain in the interaction with SKAP1 and SKAP2.** Graphs show wild type normalized NLR of two SKAP2 mutants (upper) and two SKAP1 deletion mutants (lower) with either N1-fused and C1-fused wild type HCK and HCK mutants (A) or N1-fused and C1-fused wild type LCK and LCK mutants (B) in a logarithmic scale. Wild type normalized NLR is the ratio of the NLR between a couple of SKAP2 and kinase mutants and the corresponding NLR between wild type SKAP2 and wild type kinase fused with hemi-gaussia at the same position. Similar normalization was performed for SKAP1 data. ΔLink2AB: Deletion of motif A and B in interdomain Link2; 3YF: Mutation to phenylalanine of the three tyrosine residues located at position 75, 237, and 261 of SKAP2; ΔPHLink2: Deletion of the PH domain and interdomain Link2: DIMPH: N-terminal part of SKAP2 containing the DIM and the PH domains; Y501F HCK and Y505F LCK mutations inactivate the intramolecular binding to the SH2 domain in contrast to R150A HCK and R154A LCK mutations, which inactivate the binding capacity of SH2 domain to its phosphorylated motif.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g006

experiment for SKAP2 (upper) or SKAP1 (lower) challenged against HCK (A) or LCK(B) mutants. First, we will analyze HCK data (Fig 6A). For SKAP2 (upper), the activation of the HCK kinase shown by comparing results of wildtype HCK to R150A mutant increases the interaction between SKAP2 and HCK independently of the localization of the hemi-gaussia fusion. Similar level of interaction of R150A HCK mutant with either wildtype SKAP2 or with its ΔLink2AB mutant supports the view that no other motif than the SH2/ΔLink2AB module affects the signal. The role of the SH2 domain obtained by comparing results of Y501F and R150A HCK mutants shows that its effect on the binding to SKAP2 is moderate (N1-fused data) or null (C1-fused data). Study of ΔLink2AB mutant confirms that results and suggests that the moderate effect shown in N1-fused Y501F data is due to an indirect effect. The study of 3YF SKAP2 mutant suggests also that the glutamic acid-rich region plays a role since differences of signal between ΔLink2AB mutant and 3YF one was present for wildtype HCK and N1-fused Y501F HCK mutant. However, we cannot exclude a role of Y75F mutation. Similar analysis for SKAP1 (Fig 6A lower) shows different results. If activation of HCK has a similar effect than for SKAP2, the role of the SH2 domain was completely different. Similar level of interaction of R150A HCK mutant with either wildtype SKAP1 or with its ΔLink2AB mutant supports the view that no other motif than the SH2/ΔLink2AB module affects the signal. The difference of PPI signal between Y501F HCK mutant with either wildtype SKAP1 or ΔLink2AB mutant suggests that the HCK SH2 domain has an effect. Further analysis supports both a direct binding effect on SKAP1 and an indirect one due to others interactors, which modified the interaction between SKAP1 and HCK. Note that the results for ΔPHLink2 mutant is
completely different supporting that the PH domain and/or the interdomain Link2 play other role in this interaction.

Analyzing LCK data (Fig 6B) shows a completely different situation than for HCK data. For both SKAP2 (upper) and SKAP1 (lower), LCK activation decrease the PPI signal. This difference between the two SRC kinases are probably due to a difference of binding for the second module since DIMPH SKAP2 mutant has a similar pattern than wild-type SKAP2. A direct role of the SH2 domain is clearly seen for SKAP1 and also SKAP2 even if an indirect role is possible for both. Note in both cases, that a similar level of interaction of R154A LCK mutant with either wildtype SKAP1/SKAP2 or with their corresponding ΔLink2AB mutant supports the view that no other motif than the SH2/ΔLink2AB module affects the signal.

In conclusion, this analysis supports a strong increase of HCK binding due to its activation for both SKAP1 and SKAP2 proteins in contrast to LCK where the opposite is true. The SH2 domain of HCK plays a minor role for the interaction with SKAP2 in contrast to SKAP1. For LCK, its SH2 domain seems to play a similar role for both adaptors. Interestingly, study of DIMPH SKAP2 mutant suggests that the decrease of binding of LCK after its activation is due to the DIM domain. As previously mentioned, the study of 3YF SKAP2 mutant suggests that the two conserved motifs contain more than a unique site for SH2 domain.

## Study of the three FYN isoforms

The opposite effect of activation on the binding capacities between HCK and LCK kinases suggests studying the three FYN isoforms (FYNB using exon 7A, FYNT exon 7B, and FYN none of them). Fig 7 shows the result of this experiment using as reference N1-fused and C1-fused wild type FYN isoform. The pattern of wild type isoforms and their ΔLink2AB mutants is similar for SKAP1 and SKAP2 and will be described together. Wild type FYNT has a lower PPI signal in contrast to wild type FYNB, which has an upper signal. The ΔLink2AB mutants have a similar PPI signal than their wild type isoforms for C1 data in contrast to N1 data where this signal is lower. DIMPH SKAP2 mutant shows a similar pattern supporting than the PPI signal of the wild type isoform depends mainly of the second module and the DIM domain. For N1-fused data, an effect of the two conserved motifs is also detected. As previously reported, the ΔPHLink2 SKAP1 mutant has a different pattern than the ΔLink2AB SKAP1 mutant supporting than other factors than the two conserved motifs influence the PPI signal. The different binding capacity of these isoforms are due to exon 7 coding for the C-terminal part of the SH2 domain (fourteen last aa), the hinge, and the N-terminal part of the tyrosine-kinase catalytic domain (seventeen first aa).

## Discussion

The characterization of the SKAP1 modular architecture, which recognizes SRC kinases allowed us to determine that (1) two modules control these interactions; (2) each module has a different effect depending on SKAP1 and SKAP2 [27]. If a general pattern might be obtained for all SRC kinases, each one has also its own characteristics. This analysis was possible because the luciferase complementation test had a linear response over the entire range studied, irrespective of the significance threshold. This linearity was assessed by using the Differential Interaction Scatterplots and the scatterplot of barycenter. These results illustrate the interest of luciferase complementation assay to dissect how proteins of the same family interact with other proteins. Understanding how evolution shapes such interactions has important consequences.

By comparing the present results to previous ones on SKAP2 [27], we showed that both SKAP1 and SKAP2 modular architectures are globally similar, but exhibit differences mainly
completely different supporting that the PH domain and/or the interdomain Link2 play other role in this interaction.

Analyzing LCK data (Fig 6B) shows a completely different situation than for HCK data. For both SKAP2 (upper) and SKAP1 (lower), LCK activation decrease the PPI signal. This difference between the two SRC kinases are probably due to a difference of binding for the second module since DIMPH SKAP2 mutant has a similar pattern than wild-type SKAP2. A direct role of the SH2 domain is clearly seen for SKAP1 and also SKAP2 even if an indirect role is possible for both. Note in both cases, that a similar level of interaction of R154A LCK mutant with either wildtype SKAP1/SKAP2 or with their corresponding ΔLink2AB mutant supports the view that no other motif than the SH2/ΔLink2AB module affects the signal.

In conclusion, this analysis supports a strong increase of HCK binding due to its activation for both SKAP1 and SKAP2 proteins in contrast to LCK where the opposite is true. The SH2 domain of HCK plays a minor role for the interaction with SKAP2 in contrast to SKAP1. For LCK, its SH2 domain seems to play a similar role for both adaptors. Interestingly, study of DIMPH SKAP2 mutant suggests that the decrease of binding of LCK after its activation is due to the DIM domain. As previously mentioned, the study of 3YF SKAP2 mutant suggests that the two conserved motifs contain more than a unique site for SH2 domain.

## Study of the three FYN isoforms

The opposite effect of activation on the binding capacities between HCK and LCK kinases suggests studying the three FYN isoforms (FYNB using exon 7A, FYNT exon 7B, and FYN none of them). Fig 7 shows the result of this experiment using as reference N1-fused and C1-fused wild type FYN isoform. The pattern of wild type isoforms and their ΔLink2AB mutants is similar for SKAP1 and SKAP2 and will be described together. Wild type FYNT has a lower PPI signal in contrast to wild type FYNB, which has an upper signal. The ΔLink2AB mutants have a similar PPI signal than their wild type isoforms for C1 data in contrast to N1 data where this signal is lower. DIMPH SKAP2 mutant shows a similar pattern supporting than the PPI signal of the wild type isoform depends mainly of the second module and the DIM domain. For N1-fused data, an effect of the two conserved motifs is also detected. As previously reported, the ΔPHLink2 SKAP1 mutant has a different pattern than the ΔLink2AB SKAP1 mutant supporting than other factors than the two conserved motifs influence the PPI signal. The different binding capacity of these isoforms are due to exon 7 coding for the C-terminal part of the SH2 domain (fourteen last aa), the hinge, and the N-terminal part of the tyrosine-kinase catalytic domain (seventeen first aa).

## Discussion

The characterization of the SKAP1 modular architecture, which recognizes SRC kinases allowed us to determine that (1) two modules control these interactions; (2) each module has a different effect depending on SKAP1 and SKAP2 [27]. If a general pattern might be obtained for all SRC kinases, each one has also its own characteristics. This analysis was possible because the luciferase complementation test had a linear response over the entire range studied, irrespective of the significance threshold. This linearity was assessed by using the Differential Interaction Scatterplots and the scatterplot of barycenter. These results illustrate the interest of luciferase complementation assay to dissect how proteins of the same family interact with other proteins. Understanding how evolution shapes such interactions has important consequences.

By comparing the present results to previous ones on SKAP2 [27], we showed that both SKAP1 and SKAP2 modular architectures are globally similar, but exhibit differences mainly
<figure><img src="image_11.png" /><figcaption>**Fig 7.** The level of the PPI signal with SKAP1 and SKAP2 is lower for FYNT isoform than for FYNB isoform. Graphs show FYN normalized NLR of SKAP2 (left) and SKAP1 (right) for the three wild type FYN isoforms and their mutants. N1 data are normalized according to N1-fused wild type FYN and C1 data according to C1-fused wild type FYN. See Fig 6 caption for details.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g007

between C1-fused SRC kinases data and N1-fused ones as shown for both wildtype proteins and mutants (S4 Fig in S1 File, [27]). Refer also to this previous article to have an extensive discussion of the interest to study both localization of hemi-gaussia fusion for SRC kinases. The DIM domain of SKAP2 alone binds efficiently some C1-fused SRC kinases in contrast to SKAP1. In the present work, we confirmed these results showing that both SKAP1 modules, the DIM domain and the Link2 interdomain, are necessary to bind SRC kinases independently of the fusion localization. We also showed that the SH3 domain of SKAP1 does not have a global effect on these interactions in contrast to that of SKAP2 [27]. In both proteins, this SH3 domain has an indirect effect due to its binding capacity to interactors, which probably affects the localization of these adaptors and their interaction with SRC kinases. In the present paper, we showed also that each conserved motif located in the Link2 interdomain recognizes differently SRC kinases with a conserved pattern between SKAP1 and SKAP2 (Fig 5). A principal component analysis [65] confirmed that two different patterns due to either the SH3 domain or the two motifs in the Link2 interdomain, are detected on this C1-fused kinases data of SKAP2 (Supplemental data). A similar result was obtained on the C1-fused kinases data of SKAP1, supporting that the pattern due to the SH3 domain is not lost in SKAP1. It has not been detected only because its global effect is small or null. Note that this result explains that the pattern of ΔSH3 SKAP1 mutant is so different from that of the W333K SKAP1 mutant even if there global effect is similar (Fig 3). However, each deletion of one of the two conserved motifs located in Link2 interdomain proteins induces a similar pattern between N1-fused and C1-fused SRC kinases in both SKAP1 and SKAP2 (Fig 5 and principal component analysis). These results support that for both proteins each module recognizes differently each kinase independently to the other module and that these patterns are globally conserved since their duplication.

To find out more, search for conserved motifs in the interdomain Link2 common to both proteins allows to define two motifs composed of an SH2 binding site followed or forwards by a glutamic-acid rich region [70]. These two motifs play non-redundant and complementary role for binding to SRC kinases in both proteins (Fig 5, S4 Table in S2 File). Interestingly, each of these motifs recognizes similarly some SRC kinases (HCK, FYN, LCK, LYN, and SRMS) for both adaptors. There are some exceptions as BLK and FRK, which are specifically recognized
<figure><img src="image_10.png" /><figcaption>**Fig 7.** The level of the PPI signal with SKAP1 and SKAP2 is lower for FYNT isoform than for FYNB isoform. Graphs show FYN normalized NLR of SKAP2 (left) and SKAP1 (right) for the three wild type FYN isoforms and their mutants. N1 data are normalized according to N1-fused wild type FYN and C1 data according to C1-fused wild type FYN. See Fig 6 caption for details.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0296230.g007

between C1-fused SRC kinases data and N1-fused ones as shown for both wildtype proteins and mutants (S4 Fig in S1 File, [27]). Refer also to this previous article to have an extensive discussion of the interest to study both localization of hemi-gaussia fusion for SRC kinases. The DIM domain of SKAP2 alone binds efficiently some C1-fused SRC kinases in contrast to SKAP1. In the present work, we confirmed these results showing that both SKAP1 modules, the DIM domain and the Link2 interdomain, are necessary to bind SRC kinases independently of the fusion localization. We also showed that the SH3 domain of SKAP1 does not have a global effect on these interactions in contrast to that of SKAP2 [27]. In both proteins, this SH3 domain has an indirect effect due to its binding capacity to interactors, which probably affects the localization of these adaptors and their interaction with SRC kinases. In the present paper, we showed also that each conserved motif located in the Link2 interdomain recognizes differently SRC kinases with a conserved pattern between SKAP1 and SKAP2 (Fig 5). A principal component analysis [65] confirmed that two different patterns due to either the SH3 domain or the two motifs in the Link2 interdomain, are detected on this C1-fused kinases data of SKAP2 (Supplemental data). A similar result was obtained on the C1-fused kinases data of SKAP1, supporting that the pattern due to the SH3 domain is not lost in SKAP1. It has not been detected only because its global effect is small or null. Note that this result explains that the pattern of ΔSH3 SKAP1 mutant is so different from that of the W333K SKAP1 mutant even if there global effect is similar (Fig 3). However, each deletion of one of the two conserved motifs located in Link2 interdomain proteins induces a similar pattern between N1-fused and C1-fused SRC kinases in both SKAP1 and SKAP2 (Fig 5 and principal component analysis). These results support that for both proteins each module recognizes differently each kinase independently to the other module and that these patterns are globally conserved since their duplication.

To find out more, search for conserved motifs in the interdomain Link2 common to both proteins allows to define two motifs composed of an SH2 binding site followed or forwards by a glutamic-acid rich region [70]. These two motifs play non-redundant and complementary role for binding to SRC kinases in both proteins (Fig 5, S4 Table in S2 File). Interestingly, each of these motifs recognizes similarly some SRC kinases (HCK, FYN, LCK, LYN, and SRMS) for both adaptors. There are some exceptions as BLK and FRK, which are specifically recognized
by motif A of SKAP2. These results support that such evolutionary conservation of two motifs for both proteins is partially due to their non-redundant role. A second possible mechanism may be a cooperative effect between the two motifs. Interestingly, ΔLink2AB deletion induces a greater decrease of the binding to SRC kinases than the addition of ΔLink2A and ΔLink2B effect at least for N1-fused SRC kinases with both SKAP1 and SKAP2 (S5 Table in S2 File). Motifs recognized by SH2 domain have been subject to a lot of work such as those for SRC kinases [72,73]. The two motifs in Link2 interdomain are characteristic of a binding motif to their SH2 domain even if motif B is the more conformed. In our knowledge, a glutamic acid rich region has not been associated to such interactions. Note also, that other putative signals such as a Y based sorting signal or an Atg8 protein family binding site overlapped these two regions in which tyrosine residues should be unphosphorylated [70]. We cannot exclude that these motifs play an indirect role. Lastly, we showed that the interdomain Link1 affects these interactions only in the case of SKAP2 (Fig 3 and Table 3). Search for a conserved motif specific to SKAP2 was able to detect a conserved phosphorylation site for Casein kinase 1, which deletion increases the interaction with SRC kinases by a direct or indirect binding effect (Fig 4). We have to study the two mutations Serin to Alanin, and Serin to Prolin since the first classical one (S->A) preserves the α-helix predicted by I-Tasser. This analysis confirms that only the two serins have an effect without any effect of an α-helix if it exists. These results also support that for both proteins each motif located in Link2 interdomain recognizes differently each kinase independently to the other motif and that these patterns are globally conserved since their duplication.

The respective effect of the binding capacity of the SH2 domain and kinase activation was studied for HCK and LCK kinases interacting either with SKAP1 or SKAP2 (Fig 6). In most cases, a direct effect of the SH2 domain was detected except for SKAP2 interacting with the HCK kinase. In contrast, an opposite effect of kinase activation was observed for HCK and LCK kinases. Interestingly, kinase activation increased the binding of HCK kinase, which is mainly expressed in monocytes/macrophages lineage and neutrophils as SKAP2, in which the second module and the DIM domain play a higher role than in SKAP1. In contrast, kinase activation decreased the binding of LCK kinase, which is mainly expressed in T lymphocytes as SKAP1. This suggests that 1) the effect of kinase activation mainly depends on the second module and the DIM domain 2) a convergent evolution between SRC kinases and SKAP1 or SKAP2 might occur. The first part of the hypothesis is supported by previous data showing that the activation of some SRC kinases increases their binding to SKAP2 and to DIM domains even if the subcellular localization also plays a role [27]. A similar pattern between DIMPH SKAP2 mutant and wild type SKAP2 interacting with wild type LCK and LCK mutants (Fig 6) also supports this hypothesis. The hypothesis of convergent evolution suggests that the binding capacity of the DIM domain alone as in SKAP2 is associated with an increase of binding by kinase activation. This property is absent for kinases not interacting with SKAP2 as suggested by the study of LCK. Study of the three FYN isoforms supports this hypothesis (Fig 7): FYNT isoform expressed in T lymphocytes as SKAP1 has its binding to SKAP1 and SKAP2 decreasing compare to FYN isoform in contrast to FYNB isoform expressed in the brain as SKAP2. Study of ΔLink2AB mutants and of DIMPH SKAP2 mutant agrees with that hypothesis. Direct interactions between SKAP2 and HCK in monocytes/macrophages lineage [10] or FYN in the brain [9,74] have been reported in the literature as well as between SKAP1 and FYN [4,14,16] in T lymphocytes. However, this hypothesis has to be confirmed by further works.

In conclusion, this work dissected the SKAP1 modular organization recognizing SRC kinases defining two modules common also to SKAP2 modular organization but with different properties. This analysis was performed at the steady state where the knowledge of these interactions is particularly poor. The strength of the protein-protein interaction depends mainly of
by motif A of SKAP2. These results support that such evolutionary conservation of two motifs for both proteins is partially due to their non-redundant role. A second possible mechanism may be a cooperative effect between the two motifs. Interestingly, ΔLink2AB deletion induces a greater decrease of the binding to SRC kinases than the addition of ΔLink2A and ΔLink2B effect at least for N1-fused SRC kinases with both SKAP1 and SKAP2 (S5 Table in S2 File). Motifs recognized by SH2 domain have been subject to a lot of work such as those for SRC kinases [72,73]. The two motifs in Link2 interdomain are characteristic of a binding motif to their SH2 domain even if motif B is the more conformed. In our knowledge, a glutamic acid rich region has not been associated to such interactions. Note also, that other putative signals such as a Y based sorting signal or an Atg8 protein family binding site overlapped these two regions in which tyrosine residues should be unphosphorylated [70]. We cannot exclude that these motifs play an indirect role. Lastly, we showed that the interdomain Link1 affects these interactions only in the case of SKAP2 (Fig 3 and Table 3). Search for a conserved motif specific to SKAP2 was able to detect a conserved phosphorylation site for Casein kinase 1, which deletion increases the interaction with SRC kinases by a direct or indirect binding effect (Fig 4). We have to study the two mutations Serin to Alanin, and Serin to Prolin since the first classical one (S->A) preserves the α-helix predicted by I-Tasser. This analysis confirms that only the two serins have an effect without any effect of an α-helix if it exists. These results also support that for both proteins each motif located in Link2 interdomain recognizes differently each kinase independently to the other motif and that these patterns are globally conserved since their duplication.

The respective effect of the binding capacity of the SH2 domain and kinase activation was studied for HCK and LCK kinases interacting either with SKAP1 or SKAP2 (Fig 6). In most cases, a direct effect of the SH2 domain was detected except for SKAP2 interacting with the HCK kinase. In contrast, an opposite effect of kinase activation was observed for HCK and LCK kinases. Interestingly, kinase activation increased the binding of HCK kinase, which is mainly expressed in monocytes/macrophages lineage and neutrophils as SKAP2, in which the second module and the DIM domain play a higher role than in SKAP1. In contrast, kinase activation decreased the binding of LCK kinase, which is mainly expressed in T lymphocytes as SKAP1. This suggests that 1) the effect of kinase activation mainly depends on the second module and the DIM domain 2) a convergent evolution between SRC kinases and SKAP1 or SKAP2 might occur. The first part of the hypothesis is supported by previous data showing that the activation of some SRC kinases increases their binding to SKAP2 and to DIM domains even if the subcellular localization also plays a role [27]. A similar pattern between DIMPH SKAP2 mutant and wild type SKAP2 interacting with wild type LCK and LCK mutants (Fig 6) also supports this hypothesis. The hypothesis of convergent evolution suggests that the binding capacity of the DIM domain alone as in SKAP2 is associated with an increase of binding by kinase activation. This property is absent for kinases not interacting with SKAP2 as suggested by the study of LCK. Study of the three FYN isoforms supports this hypothesis (Fig 7): FYNT isoform expressed in T lymphocytes as SKAP1 has its binding to SKAP1 and SKAP2 decreasing compare to FYN isoform in contrast to FYNB isoform expressed in the brain as SKAP2. Study of ΔLink2AB mutants and of DIMPH SKAP2 mutant agrees with that hypothesis. Direct interactions between SKAP2 and HCK in monocytes/macrophages lineage [10] or FYN in the brain [9,74] have been reported in the literature as well as between SKAP1 and FYN [4,14,16] in T lymphocytes. However, this hypothesis has to be confirmed by further works.

In conclusion, this work dissected the SKAP1 modular organization recognizing SRC kinases defining two modules common also to SKAP2 modular organization but with different properties. This analysis was performed at the steady state where the knowledge of these interactions is particularly poor. The strength of the protein-protein interaction depends mainly of
properties of SRC kinases and not of those of SKAP1/SKAP2. Interestingly, the two SRC kinases, HCK and FYNB, known to interact with SKAP2 have a higher strength of PPI than the SRC kinases that do not interact with this protein (LCK, and FYNT). This result suggests that this property is sufficient to explain what happens later during cell activation. As this is not the case for SKAP1, other events must occur, specific or not to T lymphocytes. Comparison with that of SKAP2 allowed us to define the role of different conserved motifs in these interactions, which had interesting consequences on how evolution shapes these protein-protein interaction. Lastly, study of HCK and LCK kinases as well as FYN isoforms suggests an interesting hypothesis of coevolution between the properties of these kinases and those of SKAP1 and SKAP2 adaptors.

**Supporting information**

**S1 Raw data.**

(XLSX)

**S1 File. A docx file containing the four supplemental figures and the principal component analysis.**

(DOCX)

**S2 File.** An excel file containing S1 Table Sequence of new primers used in this study; S2 Table List of the 38 species with a SKAP2 protein sequence; S3 Table Role of the different domains and interdomains of SKAP1 and SKAP2 in the interaction with SRC kinases; S4 Table Statistical analysis of deletion of motif A, and B alone or in combination in SKAP1 (upper part) and SKAP2 (lower part); S5 Table Search for an cooperative effect between motif A and B.

(XLSX)

**Acknowledgments**

We thank Patricia Cassonnet, and Mélanie Dos Santos for technical supports, Chen Kuang-Yu, Tim Krischuns and Guillaume Dugied for regularly providing HEK 293T cell culture, and the Dana-Farber Cancer Institute (Boston, MA, USA) for providing us with ORFs.

**Author Contributions**

**Conceptualization:** Jean-François Bureau.

**Formal analysis:** Jean-François Bureau.

**Funding acquisition:** Anavaj Sakuntabhai.

**Investigation:** Laurine Levillayer, Camille Brighelli, Jean-François Bureau.

**Writing – original draft:** Jean-François Bureau.

**Writing – review & editing:** Caroline Demeret.

**References**

1. Abbasi AA. Diversification of four human HOX gene clusters by step-wise evolution rather than ancient whole-genome duplications. Dev Genes Evol. 2015; 225(6):353–7. Epub 2015/10/21. https://doi.org/10.1007/s00427-015-0518-z PMID: 26481129.

2. Kliche S, Breitling D, Togni M, Pusch R, Heuer K, Wang X, et al. The ADAP/SKAP55 signaling module regulates T-cell receptor-mediated integrin activation through plasma membrane targeting of Rap1. Mol
properties of SRC kinases and not of those of SKAP1/SKAP2. Interestingly, the two SRC kinases, HCK and FYNB, known to interact with SKAP2 have a higher strength of PPI than the SRC kinases that do not interact with this protein (LCK, and FYNT). This result suggests that this property is sufficient to explain what happens later during cell activation. As this is not the case for SKAP1, other events must occur, specific or not to T lymphocytes. Comparison with that of SKAP2 allowed us to define the role of different conserved motifs in these interactions, which had interesting consequences on how evolution shapes these protein-protein interaction. Lastly, study of HCK and LCK kinases as well as FYN isoforms suggests an interesting hypothesis of coevolution between the properties of these kinases and those of SKAP1 and SKAP2 adaptors.

**Supporting information**

**S1 Raw data.**

(XLSX)

**S1 File. A docx file containing the four supplemental figures and the principal component analysis.**

(DOCX)

**S2 File.** An excel file containing S1 Table Sequence of new primers used in this study; S2 Table List of the 38 species with a SKAP2 protein sequence; S3 Table Role of the different domains and interdomains of SKAP1 and SKAP2 in the interaction with SRC kinases; S4 Table Statistical analysis of deletion of motif A, and B alone or in combination in SKAP1 (upper part) and SKAP2 (lower part); S5 Table Search for an cooperative effect between motif A and B.

(XLSX)

**Acknowledgments**

We thank Patricia Cassonnet, and Mélanie Dos Santos for technical supports, Chen Kuang-Yu, Tim Krischuns and Guillaume Dugied for regularly providing HEK 293T cell culture, and the Dana-Farber Cancer Institute (Boston, MA, USA) for providing us with ORFs.

**Author Contributions**

**Conceptualization:** Jean-François Bureau.

**Formal analysis:** Jean-François Bureau.

**Funding acquisition:** Anavaj Sakuntabhai.

**Investigation:** Laurine Levillayer, Camille Brighelli, Jean-François Bureau.

**Writing – original draft:** Jean-François Bureau.

**Writing – review & editing:** Caroline Demeret.

**References**

1. Abbasi AA. Diversification of four human HOX gene clusters by step-wise evolution rather than ancient whole-genome duplications. Dev Genes Evol. 2015; 225(6):353-7. Epub 2015/10/21. https://doi.org/10.1007/s00427-015-0518-z PMID: 26481129.

2. Kliche S, Breitling D, Togni M, Pusch R, Heuer K, Wang X, et al. The ADAP/SKAP55 signaling module regulates T-cell receptor-mediated integrin activation through plasma membrane targeting of Rap1. Mol
Cell Biol. 2006; 26(19):7130-44. https://doi.org/10.1128/MCB.00331-06 PMID: 16980616; PubMed Central PMCID: PMC1592884.

3. Timms JF, Swanson KD, Marie-Cardine A, Raab M, Rudd CE, Schraven B, et al. SHPS-1 is a scaffold for assembling distinct adhesion-regulated multi-protein complexes in macrophages. Curr Biol. 1999; 9 (16):927-30. https://doi.org/10.1016/s0960-9822(99)80401-1 PMID: 10469599.

4. Marie-Cardine A, Verhagen AM, Eckerskorn C, Schraven B. SKAP-HOM, a novel adaptor protein homologous to the FYN-associated protein SKAP55. FEBS Lett. 1998; 435(1):55-60. https://doi.org/10.1016/s0014-5793(98)01040-0 PMID: 9755858.

5. Alenghat FJ, Baca QJ, Rubin NT, Pao LI, Matozaki T, Lowell CA, et al. Macrophages require Skap2 and Siralpha for integrin-stimulated cytoskeletal rearrangement. J Cell Sci. 2012; 125(Pt 22):5535-45. https://doi.org/10.1242/jcs.111260 PMID: 22976304; PubMed Central PMCID: PMC3561861.

6. Nguyen GT, Shaban L, Mack M, Swanson KD, Bunnell SC, Sykes DB, et al. SKAP2 is required for defense against K. pneumoniae infection and neutrophil respiratory burst. Elife. 2020; 9. Epub 2020/05/01. https://doi.org/10.7554/eLife.56656 PMID: 32352382; PubMed Central PMCID: PMC7250567.

7. Nguyen GT, Xu S, Adams W, Leong JM, Bunnell SC, Mansour MK, et al. Neutrophils require SKAP2 for reactive oxygen species production following C-type lectin and Candida stimulation. iScience. 2021; 24 (8):102871. Epub 2021/08/14. https://doi.org/10.1016/j.isci.2021.102871 PMID: 34386732; PubMed Central PMCID: PMC8346660.

8. Shaban L, Nguyen GT, Mecsas-Faxon BD, Swanson KD, Tan S, Mecsas J. Yersinia pseudotuberculosis YopH targets SKAP2-dependent and independent signaling pathways to block neutrophil antimicrobial mechanisms during infection. PLoS Pathog. 2020; 16(5):e1008576. Epub 2020/05/12. https://doi.org/10.1371/journal.ppat.1008576 PMID: 32392230; PubMed Central PMCID: PMC7241846.

9. Ghelman J, Grewing L, Windener F, Albrecht S, Zarbock A, Kuhlmann T. SKAP2 as a new regulator of oligodendroglial migration and myelin sheath formation. Glia. 2021; 69(11):2699-716. Epub 2021/07/30. https://doi.org/10.1002/glia.24066 PMID: 34324225.

10. Kouroku Y, Soyama A, Fujita E, Urase K, Tsukahara T, Momoi T. RA70 is a src kinase-associated protein expressed ubiquitously. Biochem Biophys Res Commun. 1998; 252(3):738-42. https://doi.org/10.1006/bbrc.1998.9637 PMID: 9837776.

11. Raab M, Strebhardt K, Rudd CE. Immune adaptor protein SKAP1 (SKAP-55) forms homodimers as mediated by the N-terminal region. BMC Res Notes. 2018; 11(1):869. Epub 2018/12/14. https://doi.org/10.1186/s13104-018-3976-3 PMID: 30522503; PubMed Central PMCID: PMC6282339.

12. Swanson KD, Tang Y, Ceccarelli DF, Poy F, Sliwa JP, Neel BG, et al. The Skap-hom dimerization and PH domains comprise a 3'-phosphoinositide-gated molecular switch. Mol Cell. 2008; 32(4):564-75. https://doi.org/10.1016/j.molcel.2008.09.022 PMID: 19026786; PubMed Central PMCID: PMC2628593.

13. Witte A, Meineke B, Sticht J, Philipsen L, Kuropka B, Muller AJ, et al. D120 and K152 within the PH Domain of T Cell Adapter SKAP55 Regulate Plasma Membrane Targeting of SKAP55 and LFA-1 Affinity Modulation in Human T Lymphocytes. Mol Cell Biol. 2017; 37(7). Epub 2017/01/06. https://doi.org/10.1128/MCB.00509-16 PMID: 28052935; PubMed Central PMCID: PMC5359422.

14. Wu L, Yu Z, Shen SH. SKAP55 recruits to lipid rafts and positively mediates the MAPK pathway upon T cell receptor activation. J Biol Chem. 2002; 277(43):40420-7. Epub 2002/08/13. https://doi.org/10.1074/jbc.M206023200 PMID: 12171928.

15. Moog-Lutz C, Peterson EJ, Lutz PG, Eliason S, Cave-Riant F, Singer A, et al. PRAM-1 is a novel adaptor protein regulated by retinoic acid (RA) and promyelocytic leukemia (PML)-RA receptor alpha in acute promyelocytic leukemia cells. J Biol Chem. 2001; 276(25):22375-81. Epub 2001/04/13. https://doi.org/10.1074/jbc.M011683200 PMID: 11301322.

16. Liu J, Kang H, Raab M, da Silva AJ, Kraeft SK, Rudd CE. FYB (FYN binding protein) serves as a binding partner for lymphoid protein and FYN kinase substrate SKAP55 and a SKAP55-related protein in T cells. Proc Natl Acad Sci U S A. 1998; 95(15):8779-84. https://doi.org/10.1073/pnas.95.15.8779 PMID: 9671755; PubMed Central PMCID: PMC21153.

17. Bureau JF, Cassonnet P, Grange L, Dessapt J, Jones L, Demeret C, et al. The SRC-family tyrosine kinase HCK shapes the landscape of SKAP2 interactome. Oncotarget. 2018; 9(17):13102-15. Epub 2018/03/24. https://doi.org/10.18632/oncotarget.24424 PMID: 29568343; PubMed Central PMCID: PMC5862564.

18. Rolland T, Tasan M, Charoteaux B, Pevzner SJ, Zhong Q, Sahni N, et al. A proteome-scale map of the human interactome network. Cell. 2014; 159(5):1212-26. https://doi.org/10.1016/j.cell.2014.10.050 PMID: 25416956; PubMed Central PMCID: PMC4266588.

19. Ophir MJ, Liu BC, Bunnell SC. The N terminus of SKAP55 enables T cell adhesion to TCR and integrin ligands via distinct mechanisms. J Cell Biol. 2013; 203(6):1021-41. https://doi.org/10.1083/jcb.201305088 PMID: 24368808; PubMed Central PMCID: PMC3871428.
Cell Biol. 2006; 26(19):7130-44. https://doi.org/10.1128/MCB.00331-06 PMID: 16980616; PubMed Central PMCID: PMC1592884.

3. Timms JF, Swanson KD, Marie-Cardine A, Raab M, Rudd CE, Schraven B, et al. SHPS-1 is a scaffold for assembling distinct adhesion-regulated multi-protein complexes in macrophages. Curr Biol. 1999; 9 (16):927-30. https://doi.org/10.1016/s0960-9822(99)80401-1 PMID: 10469599.

4. Marie-Cardine A, Verhagen AM, Eckerskorn C, Schraven B. SKAP-HOM, a novel adaptor protein homologous to the FYN-associated protein SKAP55. FEBS Lett. 1998; 435(1):55-60. https://doi.org/10.1016/s0014-5793(98)01040-0 PMID: 9755858.

5. Alenghat FJ, Baca QJ, Rubin NT, Pao LI, Matozaki T, Lowell CA, et al. Macrophages require Skap2 and Siralpha for integrin-stimulated cytoskeletal rearrangement. J Cell Sci. 2012; 125(Pt 22):5535-45. https://doi.org/10.1242/jcs.111260 PMID: 22976304; PubMed Central PMCID: PMC3561861.

6. Nguyen GT, Shaban L, Mack M, Swanson KD, Bunnell SC, Sykes DB, et al. SKAP2 is required for defense against K. pneumoniae infection and neutrophil respiratory burst. Elife. 2020; 9. Epub 2020/05/01. https://doi.org/10.7554/eLife.56656 PMID: 32352382; PubMed Central PMCID: PMC7250567.

7. Nguyen GT, Xu S, Adams W, Leong JM, Bunnell SC, Mansour MK, et al. Neutrophils require SKAP2 for reactive oxygen species production following C-type lectin and Candida stimulation. iScience. 2021; 24 (8):102871. Epub 2021/08/14. https://doi.org/10.1016/j.isci.2021.102871 PMID: 34386732; PubMed Central PMCID: PMC8346660.

8. Shaban L, Nguyen GT, Mecsas-Faxon BD, Swanson KD, Tan S, Mecsas J. Yersinia pseudotuberculosis YopH targets SKAP2-dependent and independent signaling pathways to block neutrophil antimicrobial mechanisms during infection. PLoS Pathog. 2020; 16(5):e1008576. Epub 2020/05/12. https://doi.org/10.1371/journal.ppat.1008576 PMID: 32392230; PubMed Central PMCID: PMC7241846.

9. Ghelman J, Grewing L, Windener F, Albrecht S, Zarbock A, Kuhlmann T. SKAP2 as a new regulator of oligodendroglial migration and myelin sheath formation. Glia. 2021; 69(11):2699-716. Epub 2021/07/30. https://doi.org/10.1002/glia.24066 PMID: 34324225.

10. Kouroku Y, Soyama A, Fujita E, Urase K, Tsukahara T, Momoi T. RA70 is a src kinase-associated protein expressed ubiquitously. Biochem Biophys Res Commun. 1998; 252(3):738-42. https://doi.org/10.1006/bbrc.1998.9637 PMID: 9837776.

11. Raab M, Strebhardt K, Rudd CE. Immune adaptor protein SKAP1 (SKAP-55) forms homodimers as mediated by the N-terminal region. BMC Res Notes. 2018; 11(1):869. Epub 2018/12/14. https://doi.org/10.1186/s13104-018-3976-3 PMID: 30522503; PubMed Central PMCID: PMC6282339.

12. Swanson KD, Tang Y, Ceccarelli DF, Poy F, Sliwa JP, Neel BG, et al. The Skap-hom dimerization and PH domains comprise a 3'-phosphoinositide-gated molecular switch. Mol Cell. 2008; 32(4):564-75. https://doi.org/10.1016/j.molcel.2008.09.022 PMID: 19026786; PubMed Central PMCID: PMC2628593.

13. Witte A, Meineke B, Sticht J, Philipsen L, Kuropka B, Muller AJ, et al. D120 and K152 within the PH Domain of T Cell Adapter SKAP55 Regulate Plasma Membrane Targeting of SKAP55 and LFA-1 Affinity Modulation in Human T Lymphocytes. Mol Cell Biol. 2017; 37(7). Epub 2017/01/06. https://doi.org/10.1128/MCB.00509-16 PMID: 28052935; PubMed Central PMCID: PMC5359422.

14. Wu L, Yu Z, Shen SH. SKAP55 recruits to lipid rafts and positively mediates the MAPK pathway upon T cell receptor activation. J Biol Chem. 2002; 277(43):40420-7. Epub 2002/08/13. https://doi.org/10.1074/jbc.M206023200 PMID: 12171928.

15. Moog-Lutz C, Peterson EJ, Lutz PG, Eliason S, Cave-Riant F, Singer A, et al. PRAM-1 is a novel adaptor protein regulated by retinoic acid (RA) and promyelocytic leukemia (PML)-RA receptor alpha in acute promyelocytic leukemia cells. J Biol Chem. 2001; 276(25):22375-81. Epub 2001/04/13. https://doi.org/10.1074/jbc.M011683200 PMID: 11301322.

16. Liu J, Kang H, Raab M, da Silva AJ, Kraeft SK, Rudd CE. FYB (FYN binding protein) serves as a binding partner for lymphoid protein and FYN kinase substrate SKAP55 and a SKAP55-related protein in T cells. Proc Natl Acad Sci U S A. 1998; 95(15):8779-84. https://doi.org/10.1073/pnas.95.15.8779 PMID: 9671755; PubMed Central PMCID: PMC21153.

17. Bureau JF, Cassonnet P, Grange L, Dessapt J, Jones L, Demeret C, et al. The SRC-family tyrosine kinase HCK shapes the landscape of SKAP2 interactome. Oncotarget. 2018; 9(17):13102-15. Epub 2018/03/24. https://doi.org/10.18632/oncotarget.24424 PMID: 29568343; PubMed Central PMCID: PMC5862564.

18. Rolland T, Tasan M, Charoteaux B, Pevzner SJ, Zhong Q, Sahni N, et al. A proteome-scale map of the human interactome network. Cell. 2014; 159(5):1212-26. https://doi.org/10.1016/j.cell.2014.10.050 PMID: 25416956; PubMed Central PMCID: PMC4266588.

19. Ophir MJ, Liu BC, Bunnell SC. The N terminus of SKAP55 enables T cell adhesion to TCR and integrin ligands via distinct mechanisms. J Cell Biol. 2013; 203(6):1021-41. https://doi.org/10.1083/jcb.201305088 PMID: 24368808; PubMed Central PMCID: PMC3871428.
20. Raab M, Lu Y, Kohler K, Smith X, Strebhardt K, Rudd CE. LFA-1 activates focal adhesion kinases FAK1/PYK2 to generate LAT-GRB2-SKAP1 complexes that terminate T-cell conjugate formation. Nat Commun. 2017; 8:16001. Epub 2017/07/13. https://doi.org/10.1038/ncomms16001 PMID: 28699640; PubMed Central PMCID: PMC5510181.

21. Wang H, Rudd CE. SKAP-55, SKAP-55-related and ADAP adaptors modulate integrin-mediated immune-cell adhesion. Trends Cell Biol. 2008; 18(10):486-93. https://doi.org/10.1016/j.tcb.2008.07.005 PMID: 18760924; PubMed Central PMCID: PMC3512129.

22. Kliche S, Worbs T, Wang X, Degen J, Patzak I, Meineke B, et al. CCR7-mediated LFA-1 functions in T cells are regulated by 2 independent ADAP/SKAP55 modules. Blood. 2012; 119(3):777-85. Epub 2011/11/26. https://doi.org/10.1182/blood-2011-06-362269 PMID: 22117043.

23. Raab M, Wang H, Lu Y, Smith X, Wu Z, Strebhardt K, et al. T cell receptor "inside-out" pathway via signaling module SKAP1-RapL regulates T cell motility and interactions in lymph nodes. Immunity. 2010; 32(4):541-56. Epub 2010/03/30. https://doi.org/10.1016/j.immuni.2010.03.007 PMID: 20346707; PubMed Central PMCID: PMC3812847.

24. Raab M, Strebhardt K, Rudd CE. Immune adaptor SKAP1 acts a scaffold for Polo-like kinase 1 (PLK1) for the optimal cell cycling of T-cells. Sci Rep. 2019; 9(1):10462. Epub 2019/07/20. https://doi.org/10.1038/s41598-019-45627-9 PMID: 31320682; PubMed Central PMCID: PMC6639320.

25. Liu C, Raab M, Gui Y, Rudd CE. Multi-functional adaptor SKAP1: regulator of integrin activation, the stop-signal, and the proliferation of T cells. Front Immunol. 2023; 14:1192838. Epub 2023/06/16. https://doi.org/10.3389/fimmu.2023.1192838 PMID: 37325633; PubMed Central PMCID: PMC10264576.

26. Fujii Y, Wakahara S, Nakao T, Hara T, Ohtake H, Komurasaki T, et al. Targeting of MIST to Src-family kinases via SKAP55-SLAP-130 adaptor complex in mast cells. FEBS Lett. 2003; 540(1-3):111-6. Epub 2003/04/12. https://doi.org/10.1016/s0014-5793(03)00234-5 PMID: 12681493.

27. Levillayer L, Cassonnet P, Declercq M, Santos MD, Lebreton L, Danezi K, et al. SKAP2 Modular Organization Differently Recognizes SRC Kinases Depending on Their Activation Status and Localization. Mol Cell Proteomics. 2023; 22(1):100451. Epub 2022/11/25. https://doi.org/10.1016/j.mcpro.2022.100451 PMID: 36423812; PubMed Central PMCID: PMC9792355.

28. Ingley E. Src family kinases: regulation of their activities, levels and identification of new pathways. Biochim Biophys Acta. 2008; 1784(1):56-65. https://doi.org/10.1016/j.bbpap.2007.08.012 PMID: 17905674.

29. Robinson DR, Wu YM, Lin SF. The protein tyrosine kinase family of the human genome. Oncogene. 2000; 19(49)(49):5548-57. PubMed PMID: CCC:000165659300001. https://doi.org/10.1038/sj.onc.1203957 PMID: 11114734

30. Roskoski R, Jr., Src protein-tyrosine kinase structure, mechanism, and small molecule inhibitors. Pharmacological research. 2015; 94:9-25. PubMed PMID: Medline: https://doi.org/10.1016/j.phrs.2015.01.003 PMID: 25662515.

31. Kornev AP, Taylor SS. Dynamics-Driven Allostery in Protein Kinases. Trends Biochem Sci. 2015; 40 (11):628-47. https://doi.org/10.1016/j.tibs.2015.09.002 PMID: 26481499; PubMed Central PMCID: PMC4630092.

32. Madan LK, Welsh CL, Kornev AP, Taylor SS. The "violin model": Looking at community networks for dynamic allostery. J Chem Phys. 2023; 158(8)(8)-. PubMed PMID: CCC:000937824600017. https://doi.org/10.1063/5.0138175 PMID: 36859094

33. Meng Y, Pond MP, Roux B. Tyrosine Kinase Activation and Conformational Flexibility: Lessons from Src-Family Tyrosine Kinases. Acc Chem Res. 2017; 50(5):1193-201. Epub 2017/04/21. https://doi.org/10.1021/acs.accounts.7b00012 PMID: 28426203; PubMed Central PMCID: PMC6602532.

34. Cooper JA, MacAuley A. Potential positive and negative autoregulation of p60c-src by intermolecular autophosphorylation. Proc Natl Acad Sci U S A. 1988; 85(12):4232-6. https://doi.org/10.1073/pnas.85.12.4232 PMID: 2454466; PubMed Central PMCID: PMC280401.

35. Patwardhan P, Resh MD. Myristoylation and membrane binding regulate c-Src stability and kinase activity. Mol Cell Biol. 2010; 30(17):4094-107. Epub 2010/06/30. https://doi.org/10.1128/MCB.00246-10 PMID: 20584982; PubMed Central PMCID: PMC2937550.

36. Sato I, Obata Y, Kasahara K, Nakayama Y, Fukumoto Y, Yamasaki T, et al. Differential trafficking of Src, Lyn, Yes and Fyn is specified by the state of palmitoylation in the SH4 domain. J Cell Sci. 2009; 122(Pt 7):965-75. https://doi.org/10.1242/jcs.034843 PMID: 19258394.

37. Aakre CD, Herrou J, Phung TN, Perchuk BS, Crosson S, Laub MT. Evolving new protein-protein interaction specificity through promiscuous intermediates. Cell. 2015; 163(3):594-606. Epub 2015/10/20. https://doi.org/10.1016/j.cell.2015.09.055 PMID: 26478181; PubMed Central PMCID: PMC4623991.

38. Andreani J, Guerois R. Evolution of protein interactions: from interactomes to interfaces. Arch Biochem Biophys. 2014; 554:65-75. Epub 2014/05/24. https://doi.org/10.1016/j.abb.2014.05.010 PMID: 24853495.
20. Raab M, Lu Y, Kohler K, Smith X, Strebhardt K, Rudd CE. LFA-1 activates focal adhesion kinases FAK1/PYK2 to generate LAT-GRB2-SKAP1 complexes that terminate T-cell conjugate formation. Nat Commun. 2017; 8:16001. Epub 2017/07/13. https://doi.org/10.1038/ncomms16001 PMID: 28699640; PubMed Central PMCID: PMC5510181.

21. Wang H, Rudd CE. SKAP-55, SKAP-55-related and ADAP adaptors modulate integrin-mediated immune-cell adhesion. Trends Cell Biol. 2008; 18(10):486-93. https://doi.org/10.1016/j.tcb.2008.07.005 PMID: 18760924; PubMed Central PMCID: PMC3512129.

22. Kliche S, Worbs T, Wang X, Degen J, Patzak I, Meineke B, et al. CCR7-mediated LFA-1 functions in T cells are regulated by 2 independent ADAP/SKAP55 modules. Blood. 2012; 119(3):777-85. Epub 2011/11/26. https://doi.org/10.1182/blood-2011-06-362269 PMID: 22117043.

23. Raab M, Wang H, Lu Y, Smith X, Wu Z, Strebhardt K, et al. T cell receptor "inside-out" pathway via signaling module SKAP1-RapL regulates T cell motility and interactions in lymph nodes. Immunity. 2010; 32(4):541-56. Epub 2010/03/30. https://doi.org/10.1016/j.immuni.2010.03.007 PMID: 20346707; PubMed Central PMCID: PMC3812847.

24. Raab M, Strebhardt K, Rudd CE. Immune adaptor SKAP1 acts a scaffold for Polo-like kinase 1 (PLK1) for the optimal cell cycling of T-cells. Sci Rep. 2019; 9(1):10462. Epub 2019/07/20. https://doi.org/10.1038/s41598-019-45627-9 PMID: 31320682; PubMed Central PMCID: PMC6639320.

25. Liu C, Raab M, Gui Y, Rudd CE. Multi-functional adaptor SKAP1: regulator of integrin activation, the stop-signal, and the proliferation of T cells. Front Immunol. 2023; 14:1192838. Epub 2023/06/16. https://doi.org/10.3389/fimmu.2023.1192838 PMID: 37325633; PubMed Central PMCID: PMC10264576.

26. Fujii Y, Wakahara S, Nakao T, Hara T, Ohtake H, Komurasaki T, et al. Targeting of MIST to Src-family kinases via SKAP55-SLAP-130 adaptor complex in mast cells. FEBS Lett. 2003; 540(1-3):111-6. Epub 2003/04/12. https://doi.org/10.1016/s0014-5793(03)00234-5 PMID: 12681493.

27. Levillayer L, Cassonnet P, Declercq M, Santos MD, Lebreton L, Danezi K, et al. SKAP2 Modular Organization Differently Recognizes SRC Kinases Depending on Their Activation Status and Localization. Mol Cell Proteomics. 2023; 22(1):100451. Epub 2022/11/25. https://doi.org/10.1016/j.mcpro.2022.100451 PMID: 36423812; PubMed Central PMCID: PMC9792355.

28. Ingley E. Src family kinases: regulation of their activities, levels and identification of new pathways. Biochim Biophys Acta. 2008; 1784(1):56-65. https://doi.org/10.1016/j.bbpap.2007.08.012 PMID: 17905674.

29. Robinson DR, Wu YM, Lin SF. The protein tyrosine kinase family of the human genome. Oncogene. 2000; 19(49)(49):5548-57. PubMed PMID: CCC:000165659300001. https://doi.org/10.1038/sj.onc.1203957 PMID: 11114734

30. Roskoski R, Jr., Src protein-tyrosine kinase structure, mechanism, and small molecule inhibitors. Pharmacological research. 2015; 94:9-25. PubMed PMID: Medline: https://doi.org/10.1016/j.phrs.2015.01.003 PMID: 25662515.

31. Kornev AP, Taylor SS. Dynamics-Driven Allostery in Protein Kinases. Trends Biochem Sci. 2015; 40 (11):628-47. https://doi.org/10.1016/j.tibs.2015.09.002 PMID: 26481499; PubMed Central PMCID: PMC4630092.

32. Madan LK, Welsh CL, Kornev AP, Taylor SS. The "violin model": Looking at community networks for dynamic allostery. J Chem Phys. 2023; 158(8)(8):-. PubMed PMID: CCC:000937824600017. https://doi.org/10.1063/5.0138175 PMID: 36859094

33. Meng Y, Pond MP, Roux B. Tyrosine Kinase Activation and Conformational Flexibility: Lessons from Src-Family Tyrosine Kinases. Acc Chem Res. 2017; 50(5):1193-201. Epub 2017/04/21. https://doi.org/10.1021/acs.accounts.7b00012 PMID: 28426203; PubMed Central PMCID: PMC6602532.

34. Cooper JA, MacAuley A. Potential positive and negative autoregulation of p60c-src by intermolecular autophosphorylation. Proc Natl Acad Sci U S A. 1988; 85(12):4232-6. https://doi.org/10.1073/pnas.85.12.4232 PMID: 2454466; PubMed Central PMCID: PMC280401.

35. Patwardhan P, Resh MD. Myristoylation and membrane binding regulate c-Src stability and kinase activity. Mol Cell Biol. 2010; 30(17):4094-107. Epub 2010/06/30. https://doi.org/10.1128/MCB.00246-10 PMID: 20584982; PubMed Central PMCID: PMC2937550.

36. Sato I, Obata Y, Kasahara K, Nakayama Y, Fukumoto Y, Yamasaki T, et al. Differential trafficking of Src, Lyn, Yes and Fyn is specified by the state of palmitoylation in the SH4 domain. J Cell Sci. 2009; 122(Pt 7):965-75. https://doi.org/10.1242/jcs.034843 PMID: 19258394.

37. Aakre CD, Herrou J, Phung TN, Perchuk BS, Crosson S, Laub MT. Evolving new protein-protein interaction specificity through promiscuous intermediates. Cell. 2015; 163(3):594-606. Epub 2015/10/20. https://doi.org/10.1016/j.cell.2015.09.055 PMID: 26478181; PubMed Central PMCID: PMC4623991.

38. Andreani J, Guerois R. Evolution of protein interactions: from interactomes to interfaces. Arch Biochem Biophys. 2014; 554:65-75. Epub 2014/05/24. https://doi.org/10.1016/j.abb.2014.05.010 PMID: 24853495.
39. Defoort J, Van de Peer Y, Carretero-Paulet L. The Evolution of Gene Duplicates in Angiosperms and the Impact of Protein-Protein Interactions and the Mechanism of Duplication. Genome Biol Evol. 2019; 11(8):2292–305. Epub 2019/08/01. https://doi.org/10.1093/gbe/evz156 PMID: 31364708; PubMed Central PMCID: PMC6735927.

40. Levy ED, Pereira-Leal JB. Evolution and dynamics of protein interactions and networks. Curr Opin Struct Biol. 2008; 18(3):349–57. Epub 2008/05/02. https://doi.org/10.1016/j.sbi.2008.03.003 PMID: 18448325.

41. Nooren IM, Thornton JM. Diversity of protein-protein interactions. EMBO J. 2003; 22(14):3486–92. Epub 2003/07/11. https://doi.org/10.1093/emboj/cdg359 PMID: 12853464; PubMed Central PMCID: PMC165629.

42. Nooren IMA, Thornton JM. Structural characterisation and functional significance of transient protein-protein interactions. Journal of Molecular Biology. 2003; 325(5)(5):991–1018. PubMed PMID: CCC:000180676300011. https://doi.org/10.1016/s0022-2836(02)01281-0 PMID: 12527304

43. Cafarelli TM, Desbuleux A, Wang Y, Choi SG, De Ridder D, Vidal M. Mapping, modeling, and characterization of protein-protein interactions on a proteomic scale. Curr Opin Struct Biol. 2017; 44:201–10. Epub 2017/06/03. https://doi.org/10.1016/j.sbi.2017.05.003 PMID: 28575754.

44. Ellis JD, Barrios-Rodiles M, Colak R, Irimia M, Kim T, Calarco JA, et al. Tissue-specific alternative splicing remodels protein-protein interaction networks. Mol Cell. 2012; 46(6):884–92. Epub 2012/07/04. https://doi.org/10.1016/j.molcel.2012.05.037 PMID: 22749401.

45. Yang X, Coulombe-Huntington J, Kang S, Sheynkman GM, Hao T, Richardson A, et al. Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing. Cell. 2016; 164(4):805–17. Epub 2016/02/13. https://doi.org/10.1016/j.cell.2016.01.029 PMID: 26871637; PubMed Central PMCID: PMC4882190.

46. Yang X, Lee WH, Sobott F, Papagrigoriou E, Robinson CV, Grossmann JG, et al. Structural basis for protein-protein interactions in the 14-3-3 protein family. Proc Natl Acad Sci U S A. 2006; 103 (46):17237–42. Epub 2006/11/07. https://doi.org/10.1073/pnas.0605779103 PMID: 17085597; PubMed Central PMCID: PMC1859916.

47. Dos Santos HG, Siltberg-Liberles J. Paralog-Specific Patterns of Structural Disorder and Phosphorylation in the Vertebrate SH3-SH2-Tyrosine Kinase Protein Family. Genome Biol Evol. 2016; 8(9):2806–25. Epub 2016/08/16. https://doi.org/10.1093/gbe/evw194 PMID: 27519537; PubMed Central PMCID: PMC5630953.

48. Shah NH, Amacher JF, Nocka LM, Kuriyan J. The Src module: an ancient scaffold in the evolution of cytoplasmic tyrosine kinases. Crit Rev Biochem Mol Biol. 2018; 53(5):535–63. Epub 2018/09/06. https://doi.org/10.1080/10409238.2018.1495173 PMID: 30183386; PubMed Central PMCID: PMC6328253.

49. Taylor SS, Keshwani MM, Steichen JM, Kornev AP. Evolution of the eukaryotic protein kinases as dynamic molecular switches. Philos T R Soc B. 2012; 367(1602)(1602):2517–28. PubMed PMID: CCC:000307802800003. https://doi.org/10.1098/rstb.2012.0054 PMID: 22889904

50. Tatarova Z, Brabek J, Rosel D, Novotny M. SH3 Domain Tyrosine Phosphorylation - Sites, Role and Evolution. Plos One. 2012; 7(5)(5):- PubMed PMID: CCC:000305336300013. https://doi.org/10.1371/journal.pone.0036310 PMID: 22615764

51. Ju T, Niu W, Guo J. Evolution of Src Homology 2 (SH2) Domain to Recognize Sulfotyrosine. ACS Chem Biol. 2016; 11(9):2551–7. Epub 2016/07/19. https://doi.org/10.1021/acschembio.6b00555 PMID: 27428792.

52. Liu BA, Engelmann BW, Nash PD. The language of SH2 domain interactions defines phosphotyrosine-mediated signal transduction. FEBS Lett. 2012; 586(17):2597–605. Epub 2012/05/10. https://doi.org/10.1016/j.febslet.2012.04.054 PMID: 22569091.

53. Dionne U, Bourgault E, Dube AK, Bradley D, Chartier FJM, Dandage R, et al. Protein context shapes the specificity of SH3 domain-mediated interactions in vivo. Nature Communications. 2021; 12(1)(1):- PubMed PMID: CCC:000630475400001. https://doi.org/10.1038/s41467-021-21873-2 PMID: 33712617

54. Goodacre N, Devkota P, Bae E, Wuchty S, Uetz P. Protein-protein interactions of human viruses. Semin Cell Dev Biol. 2020; 99:31–9. Epub 2018/07/22. https://doi.org/10.1016/j.semcdb.2018.07.018 PMID: 30031213; PubMed Central PMCID: PMC7102568.

55. Peleg O, Choi JM, Shakhnovich EI. Evolution of specificity in protein-protein interactions. Biophys J. 2014; 107(7):1686–96. Epub 2014/10/09. https://doi.org/10.1016/j.bpj.2014.08.004 PMID: 25296322; PubMed Central PMCID: PMC4190652.

56. Valdar WSJ, Thornton JM. Conservation helps to identify biologically relevant crystal contacts. Journal of Molecular Biology. 2001; 313(2)(2):399–416. PubMed PMID: CCC:000171836900012. https://doi.org/10.1006/jmbi.2001.5034 PMID: 11800565
39. Defoort J, Van de Peer Y, Carretero-Paulet L. The Evolution of Gene Duplicates in Angiosperms and the Impact of Protein-Protein Interactions and the Mechanism of Duplication. Genome Biol Evol. 2019; 11(8):2292–305. Epub 2019/08/01. https://doi.org/10.1093/gbe/evz156 PMID: 31364708; PubMed Central PMCID: PMC6735927.

40. Levy ED, Pereira-Leal JB. Evolution and dynamics of protein interactions and networks. Curr Opin Struct Biol. 2008; 18(3):349–57. Epub 2008/05/02. https://doi.org/10.1016/j.sbi.2008.03.003 PMID: 18448325.

41. Nooren IM, Thornton JM. Diversity of protein-protein interactions. EMBO J. 2003; 22(14):3486–92. Epub 2003/07/11. https://doi.org/10.1093/emboj/cdg359 PMID: 12853464; PubMed Central PMCID: PMC165629.

42. Nooren IMA, Thornton JM. Structural characterisation and functional significance of transient protein-protein interactions. Journal of Molecular Biology. 2003; 325(5)(5):991–1018. PubMed PMID: CCC:000180676300011. https://doi.org/10.1016/s0022-2836(02)01281-0 PMID: 12527304

43. Cafarelli TM, Desbuleux A, Wang Y, Choi SG, De Ridder D, Vidal M. Mapping, modeling, and characterization of protein-protein interactions on a proteomic scale. Curr Opin Struct Biol. 2017; 44:201–10. Epub 2017/06/03. https://doi.org/10.1016/j.sbi.2017.05.003 PMID: 28575754.

44. Ellis JD, Barrios-Rodiles M, Colak R, Irimia M, Kim T, Calarco JA, et al. Tissue-specific alternative splicing remodels protein-protein interaction networks. Mol Cell. 2012; 46(6):884–92. Epub 2012/07/04. https://doi.org/10.1016/j.molcel.2012.05.037 PMID: 22749401.

45. Yang X, Coulombe-Huntington J, Kang S, Sheynkman GM, Hao T, Richardson A, et al. Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing. Cell. 2016; 164(4):805–17. Epub 2016/02/13. https://doi.org/10.1016/j.cell.2016.01.029 PMID: 26871637; PubMed Central PMCID: PMC4882190.

46. Yang X, Lee WH, Sobott F, Papagrigoriou E, Robinson CV, Grossmann JG, et al. Structural basis for protein-protein interactions in the 14-3-3 protein family. Proc Natl Acad Sci U S A. 2006; 103 (46):17237–42. Epub 2006/11/07. https://doi.org/10.1073/pnas.0605779103 PMID: 17085597; PubMed Central PMCID: PMC1859916.

47. Dos Santos HG, Siltberg-Liberles J. Paralog-Specific Patterns of Structural Disorder and Phosphorylation in the Vertebrate SH3-SH2-Tyrosine Kinase Protein Family. Genome Biol Evol. 2016; 8(9):2806–25. Epub 2016/08/16. https://doi.org/10.1093/gbe/evw194 PMID: 27519537; PubMed Central PMCID: PMC5630953.

48. Shah NH, Amacher JF, Nocka LM, Kuriyan J. The Src module: an ancient scaffold in the evolution of cytoplasmic tyrosine kinases. Crit Rev Biochem Mol Biol. 2018; 53(5):535–63. Epub 2018/09/06. https://doi.org/10.1080/10409238.2018.1495173 PMID: 30183386; PubMed Central PMCID: PMC6328253.

49. Taylor SS, Keshwani MM, Steichen JM, Kornev AP. Evolution of the eukaryotic protein kinases as dynamic molecular switches. Philos T R Soc B. 2012; 367(1602)(1602):2517–28. PubMed PMID: CCC:000307802800003. https://doi.org/10.1098/rstb.2012.0054 PMID: 22889904

50. Tatarova Z, Brabek J, Rosel D, Novotny M. SH3 Domain Tyrosine Phosphorylation - Sites, Role and Evolution. Plos One. 2012; 7(5)(5):- PubMed PMID: CCC:000305336300013. https://doi.org/10.1371/journal.pone.0036310 PMID: 22615764

51. Ju T, Niu W, Guo J. Evolution of Src Homology 2 (SH2) Domain to Recognize Sulfotyrosine. ACS Chem Biol. 2016; 11(9):2551–7. Epub 2016/07/19. https://doi.org/10.1021/acschembio.6b00555 PMID: 27428792.

52. Liu BA, Engelmann BW, Nash PD. The language of SH2 domain interactions defines phosphotyrosine-mediated signal transduction. FEBS Lett. 2012; 586(17):2597–605. Epub 2012/05/10. https://doi.org/10.1016/j.febslet.2012.04.054 PMID: 22569091.

53. Dionne U, Bourgault E, Dube AK, Bradley D, Chartier FJM, Dandage R, et al. Protein context shapes the specificity of SH3 domain-mediated interactions in vivo. Nature Communications. 2021; 12(1)(1):- PubMed PMID: CCC:000630475400001. https://doi.org/10.1038/s41467-021-21873-2 PMID: 33712617

54. Goodacre N, Devkota P, Bae E, Wuchty S, Uetz P. Protein-protein interactions of human viruses. Semin Cell Dev Biol. 2020; 99:31–9. Epub 2018/07/22. https://doi.org/10.1016/j.semcdb.2018.07.018 PMID: 30031213; PubMed Central PMCID: PMC7102568.

55. Peleg O, Choi JM, Shakhnovich EI. Evolution of specificity in protein-protein interactions. Biophys J. 2014; 107(7):1686–96. Epub 2014/10/09. https://doi.org/10.1016/j.bpj.2014.08.004 PMID: 25296322; PubMed Central PMCID: PMC4190652.

56. Valdar WSJ, Thornton JM. Conservation helps to identify biologically relevant crystal contacts. Journal of Molecular Biology. 2001; 313(2)(2):399–416. PubMed PMID: CCC:000171836900012. https://doi.org/10.1006/jmbi.2001.5034 PMID: 11800565
57. Scott DE, Bayly AR, Abell C, Skidmore J. Small molecules, big targets: drug discovery faces the protein-protein interaction challenge. Nat Rev Drug Discov. 2016; 15(8):533–50. Epub 2016/04/07. https://doi.org/10.1038/nrd.2016.29 PMID: 27050677.

58. Thangudu RR, Bryant SH, Panchenko AR, Madej T. Modulating Protein-Protein Interactions with Small Molecules: The Importance of Binding Hotspots. Journal of Molecular Biology. 2012; 415(2)(2):443–53. PubMed PMID: CCC:000300032500017. https://doi.org/10.1016/j.jmb.2011.12.026 PMID: 22198293

59. Cooke MP, Perlmutter RM. Expression of a novel form of the fyn proto-oncogene in hematopoietic cells. New Biol. 1989; 1(1):66–74. Epub 1989/10/01. PubMed PMID: PMID: 2488273.

60. Picard C, Gilles A, Pontarotti P, Olive D, Collette Y. Cutting edge: recruitment of the ancestral fyn gene during emergence of the adaptive immune system. J Immunol. 2002; 168(6):2595–8. Epub 2002/03/09. https://doi.org/10.4049/jimmunol.168.6.2595 PMID: 11884421.

61. Goldsmith JF, Hall CG, Atkinson TP. Identification of an alternatively spliced isoform of the fyn tyrosine kinase. Biochem Biophys Res Commun. 2002; 298(4):501–4. Epub 2002/11/01. https://doi.org/10.1016/s0006-291x(02)02510-x PMID: 12408980.

62. Lee C, Low CY, Francis PT, Attems J, Wong PT, Lai MK, et al. An isoform-specific role of FynT tyrosine kinase in Alzheimer's disease. J Neurochem. 2016; 136(3):637–50. Epub 2015/11/13. https://doi.org/10.1111/jnc.13429 PMID: 26561212.

63. Marck C. 'DNA Strider': a 'C' program for the fast analysis of DNA and protein sequences on the Apple Macintosh family of computers. Nucleic Acids Res. 1988; 16(5):1829–36. https://doi.org/10.1093/nar/16.5.1829 PMID: 2832831; PubMed Central PMCID: PMC338177.

64. Cassonnet P, Rolloy C, Neveu G, Vidalain PO, Chantier T, Pellet J, et al. Benchmarking a luciferase complementation assay for detecting protein complexes. Nat Methods. 2011; 8(12):990–2. https://doi.org/10.1038/nmeth.1773 PMID: 22127214.

65. Le S, Josse J, Husson F. FactoMineR: An R package for multivariate analysis. J Stat Softw. 2008; 25(1):1–18. https://doi.org/10.18637/jss.v025.i01 PubMed PMID: WOS:000254619400001.

66. Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 2015; 12(1):7–8. https://doi.org/10.1038/nmeth.3213 PMID: 25549265; PubMed Central PMCID: PMC4428668.

67. Armon A, Graur D, Ben-Tal N. ConSurf: An algorithmic tool for the identification of functional regions in proteins by surface mapping of phylogenetic information. Journal of Molecular Biology. 2001; 307(1)(1):447–63. PubMed PMID: CCC:000167764900034. https://doi.org/10.1006/jmbi.2000.4474 PMID: 11243830

68. Ashkenazy H, Abadi S, Martz E, Chay O, Mayrose I, Pupko T, et al. ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res. 2016; 44(W1):W344–50. Epub 2016/05/12. https://doi.org/10.1093/nar/gkw408 PMID: 27166375; PubMed Central PMCID: PMC4987940.

69. Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004; 25(13):1605–12. Epub 2004/07/21. https://doi.org/10.1002/jcc.20084 PMID: 15264254.

70. Dinkel H, Michael S, Weatheritt RJ, Davey NE, Van Roey K, Altenberg B, et al. ELM—the database of eukaryotic linear motifs. Nucleic Acids Res. 2012; 40 part 2:242–51. https://doi.org/10.1093/nar/gkr1064 PMID: 22110040; PubMed Central PMCID: PMC3245074.

71. Hornbeck PV, Zhang B, Murray B, Kornhauser JM, Latham V, Skrzypek E. PhosphoSitePlus, 2014: mutations, PTMs and recalibrations. Nucleic Acids Res. 2015; 43 part 2:212–20. https://doi.org/10.1093/nar/gku1267 PMID: 25514926; PubMed Central PMCID: PMC4383998.

72. Kaneko T, Huang HM, Zhao B, Li L, Liu HD, Voss CK, et al. Loops Govern SH2 Domain Specificity by Controlling Access to Binding Pockets. Sci Signal. 2010; 3(120)(120):1–12. PubMed PMID: CCC:000277314600001. https://doi.org/10.1126/scisignal.2000796 PMID: 20442417

73. Liu BA, Jablonowski K, Shah EE, Engelmann BW, Jones RB, Nash PD. SH2 Domains Recognize Contextual Peptide Sequence Information to Determine Selectivity. Mol Cell Proteomics. 2010; 9(11)(11):2391–404. PubMed PMID: CCC:000285055900005. https://doi.org/10.1074/mcp.M110.001586 PMID: 20627867

74. Shimamura S, Sasaki K, Tanaka M. The Src substrate SKAP2 regulates actin assembly by interacting with WAVE2 and cortactin proteins. J Biol Chem. 2013; 288(2):1171–83. https://doi.org/10.1074/jbc.M112.386722 PMID: 23161539; PubMed Central PMCID: PMC3543001.
57. Scott DE, Bayly AR, Abell C, Skidmore J. Small molecules, big targets: drug discovery faces the protein-protein interaction challenge. Nat Rev Drug Discov. 2016; 15(8):533–50. Epub 2016/04/07. https://doi.org/10.1038/nrd.2016.29 PMID: 27050677.

58. Thangudu RR, Bryant SH, Panchenko AR, Madej T. Modulating Protein-Protein Interactions with Small Molecules: The Importance of Binding Hotspots. Journal of Molecular Biology. 2012; 415(2)(2):443–53. PubMed PMID: CCC:000300032500017. https://doi.org/10.1016/j.jmb.2011.12.026 PMID: 22198293

59. Cooke MP, Perlmutter RM. Expression of a novel form of the fyn proto-oncogene in hematopoietic cells. New Biol. 1989; 1(1):66–74. Epub 1989/10/01. PubMed PMID: PMID: 2488273.

60. Picard C, Gilles A, Pontarotti P, Olive D, Collette Y. Cutting edge: recruitment of the ancestral fyn gene during emergence of the adaptive immune system. J Immunol. 2002; 168(6):2595–8. Epub 2002/03/09. https://doi.org/10.4049/jimmunol.168.6.2595 PMID: 11884421.

61. Goldsmith JF, Hall CG, Atkinson TP. Identification of an alternatively spliced isoform of the fyn tyrosine kinase. Biochem Biophys Res Commun. 2002; 298(4):501–4. Epub 2002/11/01. https://doi.org/10.1016/s0006-291x(02)02510-x PMID: 12408980.

62. Lee C, Low CY, Francis PT, Attems J, Wong PT, Lai MK, et al. An isoform-specific role of FynT tyrosine kinase in Alzheimer's disease. J Neurochem. 2016; 136(3):637–50. Epub 2015/11/13. https://doi.org/10.1111/jnc.13429 PMID: 26561212.

63. Marck C. 'DNA Strider': a 'C' program for the fast analysis of DNA and protein sequences on the Apple Macintosh family of computers. Nucleic Acids Res. 1988; 16(5):1829–36. https://doi.org/10.1093/nar/16.5.1829 PMID: 2832831; PubMed Central PMCID: PMC338177.

64. Cassonnet P, Rolloy C, Neveu G, Vidalain PO, Chantier T, Pellet J, et al. Benchmarking a luciferase complementation assay for detecting protein complexes. Nat Methods. 2011; 8(12):990–2. https://doi.org/10.1038/nmeth.1773 PMID: 22127214.

65. Le S, Josse J, Husson F. FactoMineR: An R package for multivariate analysis. J Stat Softw. 2008; 25(1):1–18. https://doi.org/10.18637/jss.v025.i01 PubMed PMID: WOS:000254619400001.

66. Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 2015; 12(1):7–8. https://doi.org/10.1038/nmeth.3213 PMID: 25549265; PubMed Central PMCID: PMC4428668.

67. Armon A, Graur D, Ben-Tal N. ConSurf: An algorithmic tool for the identification of functional regions in proteins by surface mapping of phylogenetic information. Journal of Molecular Biology. 2001; 307(1)(1):447–63. PubMed PMID: CCC:000167764900034. https://doi.org/10.1006/jmbi.2000.4474 PMID: 11243830

68. Ashkenazy H, Abadi S, Martz E, Chay O, Mayrose I, Pupko T, et al. ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res. 2016; 44(W1):W344–50. Epub 2016/05/12. https://doi.org/10.1093/nar/gkw408 PMID: 27166375; PubMed Central PMCID: PMC4987940.

69. Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004; 25(13):1605–12. Epub 2004/07/21. https://doi.org/10.1002/jcc.20084 PMID: 15264254.

70. Dinkel H, Michael S, Weatheritt RJ, Davey NE, Van Roey K, Altenberg B, et al. ELM—the database of eukaryotic linear motifs. Nucleic Acids Res. 2012; 40 part 2:242–51. https://doi.org/10.1093/nar/gkr1064 PMID: 22110040; PubMed Central PMCID: PMC3245074.

71. Hornbeck PV, Zhang B, Murray B, Kornhauser JM, Latham V, Skrzypek E. PhosphoSitePlus, 2014: mutations, PTMs and recalibrations. Nucleic Acids Res. 2015; 43 part 2:212–20. https://doi.org/10.1093/nar/gku1267 PMID: 25514926; PubMed Central PMCID: PMC4383998.

72. Kaneko T, Huang HM, Zhao B, Li L, Liu HD, Voss CK, et al. Loops Govern SH2 Domain Specificity by Controlling Access to Binding Pockets. Sci Signal. 2010; 3(120)(120):1–12. PubMed PMID: CCC:000277314600001. https://doi.org/10.1126/scisignal.2000796 PMID: 20442417

73. Liu BA, Jablonowski K, Shah EE, Engelmann BW, Jones RB, Nash PD. SH2 Domains Recognize Contextual Peptide Sequence Information to Determine Selectivity. Mol Cell Proteomics. 2010; 9(11)(11):2391–404. PubMed PMID: CCC:000285055900005. https://doi.org/10.1074/mcp.M110.001586 PMID: 20627867

74. Shimamura S, Sasaki K, Tanaka M. The Src substrate SKAP2 regulates actin assembly by interacting with WAVE2 and cortactin proteins. J Biol Chem. 2013; 288(2):1171–83. https://doi.org/10.1074/jbc.M112.386722 PMID: 23161539; PubMed Central PMCID: PMC3543001.