Methionine oxidation selectively enhances T cell reactivity against a melanoma antigen.

The impact of the peptide amino acids side-chain modifications on the immunological recognition has been scarcely explored. We investigate here the effect of methionine oxidation on the antigenicity of the melanoma immunodominant peptide 369-YMDGTMSQV-377 (YMD). Using CD8+ T cell activation assays, we found that the antigenicity of the sulfoxide form is higher when compared to the YMD peptide. This is consistent with free energy computations performed on HLA-A∗02:01/YMD/TCR complex showing that this is lowered upon oxidation, paired with a steep increase in order at atomic level. Oxidized YMD forms were identified at the melanoma cell surface by LC-MS/MS analysis. These results demonstrate that methionine oxidation in the antigenic peptides may generate altered peptide ligands with increased antigenicity, and that this oxidation may occur in vivo, opening up the possibility that high-affinity CD8+ T cells might be naturally primed in the course of melanoma progression, as a result of immunosurveillance.

NLRscape: an atlas of plant NLR proteins.

NLRscape is a webserver that curates a collection of over 80 000 plant protein sequences identified in UniProtKB to contain NOD-like receptor signatures, and hosts in addition a number of tools aimed at the exploration of the complex sequence landscape of this class of plant proteins. Each entry gathers sequence information, domain and motif annotations from multiple third-party sources but also in-house advanced annotations aimed at addressing caveats of the existing broad-based annotations. NLRscape provides a top-down perspective of the NLR sequence landscape but also services for assisting a bottom-up approach starting from a given input sequence. Sequences are clustered by their domain organization layout, global homology and taxonomic spread-in order to allow analysis of how particular traits of an NLR family are scattered within the plant kingdom. Tools are provided for users to locate their own protein of interest in the overall NLR landscape, generate custom clusters centered around it and perform a large number of sequence and structural analyses using included interactive online instruments. Amongst these, we mention: taxonomy distribution plots, homology cluster graphs, identity matrices and interactive MSA synchronizing secondary structure and motif predictions. NLRscape can be found at: https://nlrscape.biochim.ro/.

NLRexpress-A bundle of machine learning motif predictors-Reveals motif stability underlying plant Nod-like receptors diversity.

Examination of a collection of over 80,000 Plant Nod-like receptors (NLRs) revealed an overwhelming sequence diversity underlying functional specificity of pathogen detection, signaling and cooperativity. The NLR canonical building blocks-CC/TIR/RPW8, NBS and LRR-contain, however, a number of conserved sequence motifs showing a significant degree of invariance amongst different NLR groups. To identify these motifs we developed NLRexpress-a bundle of 17 machine learning (ML)-based predictors, able to swiftly and precisely detect CC, TIR, NBS, and LRR motifs while minimizing computing time without accuracy losses-aimed as an instrument scalable for screening overall proteomes, transcriptomes or genomes for identifying integral NLRs and discriminating them against incomplete sequences lacking key motifs. These predictors were further used to screen a subset of ∼34,000 regular plant NLR sequences. Motifs were analyzed using unsupervised ML techniques to assess the structural correlations hidden underneath pattern variabilities. Both the NB-ARC switch domain which admittedly is the most conserved region of NLRs and the highly diverse LRR domain with its vastly variable lengths and repeat irregularities-show well-defined relations between motif subclasses, highlighting the importance of structural invariance in shaping NLR sequence diversity. The online NLRexpress webserver can be accessed at https://nlrexpress.biochim.ro.

Effective Use of Empirical Data for Virtual Screening against APJR GPCR Receptor.

Alzheimer's disease is a neurodegenerative disorder incompatible with normal daily activity, affecting one in nine people. One of its potential targets is the apelin receptor (APJR), a G-protein coupled receptor, which presents considerably high expression levels in the central nervous system. In silico studies of APJR drug-like molecule binding are in small numbers while high throughput screenings (HTS) are already sufficiently many to devise efficient drug design strategies. This presents itself as an opportunity to optimize different steps in future large scale virtual screening endeavours. Here, we ran a first stage docking simulation against a library of 95 known binders and 3829 generated decoys in an effort to improve the rescoring stage. We then analyzed receptor binding site structure and ligands binding poses to describe their interactions. As a result, we devised a simple and straightforward virtual screening Stage II filtering score based on search space extension followed by a geometric estimation of the ligand-binding site fitness. Having this score, we used an ensemble of receptors generated by Hamiltonian Monte Carlo simulation and reported the results. The improvements shown herein prove that our ensemble docking protocol is suited for APJR and can be easily extrapolated to other GPCRs.

