Cupin: A candidate molecular structure for the Nep1-like protein family
© Cechin et al; licensee BioMed Central Ltd. 2008
Received: 20 November 2007
Accepted: 30 April 2008
Published: 30 April 2008
NEP1-like proteins (NLPs) are a novel family of microbial elicitors of plant necrosis. Some NLPs induce a hypersensitive-like response in dicot plants though the basis for this response remains unclear. In addition, the spatial structure and the role of these highly conserved proteins are not known.
We predict a 3d-structure for the β-rich section of the NLPs based on alignments, prediction tools and molecular dynamics. We calculated a consensus sequence from 42 NLPs proteins, predicted its secondary structure and obtained a high quality alignment of this structure and conserved residues with the two Cupin superfamily motifs. The conserved sequence GHRHDWE and several common residues, especially some conserved histidines, in NLPs match closely the two cupin motifs. Besides other common residues shared by dicot Auxin-Binding Proteins (ABPs) and NLPs, an additional conserved histidine found in all dicot ABPs was also found in all NLPs at the same position.
We propose that the necrosis inducing protein class belongs to the Cupin superfamily. Based on the 3d-structure, we are proposing some possible functions for the NLPs.
More than 10 years ago, a 24-kD necrosis and ethylene inducing protein, named NEP1, capable of triggering plant cell death was purified from culture filtrates of Fusarium oxysporum. Since then, several other NEP1-like proteins (NLPs) have been identified in diverse microorganisms; including bacteria, fungi, and oomycetes . In several cases, one species have more than one copy of NLPs and it is believed that several of these copies are pseudogenes [2, 3]. NLPs constitute a family of phytotoxic proteins that contains a secretory signal sequence and are able to elicit cell death and defense responses in a large number of dicot plants (reviewed by  and ). Most species with NLPs are plant pathogens but there are exceptions, since genes encoding NLPs have been detected in fungal and bacterial species that are not known to be pathogenic.
A recently published study identified three copies of NLPs in the basidiomycete Moniliophthora perniciosa (MpNEPs). M. perniciosa, the causal agent of the witches' broom disease in Theobroma cacao, is responsible for major crop losses in the Americas. The authors observed that despite the high sequence similarity, MpNEP1 and MpNEP2 present different structural features, and MpNEP2 activity was resistant to high temperatures. They also demonstrated that these genes are differentially expressed in two different life stages of the fungus .
All NLPs contain a conserved domain called necrosis-inducing Phytophthora protein 1 (NPP1) . The current lack of knowledge about functional domains, cellular targeting or protein binding motifs in this type of proteins complicates the unveiling of the actual function of NLPs . There is an increasing interest in the determination of their function, role in plant-pathogen interactions and molecular structure . The main conserved motif GHRHDWE shows no significant similarity to any currently known protein sequence and so provides no clues to NLPs function.
The Cupin superfamily was identified by Dunwell in 1998  and is among the most functionally diverse folding described to date, comprising both enzymatic and non-enzymatic members. These include helix-turn-helix transcription factor, AraC type transcription factor, oxalate decarboxylase, auxin-binding protein, globulins, etc. Many proteins on this superfamily have functions and chemical properties related to the NLPs: Auxin-Binding Proteins (ABPs) are hormone receptors and have a great influence on plant physiology; the related oxalate oxidase is involved in pathogen activities and germin-like proteins, apoplastic, glycoproteins are remarkably protease-resistant because of their cupin fold.
According to Dunwell et al. [8, 9] the cupin domain comprises two conserved motifs, each corresponding to two β-strands, separated by a less conserved region composed of another two β-strands with an intervening variable loop. The total size of the inter motif region varies from 11 residues to ca. 50 residues. The characteristic conserved sequence in motif 1 and 2 is g(x)5hxh(x)3,4e(x)6g and g(x)5pxg(x)2h(x)3n, respectively.
Homology searches using the NCBI-Blast produces no useful results in relation to the 3d-structure because the possible candidates have such a low score that they cannot be considered viable candidates. The result is a long list of necrosis and ethylene-inducing proteins, all of which are β-sheet-rich proteins, but none with useful information with associated 3d-structure. Any attempt to find other similar proteins based on their 1d-structure (sequence) to NLP protein results in other NLP proteins. In this article, we propose a 3d-structure for this protein family based on: (1) 1d-structure and conserved residues, (2) the supposed catalytic center, (3) the predicted signal sequences and target location, (4) cysteine and histidine conserved residues, and (5) the predicted 2d-structures. Our computational experiments in association with experimental clues point to the Cupin superfamily as the structure of the NLPs.
Results and Discussion
Gijzen and Nürnberger classify NLPs into two groups: those containing two cysteines (type I NLPs) and those containing four (type II NLPs). Type I NLPs occur in fungi, oomycetes and bacteria while type II NLPs do not occur in oomycetes . Sequence alignment differentiated the NLPs into these two main groups.
The statistical approach presented here parallels the different levels of phytopathogenicity shown by NLPs. They affect different species at different levels of intensity, being host specific.
The result of the alignment analysis of type I and II NLPs is shown in Figure 1. The first sequence represents all type I NLPs and will be called type I NLP consensus, the second sequence represents all type II NLPs and will be called type II NLP consensus. Because these consensus sequences statistically represent type I and II NLPs, we use them to obtain secondary structure predictions and to perform local alignments with proteins with known 3d-structure. Finally, we used the type I NLP consensus to build a 3d-structure. After obtaining type I and II consensus sequences, we submitted them to the PROF program in the PredictProtein site  to obtain a secondary (2d) structure prediction (see Figure 1).
We observed that the NLP 2d-structure may be divided into 5 parts or domains: (1) a signal peptide (positions 1–25) with an α-helix; (2) a start domain (positions 25–60) with 2–3 predicted β-strands and a predicted α-helix with low confidence level; (3) a coil flanked by two cysteines, C62 and C89, called c62c89-coil; (4) a β-strand-rich region (positions 90–220), composed by 9–10 β-strands and (5) an end domain (positions 220–270) with two predicted α-helices separated by a β-strand. The predicted central β-strand-rich region in NLPs included them in the all-β SCOP class of proteins.
In this article, we are proposing that the cupin fold is a suitable template for this region (residues 90–220) of the NLPs. The signal peptide is cut away from the sequence and the other regions play a secondary role in the main structure of the NLPs.
According to the literature, the difference between type I and II NLPs are cysteines C106 and C112, present only in type II NLPs . However, we observed other differences, both in the conservation pattern of residues and in the 2d-structure. For example, while histidine H29 and aspartate D30 are conserved among type I NLPs, they are underrepresented in type II NLPs. Also, the conserved sequence DxDxDgCY (positions 56–63) and the conserved histidines H179 and H185 in type II NLPs (not present in type I NLPs) is intriguing. Concerning the 2d-structure, the main difference is β-strand 6, not predicted in type I NLPs. In order to investigate which residues are essential in NLPs, in the sense that only them (and no other) could play an specific role (function and structure), and which of them may be substituted by other compatible ones (for example, same charge, hydropathicity, α-helix or β-strand bias, etc.), we count the number of common residues in each position of the alignment and draw them in a succession of histograms. For instance, in type I NLPs the conserved motif GHRHDWE is described by: [g89%a7%k4%] [h96%y4%] [r100%] [h100%] [d92%f4%y4%] [w100%] [e100%], where the letter represent the one letter code and the number is the frequency of the aa at this position. For type II NLPs the sequence of histograms is [g87%n13%] [h100%] [r74%k13%t13%] [h100%] [d100%] [w80%f13%l7%] [e100%]. We can see above that the first glycine may be substituted by alanine in type I NLPs. This means that a small flexible residue in this position fulfills (though glycine is more suitable) the necessary role for the structure and function of the NLPs.
We submitted all sequences in positions 132–138 and 132–139 in all the 42 NLPs to the search service of the Protein Data Bank (PDB) and obtained the following list of candidates ordered by e-value (see Methods): 1vj2 (e-value = 1.0), 1f51 (2.2), 1ixm (2.3), 2ftk (2.3), 1qtr (2.6), 1wm1 (2.6), 1x2b (2.6), 1x2e (2.6), 2c0h (2.7) and 2hi0 (3.9).
The 1vj2 structure, a protein with unknown function from Thermotoga maritima (a thermophilic Eubacteria with an optimum growth temperature of 80°C) belongs to the RmlC-like cupin SCOP superfamily, and the Mainly Beta CATH class. From the 2d-structure analysis, NLPs were recognized as β-sheet-rich structures , possibly belonging to the all-β SCOP class of proteins of which the RmlC-like cupin is a superfamily. 1vj2 presents a compatible number of β-strands with those of NLPs, their position relative to conserved residues is the same and finally its sequence rhshpwe is very similar to the pattern GHRHDWE of the NLPs. 1vj2 has four histidines acting as ligands for a manganese ion, what would explain the importance of this motif in the NLPs. The other candidates, 1f51, 1ixm and 2ftk present the sequence ghsrhdwm in the middle of an α-helix and for that they were discarded. Further, all the proteins 1qtr, 1wm1, 1x2b, 1x2e, 2c0h and 2hi0 posses many α-helices intermixed by few β-strands, and were discarded too.
Statistics for NLPs and cupins.
133 lst ligand
135 2nd ligand
139 3rd ligand
193 4th ligand
type I NLPs
type II NLPs
Comparing the residue histograms of NLPs and cupins in Table 1, we can see that the sequence hrhdxe is present in most NLPs and in most 3-residue cupins. Also, 75% of all cupin sequences and 95% of all NLPs have a histidine (fourth ligand) in the second motif and at position 193, respectively. These correspond to the most important residues in the general cupin pattern and, the substitution of any of these residues will reduce the ability of the protein to hold the metal ion, as is the case in some cupins. We concluded that the first motif in the cupins with its xhxhxxx [x-]e pattern corresponds to the GHRHDWE [gh] pattern of the NLPs and the histidine h75% in the second cupin motif corresponds to H19395% in the NLPs. Certainly this correspondence must be compatible with the 2d-structure, what we will see next.
In cupins, both positions 138 and 139 (see Table 1) typically contain negatively charged residues, such as aspartate (d) or glutamate (e). However, only E139 acts as a ligand for the metal ion. Therefore, although highly conserved in type I and II NLPs, e138 must be discarded as a viable ligand candidate. h139 could act as an ion ligand in type II NLPs, but only 19% of type I NLPs presents a histidine at this position. Among all 68 cupins in , only the sequences from Pyrococcus horikoshii and Arabidopsis thaliana have a histidine at this position, ihqhdweh (GenBank gi 3256432, a hypothetical protein) and ahhhtfgh (gi 1169199, DNA-damage-repair/toleration protein), respectively. The previously obtained T. maritima 1vj2 with its sequence rhshpweh (ligands are italicized) is included in this group, too. It seems that the third histidine confers an increased stability to the binding of the metal ion necessary in the extreme temperature living conditions of P. horikoshii and T. maritima.
The second most frequent residue at position 139, asparagine (n), is found only in two cupins among the 68 in : Arachis hypogaea (gi 1168390a) pkhadadn and Bacillus subtilis (gi 2636534) ahfdaytn. Because of the lack of histidines in positions 133 and 135, these cupins probably do not bind any metal ion.
Asparagine n139 is present in 26% of the NLPs and probably do not participate in the bind of any ion, too. Additionally, the fact that many NLPs (26%) have non-charged residues at position 139 raises the question if a charged residue is necessary at this position. Many cupins (29%) have uncharged residues (v8a7l3cg) at this position showing that these cupins do not need residue 139 at all as an ion ligand. For example, pirin 1j1l (dhphrgfet...hae) uses three histidines, (h133, h135, and h193) and a glutamate (e195) as ligands for Fe2+, but not e139, though it would be available to perform this function. Moreover, several cupins do not use any 3 rd ligand at this position. Examples are isopenicillin N synthase from Aspergillus nidulans, PDB code 1bk0 (sequence whedvslit...h and ion Fe3+); clavaminate synthase 1ds1 (sequence fhtemathr...h and ion Fe2+); hypothetical protein 1jr7 (sequence lhndgtyvee...h, and ion Fe2+) and anthocyanidin synthase from Arabidopsis thaliana 1gp4 (ahtdvsaltf...h, Fe3+). Even when present, e13953% is not very conserved in cupins for a residue that should bind to a metal ion. Additionally, site-directed mutagenesis e139 → q139 in the cupin acetylacetone dioxygenase Dke1 results in increased loss of the Fe2+ ion and reduced thermal stability , but its functional characteristics remain practically unchanged. We conclude that asparagine n139 does not act as a ligand for the metal ion, resulting in 61% of all NLPs with no ligand at this position. Finally, human cysteine dioxygenase 2ic1 (see Figure 2) has just 3 histidine ligands for the Fe2+ ion (h133, h135 and h193). The third histidine in this sequence ihdhtdshc...h does not act as a metal ligand and no other residue is necessary to hold the metal ion showing that NLPs could likewise, hold a metal ion at this site.
The Role of Cysteines
Fellbrich et al.  have shown that both cysteines C62 and C89 are conserved and necessary for the NLPs function. Also, the coil between them seems to encode a glycosylation site, which, for secreted proteins means protection against proteolysis, correct folding and thermal stability. Additionally, the highly conserved glycines G76G77 seem to promote a fold exactly in the middle of this coil enabling the cysteines to come together.
The analysis of the bonding pattern among cysteines resulted in a 90% confidence level for C62 and C89 to be forming a disulfide bridge in type I and II NLPs. A search for the pattern GnxsGGL in the PDB rendered the protein 1eh6, which has a turn at s75G76, supporting the hypothesis that both cysteines are disulfide bonded.
In relation to the other two cysteines present only in type II NLPs, the program DISULFIND attributes a probability of just 30% for the bonding of C106 to any other cysteine and 0% for C112. However, from the position of these two cysteines, it is not difficult to infer that they are bonded if β A and β B form a β-sheet (see Figure 3). These two cysteines seem to enforce that these β-strands should present this conformation. It is also possible that these two cysteines might be bonded to the two conserved histidines H133 and H134 by a zinc ion, such as in the zinc finger of WRKY proteins. WRKY-proteins have a special zinc-finger motif characterized by the pattern cx4,5cx22,23hxh, and type II NLPs have a similar pattern, cx4cx19hxh. WRKY-proteins are transcript factors with up to 100 representatives in A. thaliana . For instance, the protein AtWRKY6 is associated with both senescence- and defense-related processes . The structure of the WRKY proteins may be shared by type II NLPs. However, it is less probable that they share the same function.  suggests that NLP-induced necrosis requires interaction with a target site at the extracytoplasmic side of dicot plant plasma membrane. They show that the ectopic expression of NLP in dicot plants resulted in cell death only when the protein was delivered to the apoplast. However, Bae et al. have shown that NEP1 in the plant was localized at the cell wall and cytosol. This result indicates that NEP1 can penetrate through the plasma membrane but may not be able to penetrate organelles . It has been observed that NLPs are hydrophilic and not likely to cross the plasma membrane. Furthermore, our proposed model structure based on the ABP 1lr5 has many hydrophilic residues at the surface and the hydrophobic ones are buried supporting the hypothesis that NLPs are not able to cross the plasma membrane. Additionally, the rapid response of parsley protoplasts (approximately 150 seconds) to PpNPP1 (Phytophthora parasitica NPP1) is compatible with an interaction just at the plasma membrane level .
NLP 2d and 3d-Structure
The RmlC-like cupin superfamily belongs to the SCOP double-stranded beta-helix fold. Cupins are double-stranded because they are composed of two sequences of antiparallel strands linked with short turns. If the NLPs are cupins, then there should be a correspondence between the 2d-structures of cupins and those predicted for NLPs. Cupins are formed by 8–10 β-strands called [A]BCDE-FGHI [J].
Finally, this double strand turns like an helix building up a β-barrel of two β-sheets: CHEF and BIDG with their hydrophobic residues aiming at the interior of the barrel. Inside the barrel, in the hydrophobic pocket, we find the metal ion bound to its ligand, next to the top of the barrel (the bottom is closed by the E and F β-strands, see Figure 3). Cupins presenting a catalytic activity bind their substrates on the top of the barrel close to the metal ion at the hydrophobic pocket.
Analyzed NLP sequences.
Plant pathogens organism
GenBank accession number (NLP type)
Enwinia carotovora atroseptica
CAG75986••II (NipEca) [1 0]
AAC97382•• (NEP1)  
EDK02987II, EDJ98732II, EDJ96934 and EDJ94825 
AAS45247•• (His_VdNEP) 
CAJ98683•• (BeNEPl) and CAJ98684 (BeNEP2) 
XP 386193, XP 383570II, XP 387963II and XP 391669II 
ABQ53551•• (Mp NEP1) and ABO32369•• (Mp NEP2) 
AAY43363•• (NPP1), AAY43377° (NPP1.2) and AAY43378° (NPP1.3) 
AAX12401 and AAX12403 
AAM48170•• (PsojNIP), AAM48171 and AAM48172 
AAD53944•• (PaNie234) 
GenBank accession number (NLP type)
CAC40975•II (causes hemolysis) 
EAL86241 and EAL86501II
The Inter Motif Region (IMR)
The previous analysis about the 2d-structure of NLPs and cupins shows that predictions for the region delimited by β6 and β7 (putative E and F β-strands), which corresponds to the low conserved IMR in cupins, is a difficult task. It contains 22–32 (22–28 in type II NLP consensus) residues with 11 (11 in type II NLP consensus) residues in the coil of the IMR in type I NLPs and also has the conserved sequence S88%a50%H95%g74%. Therefore, we have chosen among the 68 cupins in  those that are similar in size. The most similar cupin to the NLPs' IMR is A. thaliana gi|461453, a possible a receptor for the hormone auxin. The 3d-structure of maize Auxin Binding Protein (ABP) has been already determined (1lrh and 1lr5 in PDB, see Figure 2). It has one β-barrel domain, is a dimer in solution, has 21 residues in the IMR, and 11 in coil. ABPs are involved in cell expansion and are located in the ER lumen, the plasma membrane, and the cell wall . It is ubiquitous amongst green plants .
It would be advantageous for the necrosis proteins to have control of the auxin-response in the host, for example changes in protoplast electrophysiology. Auxin induces H+ secretion into the cell wall causing hyperpolarization of the plasma membrane in Avena coleoptile cells , electrical response in tobacco protoplasts , and K+ currents in Nicotiana tabacum guard cells .
Auxin stimulates the growth of plant cells by regulating the activity of a H+-ATPase in the plasma membrane. Proton secretion by this transport enzyme acidifies the cell walls increasing their extensibility. The internal hydrostatic pressure of the cell then extends the walls. In the interaction between M. perniciosa and T. cacao, it has been suggested an auxin-inducing phase that causes malformation and an auxin-depleting secondary phase that kills the host . Kilaru et al. reported that the increase in auxin coincided with the phase transition of the fungus. This increase could be the effect of enhanced IAA (auxin) synthesis, suppression of IAA oxidase or secretion of IAA by the pathogen . The possibility that the 3d-structure of NLPs is similar to that of the cupins and, in particular to ABPs, places the NLPs in control of plant auxin receptors and of the resulting ionic currents. Additionally, Haberlach et al. have shown that the balance of cytokinin and auxin was an important factor in maintaining or eliminating resistance in plant tissues. More specifically, they show that P. parasitica resistant N. tabacum became susceptible under high cytokinin/auxin levels . It is conceivable that NLPs could compete with the natural ABPs resulting in an apparent increase of the cytokinin/auxin ratio in the plant. Other possibility is that NLPs could be auxin oxidases. Krupasager reported that in the dikaryotic stage, Marasmius perniciosa produces no significant amount of cytokinins or auxins but auxin-inactivating enzymes such as IAA-oxidase and laccase .
Differences between monocot and dicot ABPs and NLPs. Differences between monocot and dicot ABPs and NLPs (shown in boldface).
2d-Structure and glycosylation sites ([nx(st)])
Type I NLP consensus
β B β C
ηβ D β E
β B β C
ηβ D β E
β B β C
ηβ D β E
Type II NLP consensus
β A C
Cβ B β C
β A ,r
β A α
β B β C
β F nst
β A ,i
β A α
β B β C
ηβ D β E
β F nst
Residue histograms for NLPs and ABPs.
Type I NLPs
Type II NLPs
All 32 investigated dicot ABPs and 95% of all NLPs have a conserved histidine H162 at the EF-coil. The two NLP sequences which do not have it exactly at this position, but 4 positions downstream (PiNPP1.2 and PiNPP1.3), are inactive forms; probably originated by gene duplication of the most aggressive PiNPP1 from Phytophthora infestans. These two are expressed both in the biotrophic as in the necrotrophic phases, while PiNPP1 is expressed only in the latter phase .
Fellbrich  has shown that the last 8 residues in PiNPP1 may be deleted without loss of activity but not the last 20 residues (... ntdFGd ◇ AnvPmkdgnFlt ◇ kvgnayya). The sequence AnvPmkdgnFlt coincides 100% with the conserved residues in necrotic type I NLPs GxAnxP. It is interesting to observe that the C-terminal seems to be important for NLPs as well as for ABPs. The synthetic peptide from the C-terminal of Zea mays ABP wdedcfeaak, the 15-residue N. tabacum Nt-abp1 C-terminal peptide (ywdeecyqttswkdel) and the Nt-abp1 itself have been shown to induce hyperpolarization . However, contrary to the effect of auxins, NLPs cause depolarization, alkalization of the surrounding media and K+efflux .
The histograms of 32 dicot ABPs in the C-terminal resulted in yWDEqCyqtxxKDEL. Although conserved, mutagenesis experiments have shown that the sequence kdel is not important for the activity of ABPs and may be deleted and it is related to the two negative residues DE. Since NLPs do not have such a sequence, it is conceivable that NLPs compete with ABPs causing the previous discussed effects.
Besides ABPs, another candidate similar in terms of IMR size is Bacillus subtilis gi|2635598, a hypothetical protein with no determined structure similar to the human cysteine dioxygenase 2ic1. It is a monomer with an IMR size of 23 residues and a coil of 11 residues.  has aligned 10 cysteine dioxygenases of different organisms and the 100% conserved and functionally important residues are capitalized in Figure 2 for the sake of comparison with those in NLPs. We observe that some of these residues are among the most conserved in the necrotic type I NLPs: For instance, Y101 and R103 in β A , H133, H135 and H193 (as ion ligands) and D136.
Although an immunoassay study used for the detection of sugars in glycoconjugates did not reveal a carbohydrate moiety in PpNPP1 , glycosylation sites nxs are present in 64% of all necrotic type I NLPs (x = t6v) in the c62c89-coil. The glycosylation occurs at asparagine (n) residues in the so called nx(st) sequon and the efficiency of this process depends on the residue x.
Glycosylation is important in most cell-surface and secreted proteins and is often critical for the interaction with other subunits at the cell surface (recognition), protection against proteolytic attack, protein solubility and thermostability. For instance, P. infestans has evolved an arsenal of protease inhibitors to overcome the action of plant proteases .
The probability of glycosylation (see section Methods) of the sequon ntsg varies between 44 and 62%. For instance, PiNPP1.1 is 44%, PsojNIP 50%, PpNPP1 54% and BeNEP1 62%. Furthermore, MpNEP1 and MpNEP2 are predicted not to be glycosylated because of the bulky tryptophan (w) in their sequon (nws). His VdNEP and BeNEP2 have no sequon.
Additionally, an important difference between monocot and dicot ABPs  pointing out that ABPs and NLPs share the same 3d-structure is that dicot ABPs have a glycosylation site nis next to the beginning of the protein, that is not present in monocot ABPs (dis). This site is also present in several type I NLPs (see Table 3).
The 3d-structure of the NLPs remains to be determined experimentally. However, in this paper we presented several evidences indicating that they belong to the Cupin superfamily. Using a cupin template and the type I NLP consensus we were able to calculate a prediction for the 3d-structure of the β-strand rich portion (positions 90–220 in Figure 1) which presented stability for 3nS under molecular dynamics. This structure presents the classical signature of a cupin protein. Furthermore, the prediction of the structure of the upstream coil bordered by two cysteines remains to be addressed. Cysteines in the upstream coil of type I NLPs are disulfide bonded, simplifying the problem by removing this sequence from the analysis. However, the right positioning of this free coil becomes a new problem, which is more complex than the first one. Is it free to move around or does it make part of the β-sheets? Its glycosylation site points to some role in anchoring the protein on the cell surface.
Several 2d-structure predictions software agree with a central β-strand rich portion anked by α-helices making highly probable that the central part of the protein belongs to the all-β SCOP Class.
The conserved pattern of cysteines points to two different NLP types: type I NLPs containing two cysteines upstream of the β-barrel and type II NLPs, containing two additional cysteines in the left border of the β-barrel. Predictions show that the first two are bonded together, but not the other two. The similarity between the pattern of cysteines in the WRKY proteins and the type II NLPs is so high, that would not be surprising if the cysteines participate in metal biding in a zinc-finger conformation. Also, WRKY proteins participate in the pathogen detection system of the plant, an interesting "coincidence" for a protein involved in phytopathogenic activities. Unfortunately, we have only part of the WRKY protein structure, the zinc-finger part (C-terminal).
Figure 1 shows clearly the biological role played by the residues of the GHRHDWE motif holding a putatitive ion. For this purpose, a necessary 3 rd histidine is found downstream in the structure (H193) in the exact position according to the cupin structure and is also present in 95% of the NLPs analyzed. Only two NLPs do not have it, and they do not present necrotic activity. The relative position and number of histidines in the structure point to a metal ion containing protein. In addition, the exact metal ion remains to be determined but our study points to manganese or zinc.
The Cupin superfamily is very large with diverse functions. Many of these functions depend on a glutamate residue, which does not seem to be present in NLPs. When this residue is present, cupins are involved in enzymatic activities such as oxalate oxidase, decarboxylase, dioxygenases, etc. It remains to be determined experimentally if NLPs have some catalytic activity involving oxalate. In any case, the [Ca2+] cyt levels, ROS (H2O2) and oxalate are all intermixed in pathogen defense and sensing. Also, lignin processing is highly dependent on oxidases and peroxidases (cupins). It has been shown that germins function as oxalate oxidases (conversion of oxalate to CO2 and H2O2) and superoxide dismutase ( + 2H+ → H2O2+O2) . Any interference in such activities would be advantageous for the fungus because of the correlation between H2O2 and [Ca2+] cyt . Even with no catalytic activity, the fold is resistant to oxidation, a characteristic necessary for oxidases, decarboxylases and peroxidases.
A related class of cupins is the auxin binding proteins, which do not show catalytic activity but work as signal transducers in plant cells. They have many common structural features and conserved residues in relation to the NLPs. In this respect, the most remarkable feature is the conservation of histidine H162, present in 95% of the NLPs and in almost all dicot ABPs and is related to the fact that NLPs attack only dicot plants. Also, the way ABPs transmit the information to the cell is intimately related to ion channels. Correspondingly, the first plant reaction to the NLPs is an increase of ionic currents causing elevation of [Ca2+] cyt . The elevation of H2O2 (and other ROS species) is upstream and downstream of the [Ca2+] cyt elevation, and both [Ca2+] cyt elevation and H2O2 are known by their roles in senescence and necrotic activities. Besides all these molecular and functional clues, for example, the fungus M. perniciosa causes the formation of witches' broom on T. cacao . This reaction is typical of diverse plants in reaction to biotic stresses in its early phase and is related to the cytokinin/auxin balance.
There are several cupin candidates with a compatible IMR size, but only dicot ABPs seem to present some compatibility in relation to size and presence of the conserved histidine in the IMR of the NLPs and predicted glycosylation sites. Certainly, any auxin-like activity over the plant would be advantageous for the fungus (as discussed above). Last of all, as discussed by Fellbrich , NLPs are dependent on both, C and N-terminals of the protein for its activity, a feature shared by ABPs.
In order to obtain a non-redundant and representative set of NLPs, all sequences analyzed have a pairwise distance greater than 10%. The organisms and sequences include 42 NLPs (see Table 2).
These sequences are obtained computing the most frequent residues in the type I and II NLPs. If, in the alignment of the sequences, a gap is introduced in more than 50% of them, then the respective position is removed from the sequence. Therefore, type I and II NLP consensuses are the consensus of ≥ 50% of the sequences. Furthermore, residues present in more than 85% of all sequences are in boldface and capitalized (see Figure 1). For example, the sequence GHRHDWE occurs in ≥ 85% of all type I NLPs, indicating its important role for the protein function . Residues present in more than 70% of all sequences are in boldface (see Figure 1). For example, the same sequence GHrHDwE shows lower conservation in type II NLPs, as residues r and w occur in < 85% of type II NLP sequences.
We have performed an alignment (ClustalW with gapopen = 0) of type I NLPs for which we have experimental evidence of necrotic activity (denoted with one • and two filled circles ••) as shown in Table 2 and we denote the residues conserved in all necrotic sequences with an asterisk in Figure 1 (see alignment results in ).
The secondary structure predictions were performed using the PROF program in the PredictProtein site .
Searches in the PDB
The 10 best candidates were obtained by submitting all sequences in positions 132–138 and 132–139 in the 42 NLPs to the search service of the Protein Data Bank (PDB) site [48, 49] (Search Tool = Fasta). Each sequence generates a list of candidate proteins, which can be ordered by the e-value. The search for local sequences in the PDB was performed using Sequence features = motif.
Disulfide bonds and glycosylation analysis
The analysis of the bonding pattern among cysteines were performed with the DISULFIND program . We also used the DiANNA 1.1 disulfide bond prediction program  for the analysis of the BeNEP2. In order to determine the probability of glycosylation of the sequons in the c62c89-coil, we submitted all type I NLP sequences to the NetNGlyc 1.0 Server .
We used ClustalW with default parameters and gapopen = 0.0 to align the sequences shown in the Figure 2. Note that for the alignment, we used only the β-strand-rich region (β-barrel domain), positions 90–220. The 3d-structures were constructed with the molecular modelling program spdbv  using as template the monomer of the 1lr5 protein. First, we used Rasmol to cut the β-rich region of the 1lr5 protein. Then, we mutated the residues in the 1lr5 sequence to those residues in the type I NLP consensus, introduced the necessary residues where the alignment resulted in gaps for the 1lr5 sequence and deleted the residues in the cases the alignment produced gaps in the type I NLP consensus. After each step, we executed a minimization to fit the new residues to the old structure. We took the care to introduce and delete residues only in the coil regions of the protein. The resulting structure was minimized with GROMACS with a steepest descent algorithm until the energy of 0.01 kJ/mol was reached. Then, we performed a simulated annealing from 0 to 30 K during 100 pS in vacuum. The force field used was opls-aa/l. Finally, a simulated annealing of 100 pS from 0 to 10 K and 3 nS from 10 K to 273 K with molecular dynamics in explicit solvent (spc water model) was performed.
The authors thank Francisco Javier Medrano Martin for reviewing the manuscript. N. Lemke also would like to thank A. S. K. Braz for usefull discussions and suggestions. Work partially supported by CNPq (research grants 506414/2004-3 and 474278/2006-9) and FAPESP (research grant 2007/02827-9).
- Bailey BA: Purification of a protein from culture filtrates of Fusarium oxysporum that induces ethylene and necrosis in leaves of Erythroxylum coca. Phytopathology. 1995, 85: 1250-1255. 10.1094/Phyto-85-1250.View ArticleGoogle Scholar
- Gijzen M, Nürnberger T: Nep1-like proteins from plants pathogens: Recruitment and diversification of the NPP1 Domain Across Taxa. Phytochemistry. 2006, 16 (67): 1800-1807. 10.1016/j.phytochem.2005.12.008.View ArticleGoogle Scholar
- Tyler BM, Tripathy S, Zhang X, Dehal P, Jiang RH, Aerts A, Arredondo FD, Baxter L, Bensasson D, Beynon JL, Chapman J, Damasceno CMB, Dorrance AE, Dou D, Dickerman AW, Dubchak IL, Garbelotto M, Gijzen M, Gordon SC, Govers F, Grunwald NJ, Huang W, Ivors KL, Jones RW, Kamoun S, Krampis K, Lamour KH, Lee MK, McDonald WH, Medina M, Meijer HJG, Nordberg EK, Maclean DJ, Ospina-Giraldo MD, Morris PF, Phuntumart V, Putnam NH, Rash S, Rose JKC, Sakihama Y, Salamov AA, Savidor A, Scheuring CF, Smith BM, Sobral BWS, Terry A, Torto-Alalibo TA, Win J, Xu Z, Zhang H, Grigoriev IV, Rokhsar DS, Boore JL: Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis. Science. 2006, 313: 1261-1266. 10.1126/science.1128796.PubMedView ArticleGoogle Scholar
- Pemberton CL, Salmond GPC: The Nep1-like proteins – a growing family of microbial elicitors of plant necrosis. Mol Plant Pathol. 2004, 5 (4): 353-359. 10.1111/j.1364-3703.2004.00235.x.PubMedView ArticleGoogle Scholar
- Garcia O, Macedo J, Tibúrcio R, Zaparoli G, Rincones J, Bittencourt L, Ceita G, Micheli F, Gesteira A, Mariano A, Schiavinato M, Medrano F, Meinhardt L, Pereira G, Cascardo J: Characterization of necrosis and ethylene-inducing proteins (NEP) in the basidiomycete Moniliophthora perniciosa, the causal agent of witches' broom in Theobroma cacao. Mycol Res. 2007, 111: 443-455. 10.1016/j.mycres.2007.01.017.PubMedView ArticleGoogle Scholar
- Fellbrich G, Romanski A, Varet A, Blume B, Brunner F, Engelhardt S, Felix G, Kemmerling B, Krzymowska M, Nürnberger T: NPP1, a Phytophthora-associated trigger of plant defense in parsley and Arabidopsis. Plant J. 2002, 32 (3): 375-390. 10.1046/j.1365-313X.2002.01454.x.PubMedView ArticleGoogle Scholar
- Dunwell JM: Cupins: a new superfamily of functionally-diverse proteins that include germins and plant seed storage proteins. Biotechnol Genet Eng Rev. 1998, 15: 1-32.PubMedView ArticleGoogle Scholar
- Dunwell JM, Khuri S, Gane PJ: Microbial relatives of the seed storage proteins of higher plants: conservation of structure, and diversification of function during evolution of the cupin superfamily. Microbiol Mol Biol Rev. 2000, 64: 153-179. 10.1128/MMBR.64.1.153-179.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Dunwell JM, Culham A, Carter CE, Sosa-Aguirre CR, Goodenough PW: Evolution of functional diversity in the cupin superfamily. Trends Biochem Sci. 2001, 26: 740-745. 10.1016/S0968-0004(01)01981-8.PubMedView ArticleGoogle Scholar
- Pemberton CL, Whitehead NA, Sebaihia M, Bell KS, Hyman LJ, Harris SJ, Matlin AJ, Robson ND, Birch PRJ, Carr JP, Toth IK, Salmond GPC: Novel quorum-sensing-controlled genes in Erwinia carotovora subsp. carotovora: identification of a fungal elicitor homologue in a soft-rotting bacterium. Mol Plant Microbe Interact. 2005, 18: 343-353. 10.1094/MPMI-18-0343.PubMedView ArticleGoogle Scholar
- Jennings JC, Apel-Birkhold PC, Bailey BA, Anderson JD: Induction of ethylene biosynthesis and necrosis in weed leaves by a Fusarium oxysporum protein. Weed Sci. 2000, 48: 7-14. 10.1614/0043-1745(2000)048[0007:IOEBAN]2.0.CO;2.View ArticleGoogle Scholar
- Bailey BA, Apel-Birkhold PC, Luster DG: Expression of NEP1 by Fusarium oxysporum f. sp. erythroxyli after gene replacement and overexpression using polyethylene glycol-mediated transformation. Phytopathology. 2002, 92 (8): 833-841. 10.1094/PHYTO.2002.92.8.833.PubMedView ArticleGoogle Scholar
- Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu JR, Pan H, Read ND, Lee YH, Carbone I, Brown D, Oh YY, Donofrio N, Jeong JS, Soanes DM, Djonovic S, Kolomiets E, Rehmeyer C, Li W, Harding M, Kim S, Lebrun MH, Bohnert H, Coughlan S, Butler J, Calvo S, Ma LJ, Nicol R, Purcell S, Nusbaum C, Galagan JE, Birren BW: The genome sequence of the rice blast fungus Magnaporthe grisea . Nature. 2005, 434 (5): 980-986. 10.1038/nature03449.PubMedView ArticleGoogle Scholar
- Wang JY, Cai Y, Gou JY, Mao YB, Xu YH, Jiang WH, Chen XY: VdNEP, an Elicitor from Verticillium dahliae, induces cotton plant wilting. Appl Environ Microbiol. 2004, 70: 4989-4995. 10.1128/AEM.70.8.4989-4995.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Staats M, van Baarlen P, Schouten A, van Kan JAL, Bakker FT: Positive selection in phytotoxic protein-encoding genes of Botrytis species. Fungal Genet Biol. 2007, 44: 52-63. 10.1016/j.fgb.2006.07.003.PubMedView ArticleGoogle Scholar
- Han Y, Kim M, Lee S, Yun S, Lee Y: A novel F-box protein involved in sexual development and pathogenesis in Gibberella zeae. Mol Microbiol. 2007, 63: 768-779. 10.1111/j.1365-2958.2006.05557.x.PubMedView ArticleGoogle Scholar
- Kanneganti TD, Huitema E, Cakir C, Kamoun S: Synergistic interactions of the plant cell death pathways induced by Phytophthora infestans Nep1-like protein PiNPP1.1 and INF1 elicitin. Mol Plant Microbe Interact. 2006, 19: 854-863. 10.1094/MPMI-19-0854.PubMedView ArticleGoogle Scholar
- Bae H, Bowers JH, Tooley PW, Bailey BA: NEP1 orthologs encoding necrosis and ethylene inducing proteins exist as a multigene family in Phytophthora megakarya, causal agent of black pod disease on cacao. Mycol Res. 2005, 109 (12): 1373-1385. 10.1017/S0953756205003941.PubMedView ArticleGoogle Scholar
- Qutob D, Kamoun S, Gijzen M: Expression of a Phytophthora sojae necrosis-inducing protein occurs during transition from biotrophy to necrotrophy. Plant J. 2002, 32 (3): 361-373. 10.1046/j.1365-313X.2002.01439.x.PubMedView ArticleGoogle Scholar
- Veit S, Worle JM, Nurnberger T, Koch W, Seitz HU: A novel protein elicitor (PaNie) from Pythium aphanidermatum induces multiple defense responses in carrot, arabidopsis, and tobacco. Plant Physiol. 2001, 127: 832-841. 10.1104/pp.127.3.832.PubMedPubMed CentralView ArticleGoogle Scholar
- Jores J, Appel B, Lewin A: Cloning and molecular characterization of a unique hemolysin gene of Vibrio pommerensis sp. nov.: development of a DNA probe for the detection of the hemolysin gene and its use in identification of related Vibrio spp. from the Baltic Sea. FEMS Microbiol Lett. 2003, 229 (2): 223-229. 10.1016/S0378-1097(03)00843-7.PubMedView ArticleGoogle Scholar
- Rost B, Yachdav G, Liu J: The PredictProtein Server. Nucleic Acids Res. 2004, 32: 321-326. 10.1093/nar/gkh377.View ArticleGoogle Scholar
- Dunwell JM, Purvis A, Khuri S: Cupins: the most functionally diverse protein superfamily?. Phytochemistry. 2004, 65: 7-17. 10.1016/j.phytochem.2003.08.016.PubMedView ArticleGoogle Scholar
- Straganza GD, Eggera S, Aquinoc G, D'Auriac S, Nidetzkya B: Exploring the cupin-type metal-coordinating signature of acetylacetone dioxygenase Dke1 with site-directed mutagenesis: Catalytic reaction profile and Fe2+ binding stability of Glu-69 Gln mutant. Journal of Molecular Catalysis B, Enzymatic. 2006, 39: 171-178. 10.1016/j.molcatb.2006.01.019.View ArticleGoogle Scholar
- Eulgem T, Rushton PJ, Robatzek S, Somssich IE: The WRKY superfamily of plant transcription factors. Trends Plant Sci. 2000, 5 (5): 1360-1385. 10.1016/S1360-1385(00)01600-9.View ArticleGoogle Scholar
- Robatzek S, Somssich IE: A new member of the Arabidopsis WRKY transcription factor family, AtWRKY6, is associated with both senescence- and defence-related processes. Plant J. 2001, 28 (2): 123-133. 10.1046/j.1365-313X.2001.01131.x.PubMedView ArticleGoogle Scholar
- Qutob D, Kemmerling B, Brunner F, Kufner I, Engelhardt S, Gust AA, Luberacki B, Seitz HU, Stahl D, Rauhut T, Glawischnig E, Schween G, Lacombe B, Watanabe N, Lam E, Schlichting R, Scheel D, Nau K, Dodt G, Hubert D, Gijzen M, Nurnberger T: Phytotoxicity and Innate Immune Responses Induced by Nep1-Like Proteins. Plant Cell. 2006, 18 (12): 3721-3744. 10.1105/tpc.106.044180.PubMedPubMed CentralView ArticleGoogle Scholar
- Bae H, Kim MS, Sicher RC, Bae HJ, Bailey BA: Necrosis- and ethylene-inducing peptide from Fusarium oxysporum induces a complex cascade of transcripts associated with signal transduction and cell death in Arabidopsis. Plant Physiol. 2006, 141: 1056-1067. 10.1104/pp.106.076869.PubMedPubMed CentralView ArticleGoogle Scholar
- Ferre F, Clote P: DiANNA: a web server for disulfide connectivity prediction. Nucleic Acids Res. 2005, 33: 230-232. 10.1093/nar/gki412.View ArticleGoogle Scholar
- Jones AM, Herman EM: KDEL-containing auxin-binding protein is secreted to the plasma membrane and cell wall. Plant Physiol. 1993, 101 (2): 595-606.PubMedPubMed CentralGoogle Scholar
- Woo E, Marshall J, Bauly J, Chen J, Venis M, Napier RM, Pickersgill RW: Crystal structure of auxin-binding protein 1 in complex with auxin. EMBO J. 2002, 21: 2877-2885. 10.1093/emboj/cdf291.PubMedPubMed CentralView ArticleGoogle Scholar
- Cleland RE, Prins HBA, Harper JR, Higinbotham N: Rapid hormone-induced hyperpolarization of the oat coleoptile transmembrane potential. Plant Physiol. 1977, 59: 395-397.PubMedPubMed CentralView ArticleGoogle Scholar
- Warwicker J: Modelling of auxin-binding protein 1 suggests that its C-terminus and auxin could compete for a binding site that incorporates a metal ion and tryptophan residue 44. Planta. 2001, 212: 343-347. 10.1007/s004250000442.PubMedView ArticleGoogle Scholar
- Bauly JM, Sealy IM, Macdonald H, Brearley J, Dröge S, Hillmer S, Robinson DG, Venis MA, Blatt MR, Lazarus CM, Napier RM: Overexpression of auxin-binding protein enhances the sensitivity of guard cells to auxin. Plant Physiol. 2000, 124: 1229-1238. 10.1104/pp.124.3.1229.PubMedPubMed CentralView ArticleGoogle Scholar
- Evans HC: Pleomorphism in Crinipellis perniciosa, causal agent of witches' broom disease of cocoa. Trans Br Mycol Soc. 1980, 74: 515-523.View ArticleGoogle Scholar
- Kilaru A, Bailey BA, Hasenstein KH: Moniliophthora perniciosa produces hormones and alters endogenous auxin and salicylic acid in infected cocoa leaves. FEMS Microbiol Lett. 2007, 274: 238-244. 10.1111/j.1574-6968.2007.00837.x.PubMedView ArticleGoogle Scholar
- Haberlach GT, Budde AD, Sequeira L, Helgeson JP: Modification of disease resistance of tobacco callus tissues by cytokinins. Plant Physiol. 1978, 62: 522-525.PubMedPubMed CentralView ArticleGoogle Scholar
- Krupasager V, Sequeira L: Auxin destruction by Marasmius perniciosus. Am J Bot. 1969, 56: 390-397. 10.2307/2440814.View ArticleGoogle Scholar
- David K, Carnero-Diaz E, Leblanc N, Monestiez M, Grosclaude J, Perrot-Rechenmann C: Conformational dynamics underlie the activity of the auxin-binding protein, Nt-abp1. J Biol Chem. 2001, 276 (37): 34517-34523. 10.1074/jbc.M102783200.PubMedView ArticleGoogle Scholar
- Jennings JC, Apel-Birkhold PC, Mock NM, C J, Baker JDA, Bailey BA: Induction of defense responses in tobacco by the protein Nep1 from Fusarium oxysporum. Plant Sci. 2001, 161: 891-899. 10.1016/S0168-9452(01)00483-6.View ArticleGoogle Scholar
- Simmons CR, Liu Q, Huang Q, Hao Q, Begley TP, Karplus PA, Stipanuk MH: Crystal structure of mammalian cysteine dioxygenase: A novel mononuclear iron center for cysteine thiol oxidation. J Biol Chem. 2006, 281 (27): 18723-18733. 10.1074/jbc.M601555200.PubMedView ArticleGoogle Scholar
- Tian M, Win J, Hoorn R, Knaap E, Kamoun S: A Phytophthora infestans cystatin-like protein targets a novel tomato papain-like apoplastic protease. Plant Physiol. 2007, 143: 364-377. 10.1104/pp.106.090050.PubMedPubMed CentralView ArticleGoogle Scholar
- Woo EJ, Dunwell JM, Goodenough PW, Marvier AC, Pickersgill RW: Germin is a manganese containing homohexamer with oxalate oxidase and superoxide dismutase activities. Nat Struct Biol. 2000, 7 (11): 1036-1040. 10.1038/80954.PubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Type I NLPs alignment. [http://adelmo.cechin.googlepages.com/typeIg0.aln]
- Type II NLPs alignment. [http://adelmo.cechin.googlepages.com/typeIIg0.aln]
- Necrotic sequences alignment. [http://adelmo.cechin.googlepages.com/tInecg0.aln]
- RCSB Protein Data Bank. [http://www.rcsb.org]
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28: 235-242. 10.1093/nar/28.1.235.PubMedPubMed CentralView ArticleGoogle Scholar
- Ceroni A, Passerini A, Vullo A, Frasconi P: DISULFIND: a disulfide bonding state and cysteine connectivity prediction server. Nucleic Acids Res. 2006, 34: 177-181. 10.1093/nar/gkl266.View ArticleGoogle Scholar
- NetNGlyc 1.0 Server. [http://www.cbs.dtu.dk/services/NetNGlyc]
- Guex N, Peitsch M: SWISS-MODEL and the Swiss-PdbViewer: An environment for comparative protein modeling. Electrophoresis. 1997, 18: 14-23. 10.1002/elps.1150181505.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.