Robosample: A rigid-body molecular simulation program based on robot mechanics.

BACKGROUND: Compared with all-atom molecular dynamics (MD), constrained MD methods allow for larger time steps, potentially reducing computational cost. For this reason, there has been continued interest in improving constrained MD algorithms to increase configuration space sampling in molecular simulations. METHODS: Here, we introduce Robosample, a software package that implements high-performance constrained dynamics algorithms, originally developed for robotics, and applies them to simulations of biomolecular systems. As in the gMolmodel package developed by Spiridon and Minh in 2017, Robosample uses Constrained Dynamics Hamiltonian Monte Carlo (CDHMC) as a Gibbs sampling move - a type of Monte Carlo move where a subset of coordinates is allowed to change. In addition to the previously described Cartesian and torsional dynamics moves, Robosample implements spherical and cylindrical joints that can be distributed along the molecule by the user. RESULTS: In alanine dipeptide simulations, the free energy surface is recovered by mixing fully flexible with torsional, cylindrical, or spherical dynamics moves. Ramachandran dynamics, where only the two key torsions are mobile, accelerate the slowest transition by an order of magnitude. We also show that simulations of a complex glycan cover significantly larger regions of the configuration space when mixed with constrained dynamics. MAJOR CONCLUSIONS: Robosample is a tool of choice for efficient conformational sampling of large biomolecules. GENERAL SIGNIFICANCE: Robosample is intended as a reliable and user-friendly simulation package for fast biomolecular sampling that does not require extensive expertise in mechanical engineering or in the statistical mechanics of reduced coordinates.

LRRpredictor-A New LRR Motif Detection Method for Irregular Motifs of Plant NLR Proteins Using an Ensemble of Classifiers.

Leucine-rich-repeats (LRRs) belong to an archaic procaryal protein architecture that is widely involved in protein-protein interactions. In eukaryotes, LRR domains developed into key recognition modules in many innate immune receptor classes. Due to the high sequence variability imposed by recognition specificity, precise repeat delineation is often difficult especially in plant NOD-like Receptors (NLRs) notorious for showing far larger irregularities. To address this problem, we introduce here LRRpredictor, a method based on an ensemble of estimators designed to better identify LRR motifs in general but particularly adapted for handling more irregular LRR environments, thus allowing to compensate for the scarcity of structural data on NLR proteins. The extrapolation capacity tested on a set of annotated LRR domains from six immune receptor classes shows the ability of LRRpredictor to recover all previously defined specific motif consensuses and to extend the LRR motif coverage over annotated LRR domains. This analysis confirms the increased variability of LRR motifs in plant and vertebrate NLRs when compared to extracellular receptors, consistent with previous studies. Hence, LRRpredictor is able to provide novel insights into the diversification of LRR domains and a robust support for structure-informed analyses of LRRs in immune receptor functioning.

Genome-wide functional analyses of plant coiled-coil NLR-type pathogen receptors reveal essential roles of their N-terminal domain in oligomerization, networking, and immunity.

The ability to induce a defense response after pathogen attack is a critical feature of the immune system of any organism. Nucleotide-binding leucine-rich repeat receptors (NLRs) are key players in this process and perceive the occurrence of nonself-activities or foreign molecules. In plants, coevolution with a variety of pests and pathogens has resulted in repertoires of several hundred diverse NLRs in single individuals and many more in populations as a whole. However, the mechanism by which defense signaling is triggered by these NLRs in plants is poorly understood. Here, we show that upon pathogen perception, NLRs use their N-terminal domains to transactivate other receptors. Their N-terminal domains homo- and heterodimerize, suggesting that plant NLRs oligomerize upon activation, similar to the vertebrate NLRs; however, consistent with their large number in plants, the complexes are highly heterometric. Also, in contrast to metazoan NLRs, the N-terminus, rather than their centrally located nucleotide-binding (NB) domain, can mediate initial partner selection. The highly redundant network of NLR interactions in plants is proposed to provide resilience to perturbation by pathogens.

Distinct Roles of Non-Overlapping Surface Regions of the Coiled-Coil Domain in the Potato Immune Receptor Rx1.

The intracellular immune receptor Rx1 of potato (Solanum tuberosum), which confers effector-triggered immunity to Potato virus X, consists of a central nucleotide-binding domain (NB-ARC) flanked by a carboxyl-terminal leucine-rich repeat (LRR) domain and an amino-terminal coiled-coil (CC) domain. Rx1 activity is strictly regulated by interdomain interactions between the NB-ARC and LRR, but the contribution of the CC domain in regulating Rx1 activity or immune signaling is not fully understood. Therefore, we used a structure-informed approach to investigate the role of the CC domain in Rx1 functionality. Targeted mutagenesis of CC surface residues revealed separate regions required for the intramolecular and intermolecular interaction of the CC with the NB-ARC-LRR and the cofactor Ran GTPase-activating protein2 (RanGAP2), respectively. None of the mutant Rx1 proteins was constitutively active, indicating that the CC does not contribute to the autoinhibition of Rx1 activity. Instead, the CC domain acted as a modulator of downstream responses involved in effector-triggered immunity. Systematic disruption of the hydrophobic interface between the four helices of the CC enabled the uncoupling of cell death and disease resistance responses. Moreover, a strong dominant negative effect on Rx1-mediated resistance and cell death was observed upon coexpression of the CC alone with full-length Rx1 protein, which depended on the RanGAP2-binding surface of the CC. Surprisingly, coexpression of the N-terminal half of the CC enhanced Rx1-mediated resistance, which further indicated that the CC functions as a scaffold for downstream components involved in the modulation of disease resistance or cell death signaling.

Hamiltonian Monte Carlo with Constrained Molecular Dynamics as Gibbs Sampling.

Compared to fully flexible molecular dynamics, simulations of constrained systems can use larger time steps and focus kinetic energy on soft degrees of freedom. Achieving ergodic sampling from the Boltzmann distribution, however, has proven challenging. Using recent generalizations of the equipartition principle and Fixman potential, here we implement Hamiltonian Monte Carlo based on constrained molecular dynamics as a Gibbs sampling move. By mixing Hamiltonian Monte Carlo based on fully flexible and torsional dynamics, we are able to reproduce free energy landscapes of simple model systems and enhance sampling of macrocycles.

Random mutagenesis of the nucleotide-binding domain of NRC1 (NB-LRR Required for Hypersensitive Response-Associated Cell Death-1), a downstream signalling nucleotide-binding, leucine-rich repeat (NB-LRR) protein, identifies gain-of-function mutations in the nucleotide-binding pocket.

Plant nucleotide-binding, leucine-rich repeat (NB-LRR) proteins confer immunity to pathogens possessing the corresponding avirulence proteins. Activation of NB-LRR proteins is often associated with induction of the hypersensitive response (HR), a form of programmed cell death. NRC1 (NB-LRR Required for HR-Associated Cell Death-1) is a tomato (Solanum lycopersicum) NB-LRR protein that participates in the signalling cascade leading to resistance to the pathogens Cladosporium fulvum and Verticillium dahliae. To identify mutations in NRC1 that cause increased signalling activity, we generated a random library of NRC1 variants mutated in their nucleotide-binding domain and screened them for the ability to induce an elicitor-independent HR in Nicotiana tabacum. Screening of 1920 clones retrieved 11 gain-of-function mutants, with 10 of them caused by a single amino acid substitution. All substitutions are located in or very close to highly conserved motifs within the nucleotide-binding domain, suggesting modulation of the signalling activity of NRC1. Three-dimensional modelling of the nucleotide-binding domain of NRC1 revealed that the targeted residues are centred around the bound nucleotide. Our mutational approach has generated a wide set of novel gain-of-function mutations in NRC1 and provides insight into how the activity of this NB-LRR is regulated.

Three-dimensional modeling and diversity analysis reveals distinct AVR recognition sites and evolutionary pathways in wild and domesticated wheat Pm3 R genes.

The Pm3 gene confers resistance against wheat powdery mildew. Studies of Pm3 diversity have shown that Pm3 alleles isolated from southern populations of wild emmer wheat located in Lebanon, Jordan, Israel, and Syria are more diverse and more distant from bread wheat alleles than alleles from the northern wild wheat populations located in Turkey, Iran, and Iraq. Therefore, southern populations from Israel were studied extensively to reveal novel Pm3 alleles that are absent from the cultivated gene pool. Candidate Pm3 genes were isolated via a polymerase chain reaction cloning approach. Known and newly identified Pm3 genes were subjected to variation analysis and polymorphic amino acid residues were superimposed on a three-dimensional (3D) model of PM3. The region of highest interspecies diversity between Triticum aestivum and T. dicoccoides lies in leucine-rich repeats (LRR) 19 to 24, whereas most intraspecies diversity in T. aestivum is located in LRR 25 to 28. Interestingly, these two regions are separated by one large LRR whose propensity for flexibility facilitates the conformation of the PM3 LRR domain into two differently structured models. The combination of evolutionary and protein 3D structure analysis revealed that Pm3 genes in wild and domesticated wheat show different evolutionary histories which might have been triggered through different interactions with the powdery mildew pathogen.

Structural determinants at the interface of the ARC2 and leucine-rich repeat domains control the activation of the plant immune receptors Rx1 and Gpa2.

Many plant and animal immune receptors have a modular nucleotide-binding-leucine-rich repeat (NB-LRR) architecture in which a nucleotide-binding switch domain, NB-ARC, is tethered to a LRR sensor domain. The cooperation between the switch and sensor domains, which regulates the activation of these proteins, is poorly understood. Here, we report structural determinants governing the interaction between the NB-ARC and LRR in the highly homologous plant immune receptors Gpa2 and Rx1, which recognize the potato cyst nematode Globodera pallida and Potato virus X, respectively. Systematic shuffling of polymorphic sites between Gpa2 and Rx1 showed that a minimal region in the ARC2 and N-terminal repeats of the LRR domain coordinate the activation state of the protein. We identified two closely spaced amino acid residues in this region of the ARC2 (positions 401 and 403) that distinguish between autoactivation and effector-triggered activation. Furthermore, a highly acidic loop region in the ARC2 domain and basic patches in the N-terminal end of the LRR domain were demonstrated to be required for the physical interaction between the ARC2 and LRR. The NB-ARC and LRR domains dissociate upon effector-dependent activation, and the complementary-charged regions are predicted to mediate a fast reassociation, enabling multiple rounds of activation. Finally, we present a mechanistic model showing how the ARC2, NB, and N-terminal half of the LRR form a clamp, which regulates the dissociation and reassociation of the switch and sensor domains in NB-LRR proteins.

RAG and HMGB1 create a large bend in the 23RSS in the V(D)J recombination synaptic complexes.

During V(D)J recombination, recombination activating gene proteins RAG1 and RAG2 generate DNA double strand breaks within a paired complex (PC) containing two complementary recombination signal sequences (RSSs), the 12RSS and 23RSS, which differ in the length of the spacer separating heptamer and nonamer elements. Despite the central role of the PC in V(D)J recombination, little is understood about its structure. Here, we use fluorescence resonance energy transfer to investigate the architecture of the 23RSS in the PC. Energy transfer was detected in 23RSS substrates in which the donor and acceptor fluorophores flanked the entire RSS, and was optimal under conditions that yield a cleavage-competent PC. The data are most easily explained by a dramatic bend in the 23RSS that reduces the distance between these flanking regions from >160 Å in the linear substrate to <80 Å in the PC. Analysis of multiple fluorescent substrates together with molecular dynamics modeling yielded a model in which the 23RSS adopts a U shape in the PC, with the spacer located centrally within the bend. We propose that this large bend facilitates simultaneous recognition of the heptamer and nonamer, is critical for proper positioning of the active site and contributes to the 12/23 rule.

Tyrosinase degradation is prevented when EDEM1 lacks the intrinsically disordered region.

EDEM1 is a mannosidase-like protein that recruits misfolded glycoproteins from the calnexin/calreticulin folding cycle to downstream endoplasmic reticulum associated degradation (ERAD) pathway. Here, we investigate the role of EDEM1 in the processing of tyrosinase, a tumour antigen overexpressed in melanoma cells. First, we analyzed and modeled EDEM1 major domains. The homology model raised on the crystal structures of human and Saccharomyces cerevisiae ER class I α1,2-mannosidases reveals that the major mannosidase domain located between aminoacids 121-598 fits with high accuracy. We have further identified an N-terminal region located between aminoacids 40-119, predicted to be intrinsically disordered (ID) and susceptible to adopt multiple conformations, hence facilitating protein-protein interactions. To investigate these two domains we have constructed an EDEM1 deletion mutant lacking the ID region and a triple mutant disrupting the glycan-binding domain and analyzed their association with tyrosinase. Tyrosinase is a glycoprotein partly degraded endogenously by ERAD and the ubiquitin proteasomal system. We found that the degradation of wild type and misfolded tyrosinase was enhanced when EDEM1 was overexpressed. Glycosylated and non-glycosylated mutants co-immunoprecipitated with EDEM1 even in the absence of its intact mannosidase-like domain, but not when the ID region was deleted. In contrast, calnexin and SEL 1L associated with the deletion mutant. Our data suggest that the ID region identified in the N-terminal end of EDEM1 is involved in the binding of glycosylated and non-glycosylated misfolded proteins. Accelerating tyrosinase degradation by EDEM1 overexpression may lead to an efficient antigen presentation and enhanced elimination of melanoma cells.

Ancient diversity of splicing motifs and protein surfaces in the wild emmer wheat (Triticum dicoccoides) LR10 coiled coil (CC) and leucine-rich repeat (LRR) domains.

In this study, we explore the diversity and its distribution along the wheat leaf rust resistance protein LR10 three-dimensional structure. Lr10 is a leaf rust resistance gene encoding a coiled coil-nucleotide-binding site-leucine-rich repeat (CC-NBS-LRR) class of protein. Lr10 was cloned and sequenced from 58 accessions representing diverse habitats of wild emmer wheat in Israel. Nucleotide diversity was very high relative to other wild emmer wheat genes (π= 0.029). The CC domain was found to be the most diverse domain and subject to positive selection. Superimposition of the diversity on the CC three-dimensional structure showed that some of the variable and positively selected residues were solvent exposed and may interact with other proteins. The LRR domain was relatively conserved, but showed a hotspot of amino acid variation between two haplotypes in the ninth repeat. This repeat was longer than the other LRRs, and three-dimensional modelling suggested that an extensive α helix structure was formed in this region. The two haplotypes also differed in splicing regulation motifs. In genotypes with one haplotype, an intron was alternatively spliced in this region, whereas, in genotypes with the other haplotype, this intron did not splice at all. The two haplotypes are proposed to be ancient and maintained by balancing selection.

C-terminus glycans with critical functional role in the maturation of secretory glycoproteins.

The N-glycans of membrane glycoproteins are mainly exposed to the extracellular space. Human tyrosinase is a transmembrane glycoprotein with six or seven bulky N-glycans exposed towards the lumen of subcellular organelles. The central active site region of human tyrosinase is modeled here within less than 2.5 Å accuracy starting from Streptomyces castaneoglobisporus tyrosinase. The model accounts for the last five C-terminus glycosylation sites of which four are occupied and indicates that these cluster in two pairs--one in close vicinity to the active site and the other on the opposite side. We have analyzed and compared the roles of all tyrosinase N-glycans during tyrosinase processing with a special focus on the proximal to the active site N-glycans, s6:N337 and s7:N371, versus s3:N161 and s4:N230 which decorate the opposite side of the domain. To this end, we have constructed mutants of human tyrosinase in which its seven N-glycosylation sites were deleted. Ablation of the s6:N337 and s7:N371 sites arrests the post-translational productive folding process resulting in terminally misfolded mutants subjected to degradation through the mannosidase driven ERAD pathway. In contrast, single mutants of the other five N-glycans located either opposite to the active site or into the N-terminus Cys1 extension of tyrosinase are temperature-sensitive mutants and recover enzymatic activity at the permissive temperature of 31°C. Sites s3 and s4 display selective calreticulin binding properties. The C-terminus sites s7 and s6 are critical for the endoplasmic reticulum retention and intracellular disposal. Results herein suggest that individual N-glycan location is critical for the stability, regional folding control and secretion of human tyrosinase and explains some tyrosinase gene missense mutations associated with oculocutaneous albinism type I.

Coiled-coil domain-dependent homodimerization of intracellular barley immune receptors defines a minimal functional module for triggering cell death.

Plants and animals have evolved structurally related innate immune sensors, designated NLRs, to detect intracellular nonself molecules. NLRs are modular, consisting of N-terminal coiled-coil (CC) or TOLL/interleukin-1 receptor (TIR) domains, a central nucleotide-binding (NB) domain, and C-terminal leucine-rich repeats (LRRs). The polymorphic barley mildew A (MLA) locus encodes CC-containing allelic immune receptors recognizing effectors of the pathogenic powdery mildew fungus. We report the crystal structure of an MLA receptor's invariant CC domain, which reveals a rod-shaped homodimer. MLA receptors also self-associate in vivo, but self-association appears to be independent of effector-triggered receptor activation. MLA CC mutants that fail to self-interact impair in planta cell death activity triggered by the CC domain alone and by an autoactive full-length MLA receptor that mimics its ATP-bound state. Thus, CC domain-dependent dimerization of the immune sensor defines a minimal functional unit and implies a role for the dimeric CC module in downstream immune signaling.

Nucleocytoplasmic distribution is required for activation of resistance by the potato NB-LRR receptor Rx1 and is balanced by its functional domains.

The Rx1 protein, as many resistance proteins of the nucleotide binding-leucine-rich repeat (NB-LRR) class, is predicted to be cytoplasmic because it lacks discernable nuclear targeting signals. Here, we demonstrate that Rx1, which confers extreme resistance to Potato virus X, is located both in the nucleus and cytoplasm. Manipulating the nucleocytoplasmic distribution of Rx1 or its elicitor revealed that Rx1 is activated in the cytoplasm and cannot be activated in the nucleus. The coiled coil (CC) domain was found to be required for accumulation of Rx1 in the nucleus, whereas the LRR domain promoted the localization in the cytoplasm. Analyses of structural subdomains of the CC domain revealed no autonomous signals responsible for active nuclear import. Fluorescence recovery after photobleaching and nuclear fractionation indicated that the CC domain binds transiently to large complexes in the nucleus. Disruption of the Rx1 resistance function and protein conformation by mutating the ATP binding phosphate binding loop in the NB domain, or by silencing the cochaperone SGT1, impaired the accumulation of Rx1 protein in the nucleus, while Rx1 versions lacking the LRR domain were not affected in this respect. Our results support a model in which interdomain interactions and folding states determine the nucleocytoplasmic distribution of Rx1.

Interface analysis of the complex between ERK2 and PTP-SL.

The activity of ERK2, an essential component of MAP-kinase pathway, is under the strict control of various effector proteins. Despite numerous efforts, no crystal structure of ERK2 complexed with such partners has been obtained so far. PTP-SL is a major regulator of ERK2 activity. To investigate the ERK2-PTP-SL complex we used a combined method based on cross-linking, MALDI-TOF analysis, isothermal titration calorimetry, molecular modeling and docking. Hence, new insights into the stoichiometry, thermodynamics and interacting regions of the complex are obtained and a structural model of ERK2-PTP-SL complex in a state consistent with PTP-SL phosphatase activity is developed incorporating all the experimental constraints available at hand to date. According to this model, part of the N-terminal region of PTP-SL has propensity for intrinsic disorder and becomes structured within the complex with ERK2. The proposed model accounts for the structural basis of several experimental findings such as the complex-dissociating effect of ATP, or PTP-SL blocking effect on the ERK2 export to the nucleus. A general observation emerging from this model is that regions involved in substrate binding in PTP-SL and ERK2, respectively are interacting within the interface of the complex.

A secreted SPRY domain-containing protein (SPRYSEC) from the plant-parasitic nematode Globodera rostochiensis interacts with a CC-NB-LRR protein from a susceptible tomato.

Esophageal gland secretions from nematodes are believed to include effectors that play important roles in plant parasitism. We have identified a novel gene family encoding secreted proteins specifically expressed in the dorsal esophageal gland of Globodera rostochiensis early in the parasitic cycle, and which contain the B30.2/SPRY domain. The secondary structure of these proteins, named the secreted SPRY domain-containing proteins (SPRYSEC), includes highly conserved regions folding into beta-strands interspersed with loops varying in sequence and in length. Mapping sequence diversity onto a three-dimensional structure model of the SPRYSEC indicated that most of the variability is in the extended loops that shape the so-called surface A in the SPRY domains. Seven of nine amino acid sites subjected to diversifying selection in the SPRYSEC are also at this surface. In both yeast-two-hybrid screening using a library from a susceptible tomato and in an in vitro pull-down assay, one of the SPRYSEC interacted with the leucine-rich repeat (LRR) region of a novel coiled-coil nucleotide-binding LRR protein, which is highly similar to members of the SW5 resistance gene cluster. Given that the tomato cultivar used is susceptible to nematode infection, this SPRYSEC could be an evolutionary intermediate that binds to a classical immune receptor but does not yet, or no longer, triggers a resistance response. Alternatively, this SPRYSEC may bind to the immune receptor to downregulate its activity.